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ABSTRACT 


Trade studies of numerous appliance concepts for advanced spacecraft galley, 
p^rlal .ygiene, housekeeping, and other areas were made by t e oe,ng 
Aerospace Company, Contract NAS 9-13965, to determine 

Space Shuttle Orbiter and Modular Space Station mission requirements. In 
conjunction with these studies, analytical models of selected appliance 
conce « not currently included in the C-189A Generalized Environmental/ 
Thermal Control and Life Support Systems (ETCLSS) Computer Program sub- 
routine library were developed. This document describes the new applij e 

subroutines with complete analytical model 

user's input instructions, and validation run results. The 

components modeled were integrated with G-189A ETCLSS models for Shutt e 

Orbiter and Modular Space Station, and results from computer runs of 

systems are presented. 


KEY WORDS 
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1.0 SUMMARY 

Trade studies of numerous spacecraft appliances (Reference 1) have 
identified the most promising appliance concepts for the Space Shuttle 
Orbiter and Modular Space Station. Future spacecraft analysis effort 
will require simulation of these appliances using the G-189A ETCLSS Com- 
puter Program (Reference 2). G-189A subroutines are available for 
simulation of some appliance concepts, while others are so elementary as 
to require no special subroutine. Subroutines which model the remainder 
of optimum appliances chosen from the trade studies have been developed 
and are described in this report. The following six new subroutines 
were written, some of which model more than one appliance component; 

CHILLR - simulates a thermally insulated refrigerator or freezer 
' locker cooled either by an externally chilled fluid or a 

self-contained refrigeration unit 

FTRAY - food warming/serving tray (Skylab-type) 

ROSMOS - reverse osmosis waste water treatment unit 

SHOWER - spacecraft v^hole body shower 

WASDRY - simulates a dishes or clothes washer, dryer, or washer/ 
dryer combination. Will also model a towel /cloth drying 
rack. 

WASTEC - dryjohn (commode) and/or urinal 

Complete user's input instructions for each subroutine are included 
in Section 3, together with the analytical model description and 
method of solution. The order of presentation of data for each 
subroutine corresponds to the format prescribed for new G-189A 
component subroutines in Reference 2. In Section 4 are presented performance 
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1.0 (Continued) 

data for each appliance component model with- constant dummy inlet and 
ambient conditions. These data are compared with test data where avail- 
able, or with analytical predictions, to verify the accuracy of the 
subroutine models. 

The appliance components were incorporated into G-189A models of the 
Space Shuttle Orbiter and Modular Space Station ECLSS to verify their 
performance and operational status within the G-189A software environ- 
ment. For the Shuttle case described in Section 5, an available steady- 
state G-189A Shuttle ECLSS model (Reference 3) was used to simulate four 
orbital phases of operation. Results are presented for the basic unmodi- 
fied Shuttle case, and for the same case except with the appliance models 
added. The computer results demonstrate that the appliance subroutines 
used are compatible with the G-189A program and accurately model the 
respective appliance components. No operational Modular Space Station 
G-189A ECLSS model is currently available, so a simplified model was 
developed and the optimum appliances selected from Reference 1 included. 

For this case a transient 10-hour simulation using typical appliance usage 
schedules, described in Section 6, was run. Complete results are plotted 
for the system as a whole and for each individual appliance component. 

Again, the results demonstrate the usage and accuracy of the new appliance 
subroutines. 

This document, in its entirety, satisfies the requirements of contract 
DRL T-968, Line Item No. 5, "Digital Computer Program Requirements." A 
computer program "Users Manual", described in DRL Line Item No. 3, was 
combined into this single document due to the duplication of a large amount 
of data. The Users Manual consists of three basic parts. The first part — 
mathematical models, solution methods, and users input instructions-- 
comprises Section 3 of this document. The second part— subroutine listings- 
comprises Appendix A. The third part is sample problem descriptions. These 
are given in Section 4, for the individual appliance components, and in 
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1.0 (Continued) 

Sections 5 and 6 for the all-up appliance subsystems within a Shuttle and 
Space Station ECLSS. Additional data input and computer results for the 
Shuttle and Space Station systems are also given in Appendixes D through G. 
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2.0 INTRODUCTION 

Trade studies of numerous spacecraft appliances., (Reference 1) have identified 
the appliance concepts which best satisfy the Space Shuttle Orbiter and 
Modular Space Station mission requirements. Future spacecraft analysis 
effort will require simulation of these appliances using the G-189A ETCLSS 
Computer Program. This program, Reference 2, provides system level ECLSS 
performance simulation by performing mass and energy balances throughout 
all the interactive components and flow loops comprising a total system. 

Use of the G-189A program requires a subroutine for each component in the 
system. These subroutines are all similar in that a standard format of 
G-189A flow and thermodynamic data form the input for each. The input 
data for a given component are taken from the output data from the upstream 
component. The subroutine then must modify these input data in a manner 
which reflects the performance of the component it models, and present the 
output data in the required G-189A format. The program allows the user to 
control or modify the solution as it progresses by calling two subroutines, 
6P0LY1 and GP0LY2, immediately prior to and following each component solu- 
tion. These routines may be used, for example, to alter fluid flow paths, 
turn components on or off, reevaluate component model data based on the 
current solution results, and compute and store parameters for later 
plotting. 

The optimum appliance concepts selected from the trade studies in Reference 1 
are shown in Tables 2-1 and 2-2 for Shuttle Orbiter and Modular Space Station. 
Some of these concepts do not require a new G-189A subroutine since (1) a 
routine is already available, (2) no thermal/mass exchange is involved, or 
(3) operation of the component is so simple it requires, only a minor addi- 
tion to the GPOLY routine logic. Appliances in this category are as follows: 

None needed 
None needed 

Simple GPOLY logic only 
required 


0 Reusable dishes, wet 
and dry wipes 

0 Vomitus collection 

0 Partial body washing, wet 
wi pes 
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2.0 (Continued) 

0 Partial body drying, dry 
wipes or electric dryer 

0 Wet shave 
0 Windup razor 

0 Toothbrush 

0 Vacuum refuse collection 
0 Tape recorder, TV 


•None needed (or a simple 
heater using G-189A routine 
ALT COM) 


GPOLY logic required for 
water usage only 

None needed (or a simple heater 
using G-189A routine ALTCOM if 
electric) 

GPOLY logic required for water 
usage only 

GPOLY logic only required, or 
G-189A routine ALTCOM for an 
electric heater 

GPOLY logic only required, or 
G-189A routine ALTCOM for an 
electric heater 


For the remaining appliances, six new G-189A subroutines have been written, 
some of which will model more than one type of appliance. These subroutine 
have been designated as G-189A component subroutines number 66 through 71, 
and are generally described as follov^s: 


Subroutine 

Number 


Subroutine 

Name 


Description 


71 


66 

69 

67 


CHILLR (simulates a thermally insulated locker 

cooled either by an externally chilled 
fluid or a self-contained refrigeration 
unit) 

* Freezer 

* Refrigerator 

FT RAY 

* Food warming/serving tray (Skylab-type) 

ROSMOS 

* Reverse osmosis waste water treatment unit 

SHOWER 

* Spacecraft whole body shower 
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Subroutine 

Name Description 

WAS DRY 

* Clothes washer 

* Clothes dryer 

* Combined clothes washer/dryer 

* Dishwa^her/dryer 

* Towel/cloth drying rack 

WASTEC , • ' 

* Dryjohn 

* Urinal 

A listing of the subroutines is given in Appendix A. They have been written 
in conventional Fortran V language and are operational on the NASA JSC 
SRU 1108 EXEC II computer system. This document presents complete user's 
input instructions, results of verification runs, and demonstration of their 
operation in all-up Shuttle Orbiter and Modular Space Station ECLSS simulation 
runs. The English system of units is a built-in feature of the G-189A program, 
both for input and output data. Therefore, all input data herein described 
were necessarily presented in. English units to be compatible with the pro- 
gram. Output data are also presented exactly as generated by the program; 
i.e., in English units. Where specific test data are compared with analytical 
results, metric units are also included. 


Subroutine 

Number 


70 


68 
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TABLE 2-1 


SHUHLE APPLIANCE CONCEPTS TRADE STUDY RESULTS 


HABITABILITY 

SUBSYSTEM 

HABITABILITY 

FUtiCTION 

"Appliance 

FUNCTION 

CONCEPT 

CHOSEN 

FIRST 

RATED 

CONCEPT 



FOOD STOfWGE 

REFRIGERATED 

Space Radiator 

Space Radiator 

Thermoelectric 


FOOD 

WARMING 

Heating Trays 

Heating Trays 

Convective Oven 

FOOD 

PREPARATION 





KANAGEMtiNT 








DISH CLEANUP 

Reusable Dishes 

Reusable Dishes 

Reusable Dishes 




and Utensils 

and Utensils 

and Disposable 


GALLEY 


with Disposable 

with Disposable 

Utensils with 


CLEAf^UP 


Wet/Dry Wipes 

Wet/Dry Wipes 

Disposable 






Wet/Dry Wipes 

' 


FECAL 


Apollo System 

Skylab System 



COLLECTION 

Dry John System 




WASTE 

URINE 


Apollo System 

Sky lab System 


COLLECTION 

COLLECTION 






VOMITUS 

Disposable Bags 

Disposable Bags 

Intimate Adaptoi 



COLLECTION 


■ 


PERSONAL 


PARTIAL BODY 

Disposable Wet 

Disposable Wet 

Skylab-Type 

IlYGIENE 

BODY 

WASHING 

Wipe 

Wipe 

Disposable 

Washcloth 


CLEANSING 







PARTIAL BODY 

Disposable Dry 

Disposable Dry 

Electric Dryer 

. 


DRYING 




• * » 


SHAVING 

Safety or 

Safety or 

Safety or 

‘ 

PERSONAL 


Windup 

Windup 

Windup 

j i i 

GROOMING 

DENTAL CARE 

Toothbrush 

Toothbrush 

Electric 




w/Dentifrice 

w/Dentifrice 

Toothbrush 


EQUIPMENT 

SURFACE 

Disposable Wet/ 

Disposable Wet/ 


i • 

CLEANUP 

WIPING 

Dry Wipes 

Dry Wipes 




MANUAL 

Disposable Trash 

Disposable Trash 

Disposable 



COLLECTION 

Bag 

Bag 

Recep tides 


REFUSE 

VACUUM • 

Skylab-Type 

Vacuum-Vented 

Skylab-Type 

HOUSEKEEPING 

MANAGEMENT 

COLLECTION 

Electric 


Electric 



REFUSE 

Storage Bin/ 

Storage Bin/ 

Vacuum Storage 



DISPOSAL 

Container 

Container 



GARMENT/LINEN 

CLOTHES WASH/ 

Disposable 

Disposable 

Mechanical 


MAINTENANCE 

DRY 

Clothes 

Clothes 

w/Clothes Line 



MUSIC 

Cassette Record 

it 

* 




Recorder 




ENTERTAINMENT 

LIBRARY 

Books 

* 

> 

OFF-DUTY 


TELEVISION 

Conmercial Type 

* 

ir 

ACTIVITIES 


GAMES 

Cards, Handball, 


* 




Etc. 




PHYSICAL 

EXERCISERS 

Exer Gym, Hand 

it 

* 


CONDITIONING 


Exerciser 




P NOT TRADED 
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TABLE 2-2 


MODULAR SPACE STATION APPLIANCE CONCEPTS TRADE STUDY RESULTS 


HABITABILITY 

SUBSYSTEM 

HABITABILITY 

FUNCTION 

APPLIANCE 
FUNCTION ' 

CONCEPT 

CHOSEN 

FIRST 

RATED 

CONCEPT 


FOOD 

MANAGEMENT 

FOOD 

STORAGE 

REFRIGERATED 

FROZEN 

Space Radiator 
Space Radiator 

Space Radiator 
Space Radiator 

Thermoelectric 

Tnerme&lectric 

FOOD 

PREPARATION 

WARMING 

Heating Trays 

Heating Trays 

Convective Oven 

GALLEY 

CLEANUP 

DISH CLEANUP 

' Water Spray 
Wash/Elec. 
Heat Dry 

Reusable Dishes 
and Disposable 
Wet/Dry Wipes 

Reusable Cups & 
Dishes - 
Disposable 
Utensils and 
Disposable Wet/ 
Dry Wipes 

PERSONAL 

HYGIENE 

HASTE 

COLLECTION 

FECAL 

COLLECTION 

URINE • 
COLLECTION 

VOMITUS 

COLLECTION 

Ory John System 
Disposable Bags 

Apollo System 
Apollo System 
Disposable Bags 

Sky lab System 

Sky lab System 

Intimate 

Adaptor 

BODY 

CLEANSING 

SHOWER 

PARTIAL BODY 
WASHING 

PARTIAL BODY 
DRYING 

Collapsible 
Reusable Wipes 

Reusable Wipes 

Collapsible 
Reusable Wipes 

Reusable Wipes 

Mechanical 

Skylab-Type 

Disposable 

Washcloths 

Disposable 
Dry Wipes 

PERSONAL 

GROOMING 

SHAVING 

HAIRCUTTING 

NAIL CARE 

DENTAL CARE 

Windup 

Razor Comb 

Vacuum 

Collection 

Manual Clipper 

Toothbrush 
w/Denti f ri ce 

Windup 

Razor Comb 

Vacuum 

Collection 

Manual Clipper 

Toothbrush 

w/Dentifrice 

Vacuum Driven 

Power Clipper 

Vacuum 

Collection 

Nail File 

Vacuum 

Collection 

Electric 

Toothbrush 


EQUIPMENT 

CLEANUP 

SURFACE 

WIPING 

Reusable Wet/ 
Dry Wipes 

Disposable Wet/ 
Dry Wipes 

Sponge 

Skylab-Type 

HOUSEKEEPING , 

REFUSE 

MANAGEMENT 

MANUAL 

COLLECTION 

VACUUM 

COLLECTION 

REFUSE 

PROCESSING 

REFUSE 

DISPOSAL 

Disposable Bags 

Skylab-Type 

(Electric) 

Compactor 
(Air Pressure) 

Storage Bin/ 
Container 

Disposable Bags 

Skylab-Type 

(Electric) 

Compactor 
(Air Pressure) 

Storage Bin/ 
Container 

Disposable 

Recepticles 

Vacuum Vented 

Compactor 

(Vacuum) 

Vacuum Storage 


GARMENT/LIMEN 

MAINTENANCE 

CLOTHES 

WASH/ORY 

Water Spray 
Agitation Plus 
Electric Dry 

Disposable 

Clothes 

Water Spray 
Agitation Plus 
Clothes Line 

OFF-DUTY 

ACTIVITIES 

ENTERTAINMENT 

— 

MUSIC 

LIBRARY 

TELEVISION 

GAMES 

- 

Casette Recorder 
Books 

Commercial Type 

Cards, Handball; 
Etc. 

it 

it 

* 

■* 

* 

* 

* 

* 


PHYSICAL 

CONDITICHIN'G 


EXERCISERS 


Exer Gym, 

Hand Exerciser 
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3.0 SUBROUTINE OPERATION AND WTH MODEL DESCRIPTIONS 

In this section are described (1) all user’s input instructions required 
to operate the subroutines, (2) the analytical math model assumed to des- 
cribe the component, and (3) the methods and logic used to obtain the 
solution. This section may serve independently as a user's manual for 
the subroutines. The symbols used in the subroutine descriptions are 
defined at the end of each section. The Fortran V listings of the sub- 
routines are given in Appendix A. The logic flow chart for each subroutine, 
and the definition of the Fortran names used in the subroutine, are also 
given at the end of each section. The description of each subroutine is in 
the format prescribed in Reference 2 for new G-189A subroutine writeups as: 


1.0 Subroutine Description 

2.0 Subroutine Data 

2.1 General Notes 

2.2 Instruction Options 

2.3 Heat Loss V-Array Data 

2.4 Steady-State K-Array Data 

2.5 Steady-State V-Array Data 

2.6 Transient K-Array Data 

2.7 Transient V-Array Data 

2.8 Extra K-Array Data 

2.9 Extra V-Array Data 

3.0 Analytical Model Description 




Subsections not 
required for a 
particular subroutine 
are omitted. 




It is assumed that the user of these subroutines is familiar with the basic 
operation of the G-189A program. Every G-189A component has two possible 
fluid flow paths through it, as shown below. One side is arbitrarily 



Primary 

fluid 

outlet 


Secondary 

fluid 

outlet 
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3.0 (Continued) 

designated as primary and the other secondary. The description of each 
component must specify which flow paths are used. For example, for the 
SHOWER component, the primary inlet is used for the air flow, the secondary 
inlet is water, and the air/water mixture is outlet on the primary side. 

The secondary outlet for the SHOWER is not used. 


All the data for every G-189A component are stored in a V-array (which is 
equivalenced to a K-array to allow integer data). During solution of a 
particular component, its individual floating point data are loaded from 
the V-array into an R-array. The component subroutine then operates on 
these data in a manner which reflects the performance of that component. 
Thus, the data input descriptions in the following sections are concerned 
primarily with the definition of each K and R-array location for each 
component. Some of these data locations are used to store input data, and 
others output data. Thus the following convention has been followed to 
define the type of data for each component reference location: 


Data Type 


Meaning 


I(R) Input data - required 

1(0) Input data - optional 

0 Output (computed) data 

I(R), 0 Initial input data required; thereafter 

data is computed and output 

Default values for the model input data are built into each subroutine to 
describe specific appliances selected in the trade studies of Reference 1. 

For the refrigerator/freezer subroutine, the default values- describe the 
Shuttle freezer kit design described in Reference 32. The user may there- 
fore simulate these specific appliances with a minimum of input; or he may 
simulate other related appliances with the same subroutines by supplying 
new, input data, which would then override the default values. Data reference 
locations for which default input data are included in the subroutine are 
designed by 1(0) in the writeups. The default input values are listed in a 
table at the end of the model description for each subroutine. 
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3,0 (Continued) 

A| major part of each appliance subroutine is the thermal model used to 
. simulate the heat transfer within the appliance and between the appliance 
and its surroundings. An equivalent electrical resistor/capacitor nodal 
network was used in each case for this purpose. These nodal models are 

shown for the various appliances in the following sections using these 
symbols; 

• Node with thermal mass 

O Steady-state node 

M Boundary node 

^WV— Thermal conduction or 

convection conductor (linear) 

— Thermal radiation conductor (nonlinear) 

— ► One-directional fluid flow 
conductor 

>■ Other heat addition 

•HI— Thermal capacitance 

•'l|— — Thermal ground 

Nodes having thermal mass are solved by equating the net heat input to the 

change in heat storage. A thermal capacitance is specified for these nodes 
as defined by the relation 

' ""=p 

c = nodal thermal capacitance 

m = nodal mass 

Cp = specific heat 

Steady-state nodes are used to model gaseous fluid or other special nodes 

having negligible thermal mass. These nodes are in thermal equilibrium 
with their surroundings; that is, their temperature is computed such that 
the heat in is equal to the heat out. Boundary node temperatures are not 
computed in the subroutines. They must be input by the user and may be 
held constant or varied during a run based on the progressing solution. ■ 


3-3 



D2-118571-1 


3.0 (Continued) 

The nodes are interconnected by thermal conductors evaluated in the following 
ways : 




solid conduction 


- 

fluid convection 

“p 

- 

fluid flow 

oAp- 

— 

radiation 


where G = thermal conductor 

k. = material thermal conductivity 
A = "window" area between nodes 
Ag = surface area of node 

A = length between nodal centers 
m = fluid mass flow rate 

= convection heat transfer coefficient 

a = Stefan-Boltzmann radiation constant 
'3- - radiation interchange factor 

The first three conductors are referred to as linear and transfer heat 
proportional to the first power temperature difference (T. - T^. )• The 
fourth conductor is radiation whi.ch transfers heat as a function of the 
fourth power temperature difference (T> - T^.^). Some conductors are 
designated as "one-directional" elements, meaning that heat is transferred 
through them in one direction only. This feature is typically used for 
fluid flow simulation, in which the stored energy travels downstream only, 
and also used for satisfying certain boundary conditions at a line of 
symmetry within a model. 

The G-189A program is designed to run either with or without a system 
pressure drop analysis. If this option is not used, the system pressures 
will remain constant around each flow loop, except where they are set by 
the user when necessary by 6P0LY logic. For problems requiring pressure 
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3.0 (Continued) 

drop anailysis, a set. of six standard pressure drop equations is provided 
within the G-189A program to model any component. These equations are 
described in the G-189A program manual, Reference 2, and are shown in 
Table 3-1 taken from that reference. They are considered adequate to 
describe the pressure drop characteristics of the new appliance subroutines. 
The methods and recommended pressure drop data are included in the following 
sections for each appliance subroutine. 
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TABLE 3-1 

G-189A PRESSURE DROP OPTIONS 


AP Option Code 


Pressure Drop Equations Solved 


K Array Data Req'd 
(PDK* Loads heet) 


0 (default) 


AP = 0 


hf L 


* AP = 


c 


2o g. 


(psf) 


where: 


<jJ 

Jri 


t 


V 

p 


= Cross sectional area (ft ) 

= Hydraulic diameter (ft) 

= Fanning friction factor (dimensionless) 

= Newton* s lav conversion factor 

(Ih - ft/lb_ - hr2) 
m r 

= Equivalent length of pipe or duct for 
. pressure drop calculations (ft) 

a Flow (lb /hr) 

3 

= Flow stream density (ib/ft ) 


The fanning friction factor, f, is calculated 
using one of the following equations which 
are selected on the basis of the fluid stream 


Reynolds number. Re 


P A. 


V Array Data Req'd 
(PDPR or PDSE Loadsheet) 



If A is not input the 
flov'^passage is assumed 
to be circular and A 
is calculated as 
follows : 
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TABLE 3-1 (Continued) 

G-189A PRESSURE DROP OPTIONS 

K Array Data Req’d V Array Data Req’d 

Pressure Drop Equations Solved (PDK* Loadsheet) (PDPR or PDSE Loadsheet 


f 


16 

Re 

0.0791 

Re 


(0<Re<2100) 

(2100<Re<100000) 


» O.OOlU + ' (l00000<Re<-) 

. Re 

*N0TE: If the fluid is gaseous and a Mach 

number >0.1 is calculated then AP 
option codes of 1 or 2 will use 
compressible flow equations and 
solve the pressure drop using a 
method similar to that described 
by Shapiro in Reference 5-1* 


» AP 



(psf) 


where: 

K = Loss coefficient (dimensionless) 

A^, g^, w, p «* As defined under AP 
option code 1 


K. 
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‘'Oils POOR 




AP Option Code 


TABLE 3-1 (Continued) 

G-189A PRESSURE DROP OPTIjONS 

Pressure Drop Equation Solved 


K Array Data Req'd 
(PDK* Loadsheet) 


AP » C 


(psf) 


CO 

I 

cx> 


vhsre: 


c = Constant coefficient (dimensionless) 

n = Exponential coefficient for flow 
(dimensionless) 

v,p “ As defined Tinder AP option code 1 




AP = AP (psf) 
where : 

= Conversion factor 

5*20U (psf/in. H 2 O) for gases 

= ll4lt.O (psf/psi) for liquids 


_ Interpolated value for oAP 

n ^ 


where : 

Cj^.D « 


CiP 


As defined under AP option 
codes 1 and k 


Input table number for 
pAP vs. w data 


Fo r gases : 

Dependent variable; 
pAP (lb - in. HgO/ft’) 

1st Independent variable: 
V (Ib/hr) 


V Array Data Req'd 
(PDPR or PDSE Loadsheet) 


C y n 


AP 

Input AP as (In. H_0) 
for gas flows or (psi) 
for liquid flows. 
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TABLE 3-1 (concluded) 

G-189A PRESSURE DROP OPTIONS 

K Array Data Req'd 

A? Option Code Pressure Drop Equation Solved (PDK* Loadsheet) 


For liquids : 

Dependent variable: 
pAP (Ib-psi/ft^) 

1st Independent variable: 
V (ib/hr) 


Input table number for 

AP vs w data 


For gases : 

Dependent variable: 

AP (in. HgO) 

1st Independent variable: 
w (ib/hr) 

For liquids : 

Dependent variable: 

AP (psi) 

let Independent variable: 
w ( Ib/hr) 


6 AP = (Interpolated value .for AP) 

where: 

C. . = As defined under AP option 

code L 


V Array Data Req'd 
(PDPR or PDSE Loadsheet) 
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3.1 CHILLR 

The CHILLR subroutine is designated as G-189A No. 71. 

3.1.1 Subroutine Description 

The CHILLR subroutine vVill simulate a refrigerator or freezer cooled either 
by an externally chilled coolant (e.g., from a radiator coolant circuit) 
or by a self-contained refrigeration unit. The model is fairly generalized 
and may be used to simulate the following configurations: 


Refrigerant 


out 


Refrigerant 
in 


^ Fan (optional) 


(a) Cooled by externally 
chilled liquid 


.Cooling 

Coils 




Refrigeration 
- Unit /■ 


Chiller 
Coolant 
in 


(b) Self-contained refrigeration 
unit with wrapped cooling coils 



.(c) Self-contained refrigeration 
unit with cooling coil interface 


(d) Self-contained refrigeration 
unit with cold plate interface 


Figure 3-1. CHILLR Component and Flow Schematic 
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3.1.1 (Continued) 

The thermal model in either case is shown in Figure 3-2. In addition to 
the thermal network shown in the figure, the locker inner walls and con- 
tents are thermally connected to the ambient surroundings using the standard 
G-189A subroutine QSURR. This routine models the heat exchange from an 
arbitrary structure to ambient via insulation, thermal shorts, and conduction/ 
convection/radiation paths. The output from the QSURR subroutine defines 
the heat loss from the internal structure to ambient, which is desig- 
nated as in Figure 3-2. 

The CHILLR subroutine has been used to simulate the Shuttle food and 
medical sample freezer kit described in Reference 32. Excellent correlation 
between the subroutine results and independent detailed freezer thermal 
analysis has been obtained, as presented in Section 4.1. The model input data 
used for that freezer design are included directly in the subroutine as 
default data, and are listed in Table 3-2. 

3.1.2 Subroutine Data 

3. 1.2.1 General Notes 

a. Only the component primary side is used. If an externally chilled 
coolant is assumed, the coolant must inlet and exit on the primary 
side. If a self-contained cooling unit is used, the primary fluid is 
assumed to cool the heat rejection device at the hot side of the unit. 

b. Any flow codes may be used, but no fluid phase change is allowed. 

c. If ambient air properties (for evaluating door opening effects) are 
taken from an associated cabin component [an option with instruction 
NSTR(4)], then the cabin component number must be included on the 
KBAS card for the refrigerator or freezer. The CHILLR routine uses 
cabin V-array locations 1 through 19. 
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(a) Thermal model of refrigerator/freezer 
locker including cooling coils 


Cooling 

coils 




(b) Thermal model of refrigerator/ 

• freezer chiller side, including 
cooling coils - assuming 
externally chilled coolant fluid 


(c) Thermal model of refrigerator/ 
freezer chiller side, including 
cooling coils - assuming 
self-contained cooling unit 


RFPTiODUCIBiLiTY OF THE 
IS POOR 


Figure 3-2. Thermal Model of Refrigerator/Freezer 
Locker (a) and Chiller (b) and (c) 
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3. 1.2.1 (Continued) 

d. The net heat loss from the refrigerator or freezer locker, locker fan, 
and door opening effects to the ambient surroundings should be added 
to the associated cabin sensible heat input. This heat loss is given 
in R{53), and should be added to cabin R-array location 66 in 6P0LY 
logic. 

3. 1.2.2 Instruction Options 

NSTR(l): Flag used to designate type of refrigerator or freezer being used 

=0 No self-contained refrigeration unit is used. Unit is cooled 
by chilled liquid inlet to primary side from external coolant 
circuit. 

=1 A self-contained refrigeration unit is used, with the heat 
rejection device cooled by the primary side flow. 

NSTR(2): Control logic for fan inside refrigerator or freezer locker 

=0 No fan is used. 

= 1 Fan is on continuously if transient case, or for fraction of 

time in R(93) if steady state. 

= 2 Fan is controlled ON/OFF with the self-contained cooling unit 
according to NSTR(3) (transient case with NSTR(l) = 1 only]. 

NSTR(3): Temperature control logic (used in transient case only). If 

a self-contained cooling unit is used, this will control the 

cooling unit ON/OFF cycle; and also the locker fan ON/OFF cycle’ 
if NSTR(2) = 2. If an external coolant circuit is used, this 
will control the fan ON/OFF cycle if NSTR(2) = 2. 
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3. 1.2.2 (Continued) 

= 0 Previous components turned ON/OFF based on user's input control 
parameter in R(113). 

= 1 Previous components turned ON/OFF based on locker internal 
temperature .limits between R(83) (low temperature to turn 
components OFF) and R(84) (high temperature for components to 

turn ON). 

= 2 Previous components cycled ON/OFF according to a timer cycle 
in R(llO) through R(112). 

NSTR(4): Identifies ambient gas properties for door opening effects 

=’0 Use cabin component NCAB (in K-array location 8) 

_ j Us0 gas properties input in R(96) through R(99) 

NSTR(5): Door opening schedule 

= 0 Input average number of times door is opened per day [R(71)] 
and mission time for first opening in R(114) 

= 1 Input arbitrary daily door opening schedule in extra V-array 
location beginning with R(125) 

NSTR(6): Flag used to identify type of refrigeration unit (used only 

if NSTR(l) = 1) 

= 0 Stirling cycle 

= 1 Vapor compression 
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3. 1.2.3 Heat Loss V-Array Data 


Reference 

Location Description 

51 Temperature of locker internal walls 
and food outer surface (®F) 

52 Summed thermal conductance from locker 
internal walls to ambient surround- 
ings, excluding conductance through 
locker fan (Btu/hr °F) 

53 Net heat loss from locker internal sur- 
face, from fan (if present), and from 
door opening to ambient surroundings 
(Btu/hr) 

54 Ambient gas temperature (°F) 

55 Thermal conductor from refrigerator 

or freezer locker external surface 
(not including fan) to ambient gas 
(Btu/hr °F) 

56 Convective heat loss from locker outer 
surface (excluding fan) to ambient 

t gas (Btu/hr) 

57 Ambient wall temperature (°F) 

58 Thermal radiation' conductor (A3) between 

locker external surface and ambient 
walls (sq ft) 

59 Radiative heat loss from refrigerator 
or freezer locker to 'ambient walls, 
excluding fan (Btu/hr) 

60 Temperature of ambient structure 
attached to locker (°F) 

61 Thermal conductor through structural 
hard-attach-points between locker 

• inner walls and attached ambient struc- 

ture, excluding conductance through 
insulation (Btu/hr °F) 

62 Thermal loss through structural hard- 
ly attach-points (Btu/hr) 

: - 

f 


Data Type 
1 ( 0 ), 0 

0 

0 

1 ( 0 ) 

1 ( 0 ) 

0 

1 ( 0 ) 

1(0) 

0 

1 ( 0 ) 

1(0) 

0 
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3. 1.2. 3 (Continued) 


Reference 

Location 

Description 

Data Type 

63 

Refrigerator or freezer locker external 
surface temperature (°F) 

0 

64 

Thermal conductor through locker 
wall insulation (Btu/hr °F) 

1(0) 

3.1, 2. 4 

Steady-State and Transient V-Array Data 



Note: Locations used for transient runs only are marked with 


an asterisk. 


Reference . 

Location Description 

Data Type 

65 

Net cooling provided to refrigerator or 
freezer locker (Btu/hr) 

0 

66 

Set point temperature; used to 
determine steady-state duty cycle for 

1(0) if steady state, 
NSTR(l) = 1 and 
NSTR(6) = 0 


Stirling cycle cooling unit ( F) 

67 

Steady state duty cycle used for _ 
self-contained cooling unit (fraction 

1(0) if NSTR(l) = 1 
0. if NSTR(6) = 0 


of time on) 


68 

Total internal locker volume including 
food (cu ft) 

1(0) 

69 

Packaged food volume (cu ft) 

1(0) 

70 

Air change per door opening 
(fraction of void volume) 

1(0) 

71 

Number of times per day door is to be 
opened 

1(0) if steady state 
or NSTR(5) = 0 

*72 

Dry food mass (lbs) 

1(0) 

*73 

Fraction of total food mass assigned 
to food outer surface node 

1(0) 

*74 

Refrigerator or freezer locker inner 
shell thermal mass (Btu/°F) 

1(0) 

75 

Food inner node temperature (°F) 

1(0), 0 
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3. 1.2. 4 (Continued) 
Reference 


Location 

“Description 

Data Type 

76 

Ratio of frozen food thermal 
conductivity to unfrozen food 
thermal conductivity 

1(0) for freezer only 

77 

Air temperature inside food 
compartment (°F) 

I(R) if NSTR(2)?«0, 0 

78 

Thermal conductor, food inner node 
to outer node, in unfrozen condition 
(Btu/hr °F) 

1(0) 

79 

Air temperature inside chilling 
compartmen,t.(°F) - 

I(R) if NSTR(2)j«0, 
0 

80 

Surface area of cooling coils (sq ft) 

I(R) if NSTR(2)/0 

81 

Surface area of food and inner locker 
walls (sq ft) 

I(R) if NSTR(2)/0 

82 

Conductive thermal conductor from 
cooling coils to surface of food 
and inner locker walls (Btu/hr °F) 

1(0) 

*83 

Food surface lower temperature limit 
for cooling unit and/or fan to turn 
off(°F) 

1(0) if NSTR(3)=1 

*84 

Food surface high temperature limit 
for cooling unit and/or fan to turn 
on (°F) 

1(0) if NSTR(3)=1 

85 

Convective heat transfer coefficient; 
air inside locker to food and walls 
inner surface and cooling coils 
(Btu/hr-sq ft-°F) 

I(R) if NSTR(2)/0 

86 

Conductive thermal conductor from fan 
to attached ambient structure 
(Btu/hr °F) 

I(R) if NSTR(2)/0 

87 

Effective fan housing temperature (°F) 

I(R) if NSTR(2)/0, 
0 

88 

Thermal radiation conductor A3- from 
fan housing to ambient walls (sq ft) 

I(R) if NSTR(2)/0 
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3. 1.2.4 (Continued) 
Reference 


Locati on 

Descriotion 

89 

Conductive thermal conductor from fan 
housing to inner locker walls and 
food surface (Btu^r °F) ° 

90 

hSnl'ln from fan 

housing to ambient gas (Btu/hr °F) 

91 

* 4'* 

(watt°r^ electrical power 

92 ‘ 

Fan motor efficiency (fraction) 

93 

Fraction of time fan is on duririq 
steady-state run ^ 

*94 • 

Combined thermal mass of food outer 
(Btu/°F)"°^^ 

*95 

Food inner node thermal mass (Btu/°F) 

96 

Inlet gas temperature when door is 
opened (°F) 

97 

Ob/L^ft)^^^^^'^^ opened 


98 

99 

100 

101 

102 

*103 

104 


.Inlet gas absolute humidity when door 
IS opened (lb water vapor/ lb dry air) 

Inlet gas specific heat when door is 
opened (Btu/lb °F) • 

Coefficient of performance of 

cooling pro- 

silould' not include 

coolant pump power 

Air circulation flow rate provided 
by locker fan (cfm) 

Constant evaporator temperature for 
self-contained refrigeration device (°F 

Cooling coils thermal mass-dry (Btu/°F) 

Cooling coils temperature (°F) 


Data Type 
I(R) if NSTR(2)/0 

I(R) if NSTR(2)/0 
I(R) if NSTR(2)/0 

I(R) if NSTR(2)/0 

I(R) for steady 
state 

0 

0 

1(0) if NSTR(4)=1 
1.(0) if NSTR(4)=1 
1(0) if NSTR(4)=1 
1(0) if NSTR(4)=1 
1(0) if NSTR(l)=l 

I(R) if NSTR(2)/0 

I(R) if NSTR(1)=1 
and NSTR(6)=1 

1 ( 0 ) 

1 ( 0 ), 0 
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3. 1.2. 4 (Continued) 


Reference 


Location 

Description 

Data Type 

105 

Heat added into locker by opening 
door on previous iteration (Btu/hr) 

0 

106 

Total cooling capacity of Stirling 
cycle cooling unit (Btu/hr) 

1(0) if NSTR(1)=1, 
NSTR(6)=0 

107 

Coolant pump power used with Stirling 
cycle cooling unit (watts) 

1(0) if NSTR(1)=1, 
NSTR(6)=0 

108 

Fraction of coolant pump power 
transferred to coolant with 
Stirling cycle unit 

1 ( 0 ) if NSTR(1)=1, 
NSTR(6)=0 

109 

Total electrical power used on 
previous iteration by self-contained 
cooling unit (watts) 

0 

*110 

Preset time for cooling unit and/or 
fan to remain on with timer (minutes) 

I(R) if NSTR(3)=2 

*111 

Preset time for cooling unit and/or fan 
to remain off with timer (minutes) 

I(R) if NSTR(3)=2 

*112 

Initial displacement of time into timer 
ON/OFF cycle (minutes) Note: The. cycle 
begins with the ON setting. 

I(R) if NSTR(3)=2 

*113 

User's control parameter to turn on 
cooling unit and/or fan. Set to 1.0 
v;hen these components are to be turned 
on, or 0.0 when they. are to be off. 

I(R) if NSTR(3)=0 

*114 

Next mission time which locker door 
shall be opened (seconds) 

I(R) if NSTR(5)=0, 
0 

*115 

Total water content (frozen and 
unfrozen) in outer food node (lb 
water/lb total food mass in outer 
food node) 

1 ( 0 ) 

*116 

Total water content (frozen and 
unfrozen) in inner food node (lb 
water/lb total food mass in inner 
food node) 

1(0) 
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3. 1.2. 4 (Continued) 
Reference 


Location 

Description 

Data Type 

*117 

Fraction of water in outer food node 
that is frozen 

I(R) only if food 
temperature initial 
ized at 32.0°F, 0 

*118 

Fraction of water in inner food node 
that is frozen 

I(R) only if food 
temperature initial 
ized at 32.0°F, 0 

*119 

Specific heat of dry food (Btu/lb °F) 

1(0) 

120 

Conductive thermal conductor from 
cooling coils to self-contained 
cooling unit structure (Btu/hr °F) 

1(0) 

121 

Cooling unit structure temperature 
attached to conductor R(120) (°F) 

1(0) 

122-124 

Not used 
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3. 1.2.5 Transient K- Array Data 


Reference 

Location Description Data Type 

16 R-array location to be used for next 1(0) if NSTR(5)=1, 

mission time locker door will be opened; 0 
equal to 125 or greater.. 


17 


This location stores cooling unit 
mode on previous iteration: 0=0FF, 
1=0N. Note: This location used to 

specify initial ON/OFF mode for 
cooling unit. 


I(R) if NSTR(3)=1, 
0 


18 Not used 

19 Number of mission days for which _ 0 

■ solution has been obtained; used in 

door opening schedule logic 


3.1. 2. 6 Extra V-Array Data (if NSTR(5)=1, there must be reserved N+1 

locations where N is the number of locker door openings per day.) 


Reference 

Location Description Data Type 


125— >124+N In these locations are stored the 
mission times (hrs) when the locker 
door is to be opened, where N is 
the number of door openings per day_. 
Each opening is automatically 
repeated every 24 hr$. 

125+N Must be set to zero to indicate the 

end of door opening time entries 


I(R) if NSTR(5)=1, 
0 


I(R) if NSTR(5)=1 


3.1.3 Analytical Model Description 
3. 1.3.1 Fan Thermal Balance 

Use of a fan to circulate the locker air is optional according to the 
input instruction NSTR(2). If the fan is omitted, all the thermal con- 
ductors connected to it are simply ignored in the model. If a fan is 
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3 . 1.3.1 (Continued) 

included, its electrical input power (qp) and motor efficiency (n) must 
be input. The heat generated within the fan motor is then given by 

djuj," (1 " n) (3.1.1) 

(See Table 3-1 for definition of symbols.) The output shaft power is 
assumed to be dissipated as heat of compression of the inside locker air 
according to the relation: 


qa = nqp 


(3.1.2) 


The fan motor/housing is assumed to be a steady-state node; that is, its 
thermal mass is neglected and it is in thermal equilibrium with its sur- 
roundings (heat in = heat out). An energy balance on the fan is therefore 
as follows: 

qp = q^ + Gj^(Tp-T^) + ^ ^'^F~'^2^ * '^*^3 (3.1.3) 

v;here the temperatures in the radiation term are in '*R. To solve 
equation (3.1.3) for fan temperature Tp, the radiation term is first 
linearized as follows: 

aG3 (Tp^-T3^) = crGg (Tp^+Tg^) {TP+T3) (TP-T3) (3.1.4) 


The effective radiation conductor Y is assumed to remain constant over 
a single time step (or steady-state iteration). Thus, it may be evaluated 
first, and the fan temperature obtained by solving equations (3.1.3) and 
(3.1.4): 


= *^F " ^a Vo ^ ^~^3 

F ■ G^ + G^ + G^ + Y 


(3.1.B) 
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3. 1.3.1 (Continued) 

I"f the fan is turned off, or if there is no fan present, the convection 
conductors G^, G^, and G^ in Figure 3-2 are assumed zero. The air 
temperatures in the food and chilling compartments are set to the food 
temperature and cooling coils temperature, respectively. When the fan is 
operating, these temperatures are computed by equating the heat into and 
out of the air nodes, as follows: 


(v\) = S, (3.1.6) 

(V\) (V'c) (3-1-7) 

Solving equations (3.1.6) and (3.1.7) for the locker air temperatures, 
one obtains: 



(G^+GJ + G^GJ^ 




G 6 
a 0 


(3.1.8) 


+ G 

0 0 a a 


G* + G 
0 a 


(3.1.9) 


3. 1.3.2 Storage Compartment Thermal Solution 

The steady-state nodes in the model (assumed to have negligible thermal 
mass) are all connected with the fan, as seen in Figure 3-2, and their 
solution was discussed in Paragraph 3. 1.3.1, The remaining nodes normally 
have thermal mass, although it is neglected when obtaining a steady-state 
solution. Their solution is obtained by a standard forward difference 
method as follows: 

"'"new ^old t ^ ^into node (3.1.10) 

(One exception to this equation is the cooling coils temperature T when a 
vapor compression unit is operating, as discussed in Section 3. 1.3. 4). For 
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3. 1.3.2 (Continued) 

a transient solution) the value of y is given by 

^transient ” mCp (3.1.ir) 

The value of mCp is simply the thermal mass of the node. The computing 
time increment At is taken to be 0.4 times the smallest of the time incre- 
ments required for computational stability, where 

= thermal mass . . 

stability sum of conductors to ground v3.l.l2; 

If the stability time increment is larger than the G-189A system model time 
step, then the system time step is used in equation (3.1.11). If the sta- 
bility time increment is smaller than the system time step, then internal 
iterations are performed until the required system time step is achieved. 

For the steady-state solution, the value used for y in equation (3.1.10) 
is given by 

^steady state ~ sum of conductors to ground (3.1.13) 

For steady-state solutions, internal iterations of these nodes using 
equations (3,1.10) and (3.1.13) are continued until the old and new temper- 
atures differ by less than 0.001°F, or until 25 iterations are made. 

3. 1.3.3 Freezing Effects 

The thermal properties of the food nodes are assumed to vary with, the 
water/ice balance. When food temperatures drop to 32° F, their temperature 
is held constant until the latent heat of their moisture is removed. Thawing 
of frozen food was considered to represent a failure condition and 
was not included in the model. The thermal conductivity of unfrozen food 
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^ 3. 1.3. 3 (Continued) 


is assumed until all the liquid water has frozen, after which 
. of frozen food is used. The thermaT mass of the food nodes 
ient case) is computed as follows: 

the conductivit;y 
(in the trans- 

‘"’^p’food ‘ ‘“p’bone-dr/ 

food 

(3.1.14) 

•* 

3. 1.3.4 Self-Contained Refrigeration Unit 



If the input instruction NSTR(l) is set to 1, a self-contained refrigeration 
unit is assumed as follows: 



Logic for two types of refrigeration units is included in the subroutine, 
depending on the input value of NSTR(6). If NSTR(6)=0, a Stirling cycle 
is assumed. This type of unit has been selected for the conceptual design 
of a Shuttle food and medical sample freezer, described in Reference 32. 

If NSTR(6)=1, a vapor compression unit is assumed, with a compression/ 
condensation/expansion/evaporation cycle. For both cases, the coefficient 
of performance of the unit, defined as follows, must be input: 

COP 5— (3.1.15) 
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3. 1.3. 4 (Continued) 


Where is the total cooling provided to the locker, and is the 
electrical power input to the unit. For the vapor compression unit, the 
coolant pump, if present, is assumed to cycle on and off with the refrigera- 
tion unit, and the electrical power q^ should include the pump power. To 
be compatible with the Shuttle freezer design, the coolant pump used with 
the Stirling cycle is assumed to be on continuously. It is therefore 
treated separately from the Stirling unit, and the pump electrical power 
should not be included in q^ . The heat rejected at the fjreezer hot side 

is given by 

qj^ = q|_ + q^ + (T— j)lp (3.1.16) 


where qp is the coolant pump power for the Stirling cycle (or zero for 
vapor compress! on) , and j is the fraction of the pump power transferred 
to the coolant. Combining equations (3.1.15) and (3.1.16), the total heat 
rejected by the cooling unit to the environment is given by 

= q^ (1 + COP) + (1 - j)qp = q|_0 + ^) + (1 - j)qp (3.1.17) 

This rejected heat is all put into the primary flow which is assumed to 
be directed to the chiller unit heat rejection device. Thus, the primary 
flow outlet temperature is given by 


= T . 

mary out 'primary in 


qp (1 + COP) + (1 - j)qp1 
(mcD 


/I 1 i«i 


For the vapor compression cycle, the cooling coils temperature T^ is assumed 
to remain constant when the unit is turned on, based on the design operating 
pressure of the working fluid. For the Stirling Cj'Cie, the rate of heat 
removed from the cooling coils is assumed to be constant when the unit is on. 

For steady state runs using a self-contained refrigeration unit, an 
equipment duty cycle is input to define the fraction of time the unit is 
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3. 1.3. 4 (Continued) 

turned on. The total cooling provided, , and electrical power, , 
are multiplied by this fraction before using in equations (3.1.T6) through 
3.1.18). For the Stirling cycle case, the duty cycle is internally adjusted 
to give the required control temperature input in R(66). These adjustments 
are made using the following relation: 

T - T 

ambient control 

^ambient “ ^food 



3. 1.3. 5 Locker Door Opening Effects 

The heat input from air exchange through the locker door will depend on 
the temperature, density, specific heat and humidity of the ambient air. 

These properties may be input directly, or a cabin component may be used from 
which to obtain the data, as specified in input instruction NSTR(4). The 
door may be opened according to an arbitrary input schedule (transient 
runs only) or at constant periodic intervals, according to input instruction 
NSTR(5). The heat input from the door is given by 

a = D V, . ffc (T .. . - T ) + aQ, .x]/At (3.1.19) 

'^door *^3 locker ^ p,' ambient ;o' ^LAT-" 

a 

where f is the fraction of locker void volume assumed to be exchanged 
with ambient air. The time At is the system computing time step for 
transient runs, and the average time ‘between door openings for steady-state 
runs. 

3. 1.3. 6 Pressure Drop Considerations 

If a self-contained cooling unit is used, the only fluid flow connections 
within the ECLSS model will be to provide a fluid to cool the heat rejection 
device. Typically, this will involve either cabin air flow or a system 


^new ~ *^old 
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3.1 .3.6 (Continued) 

coolant flowing through some type of heat exchanger. In the former case, 
no pressure drop model data would be required for the CHILLR component 
since it would not be connected within an actual fluid loop. For the 
latter case, the specific pressure drop data to be used would depend on 
the geometry of the heat exchanger device. This has not yet been designed 
for a spacecraft refrigerator or freezer, and therefore the pressure drop 
model data cannot now be specified. Similarly, if an externally chilled 
fluid were used to cool the refrigerator or freezer, the geometry of the 
cooling coils inside the locker would have to be known in order to specify 
pressure drop model data. However » it is believed in both cases that prob 
ably the G-189A pressure drop option (1) in Table 3-1 would best describe 
the device. Only the flow path cross-sectional area, hydraulic diameter, 
and effective length would be input. The program would then determine the 
flow regime (laminar, transition, or turbulent) and compute the pressure 
drop accordingly. 
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TABLE 3-2 


DEFAULT VALUES FOR CHILLR SUBROUTINE INPUT DATA 


R- Array* 

■ • • 

Default 

R-Array 

Default 

Location 

Value 

Location 

Value 

51 

- 9.5 

84 

- 7 

54 

70 

961 

70 

55 

33.5 

97 

0.075 

57 

70 

98 

0.0078 

58 

11 

99 

0.24 

60 

70 

100 

0.364 

61 

0.2 

103 

3.8 

64 

1.48 

104 

A(2) NSTR(1)=0 

66 

- 8.5 


-10 NSTR(1)=1, NSTR(6)=0 

67 

0.63 


R(102) NSTR(1)=1 , NSTR(6)=1 

68 

4.63 

106 

75 X 3.412 

69 

4.13 

107 

5 

70 

1 

108 

1 

71 

18 

115 

0.5 

72 

135 

116 

0.5 

73 

0.17 

117 

0 R(51) > 32 

74 

3.8 


1 R(51) < 32 

75 

- 8.6 

118 

0 R(75) > 32 

76 

3.5 


1 R(75) < 32 

78 

18.57 

119 

0.49 

82 

260 

120 

0.129 

83 

- 10.5 

121 

97 


K(16)=125 


*See Section 3.1.2 for definition of variables. 
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TABLE 3-3 

DEFINITION OF SYMBOLS FOR CHILLR SUBROUTINE DESCRIPTION 


C = R(103) = Thermal mass of cooling coils-dry (Btu/°F) 
c 

C^. = R(95) = Thermal mass of food inner node (Btu/°F) 

C = R(94) = Combined thermal mass of food outer surface node and inner 

° locker walls/shelves (Btu/°F) 

COP = R(IOO) = Coefficient of performance of refrigeration unit (Btu 
cooling provided/Btu input electrical power) 

C = Specific heat (Btu/lb °F) 

P 

C = Specific heat of air (Btu/lb °F) 

^a __ 

D •= R(67) = Operational duty cycle for self-contained refrigeration unit 

(fraction of time on) 

f = R(70) = Fraction of locker void volume which is exchanged with 

ambient air when door is opened 

G = R(lOl) X 6o(^jCp = Thermal conductor (m Cp) of air between 

^ chilling and storage compartments (Btu/hr °F) 

G = R(85) X R(80) = Thermal convection conductor between air in chilling 
^ compartment and cooling coils (Btu/hr °F) 

G = R(78) - unfrozen food \ ^ Thermal conductor between inner and 

^ R(78) X R(76) - frozen food f outer food nodes (Btu/hr °F) 

G = R(89) = Thermal conductor from fan motor to food outer surface and 
locker shel ves/walls (Btu/hr °F) 

G = R(85) X R(81) = Thermal convection conductor between air in food 
® compartment and food/walls outer surface (Btu/hr F) 

G = R(82) = Thermal conductor from cooling coils to food outer surface 
and locker shelves/wal Is (Btu/hr °F) 

6i,Go,Go = R(89), R(90), R(88) = Thermal conductor between fan and ambient 
.123 surroundings by conduction (Btu/hr °F), 

convection (Btu/hr °F) and radiation 
(sq ft), respectively 

j = R(108) = Fraction of coolant pump power transferred to coolant with 
Stirling cycle unit 

m = mass (lbs) 
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TABLE 3-3 (continued) 

DEFINITION OF SYMBOLS FOR CHILLR SUBROUTINE DESCRIPTION 


ill = Mass flow rate (Ib/hr) 


= Heat input to locker air from fan (Btu/hr) 

qp = Electrical power input to the self-contained refrigeration 
unit (Btu/hr) 

qp = Total electrical power input to fan motor (Btu/hr) 

qj, = Total heat rejection rate to the environment by the self-contained' 

" refrigeration unit (Btu/hr) 

q. = Cooling rate provided to the locker by the self-contained 
^ refrigeration unit (Btu/hr) 

qj^ = Heat generated within the fan motor (Btu/hr) 

qp = R(107) X 3.413 = Coolant pump input electrical power used with 
^ Stirling cycle unit (Btu/hr) 

q = Heat loss from locker internal surface and food to ambient 

surroundings through the locker walls, excluding heat transfer 
to the fan (Btu/hr) 

Q, nx = Latent heat of moisture (condensation and/or fusion) in air entering 
locker door (Btu/lb) 

T = R(79) = Air temperature inside chilling compartment (°F) 

T = R(77) = Air temperature inside food compartment (°F) 

®f 

T^ = R(104) = Cooling coils temperature (°F) 

Ti- = R(87) = Effective fan housing temperature (®F) 

T| = R(75) = Temperature of food inner node (°F) 

T = R(51) = Temperature of food outer node and locker internal shelves/ 

° walls (°F) 

T. = A(2) = Inlet fluid temperature to cooling coils (for an externally 

chilled unit) or to chiller heat rejection device (for a self- 
contained refrigeration unit); (°F) 
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TABLE 3-3 (concluded) 

DEFINITION OF SYMBOLS FOR CHILLR SUBROUTINE DESCRIPTION 

T, ,T«,T« = Temperatures of ambient surroundings connected to fan by 
^ conduction, convection, and radiation s respectively (°F 

unless otherwise specified) 

Y = Effective radiation conductor between fan and ambient wall, linearized 

according to equation (3.1.4) 

At = Computing time .step (hours) 

a = Humidity from ambient air condensed inside locker when door is 
opened (lb water vapor/lb dry air) 

Y = Multiplier to compute new temperatures, defined by equation (3.1.10) 
n = Fan motor efficiency (fraction) 

a = Stefan-Boltzmann constant (0-. 1714 x 10“® Btu/hr-sq ft-°R‘*) 

= Density of air (Ib/cu ft) 

a 
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Figure 3-3. Logic Flow Chart for CHILLR Subroutine 
(a) Main Subroutine Logic 
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Figure 3-3- Logic Flow .Chart for CHILLR Subroutine (concluded) 
(b) Fan and Cooling Unit Control Logic 
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TABLE 3-4 

DEFINITION OF CHILLR SUBROUTINE VARIABLES 

C2 = Thermal mass of cooling coils (or other cold sink); (Btu/°F) 

CAIR = Thermal mass flow rate (ifi C„) of air circulated by food locker 

fan (Btu/hr “F) P 

CFOODI = Thermal mass of inner food node (Btu/°F) 

CFOODO = Thermal mass of outer food/inner walls node (Btu/°F) 

CHILR = A logical variable which is TRUE if a self-contained refrigeration 
unit is present and turned on; otherwise, it is FALSE 

CPG = Ambient gas specific heat (Btu/lb °F) 

DELTA = Maximum temperature change of all nodes during a single internal 

iteration (°F) 

DT = Temperature change computed for each node during a single internal 
iteration (°F) 

DTH = System transient time step (hours) 

DTYME = Internal transient time step used (hours) 

FAN = A logical variable which is TRUE if there is a circulation fan 
present and turned on inside the food locker; otherwise, it is 
FALSE 

HC = Heat transfer coefficient; air to food locker walls and food 
surface (Btu/hr-sq ft-°F) 

LFLAG = Logic flag equal to 2 for steady-state case with fan off; 
othervn'se, equal to -2 

LL = Number of transient iterations required for thermal balance 

LPASS = Counter for number of internal iterations of thermal balance equations 

NFLAG = Flag used in logic to compute ambient gas properties 

PWCAB = Water vapor partial pressure in cabin gas (psia) 

PWF = Water vapor partial pressure inside locker (psia) 

Q51 = Total heat transferred into outer food node (Btu/hr) 

Q78 = Hfest transferred from food inner to outer node (Btu/hr) 

QCOOL = R(65) = Net cooling provided to locker (Btu/hr) 
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TABLE 3-4 (concluded) 

DEFINITION OF CHILLR SUBROUTINE VARIABLES 

QDOOR = Heat added during current iteration from opening door (Btu/hr) 

QFAN - Heat from food locker fan which is dissipated within the fan 
housing (6tu/hr) 

QFANA = Heat from food locker fan work of compression of the gas within 
the locker (Btu/hr) 

QFANF = Heat transfer from fan housing to locker inner walls and food 
outer surface (Btu/hr) 

QFREEZ = Latent heat required to freeze all the liquid water in food 
inner or outer nodes (Btu) 

QLAT = Latent heat of moisture (condensation and/or fusion) in air 
entering locker door (Btu/lb) 

QREJ = Heat rejected by self-contained cooling unit to cabin gas (Btu/hr) 

R51 = Variable used for temporary storage of R(51) during thermal balance 

R98F = Absolute humidity of gas inside locker (lb water vapor/lb dry air) 

RHOG = Ambient gas density (Ib/cu ft) 

SG2 = Effective sum of thermal conductors to cooling coils (or other 
cold sink); (Btu/hr-°F) 

SG51 = Effective sum of thermal conductors to outer food node (Btu/hr-°F) 

SG75 = Effective sum of thermal conductors to inner food node (Btu/hr-°F) 

TC = Temperature of cooling coils (or other cold sink); (°F) 

TG = Ambient gas temperature (°F) 

X ' / = Multiplier in steady-state heat balance equation (3.1,13) 

cc = Difference in absolute humidity between ambient air and air 
inside chiller locker (lb water vapor/lb dry air) 


3-37 


D2-11857\-' 


3.2 FTRAY 

The FTRAY subroutine is designated as 6-189A No. 66. 


3.2.1 Subroutine Description 

This subroutine simulates the performance of a Skylab-type food warming/ 
serving tray. A typical food tray, shown in Figure 3-4, had eight recessed 
food cavities , .of which three had embedded thermostatically controlled 


electrical resistance heaters to 
warm the food. The thermal model 
used for a single food warming 
cavity is shown in Figure 3-5. 
The cavity is assumed cylindrical 
and is subdivided into five food 
nodes of equal volume. In addi- 
tion to the thermal network shown 
in Figure 3-5, the food is 
thermally connected to the ambi- 
ent surroundings using the G-189A 
subroutine QSURR. 



Tray 


3.2.2 Subroutine Data 

3. 2.2.1 General Notes 

a. Cabin air flow over the food 'tray component is optional. If an upstream 
component is used to provide air flow, it must be connected to the food 
tray primary side. All flow codes are allowed. 

b. The total heat lost from all food trays to the ambient surroundings is 
given in R(53). This. is- the sum of convection [R(56)], radiation [R(59)], 

■ • and conducti■pn^[•R(.62')]. If air flow is provided by an upstream component, 
the convective'heat ■[R(56)] is transferred to that flow. All heat loss 
not rejected to the primary side air flow should be added to the cabin 
sensible heat load [cabin component R-array location R(66)] in GPOLYl. 
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3. 2. 2.1 (Continued) 

c. Any number of heated food cavities may be simulated by a single G-189A 
component using FTRAY. The routine determines the performance of a 
single food cavity and assumes all others are identical. If some 
food trays have different input data or time schedules, they must be 
simulated by separate G-189A components. 

3. 2.2.2 Instruction Options 

NSTR(l): Initial conditions for the five food nodes 

= 0 Assume properties for all five nodes initially equal, and input 
temperature and v^ater content for node 1 only in R(71) and R(96) 

= 1 Input initial temperature and water content for all five nodes 
in R(71) through R(75) and R(96) through R(IOO) 

3. 2. 2.3 Heat Loss V-Array Data 

Reference 

Location Description Data Type 


51 

Average temperature of the five food 
nodes (°F) 

O' 

52 

Summed thermal conductance from one 
food cavity to ambient surroundings 
(Btu/hr °F) 

0 

53 

Total heat loss from all heated food 
cavities to surroundings (Btu/hr) 

0 

54 

Ambient gas temperature (°F) 

1(0) 

55 

Thermal convection conductor from 
one food cavity to ambient gas 
(Btu/hr °F) 

1(0) 

56 

Convective heat loss from all heated 
food cavities to ambient gas (Btu/hr) 

0 

57 

Ambient wall temperature (°F) 

1(0) 
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3. 2. 2. 3 (Continued) 


Reference 

Location 

Description 


Data Type 

58 

Thermal radiation conductor (A3) between 
one food cavity and ambient walls 
(sq ft) 

1(0) 


59 

Radiative heat loss from all heated 
food cavities to ambient walls (Btu/hr) 

0 


60 

Temperature of ambient structure 
connected to food node (°F) 

1(0) 


61 

Thermal conductor from one food node to 
attached ambient structure (Btu/hr °F) 

1(0) 


62 

Thermal heat loss from all heated 
food cavities to attached ambient 
structure (Btu/hr) 

0 


3. 2. 2.4 

Steady-State K-Array Data 



Riiference 

Location 

Description 


Data Type 

16 

Not used 



17 

Not used 



18 

Ice/water transition indicator. 

0 = water in all food nodes is 
frozen. 1 = some food nodes are in 
ice/water transition. 2 = water in 
all food nodes is liquid. 

0 


3. 2. 2. 5 

Steady-State V-Array Data 



Reference ... 

Location Description 


Data Type 

65 

Total heater power for all heated 
food cavities (Btu/hr) 

0 


66 

Number of food cavities to be heated 

I(R) 


67 

Single food cavity radius (in). 

1(0) 
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3. 2. 2. 5 

(Continued) 



Reference 
Ldcati on 

Description 


Data Type 

68 

Single food cavity depth (in). 


1(0) 

69 

Single food cavity average heater 
power (watts) 


1(0) 

70 

Maximum sum of thermal conductors 
any one food node (Btu/hr °F) 

to 

0 

71 

Temperature at food node 1 (°F) 


1(0). 0 

72 

73 

74 

75 ■ 

Temperature at food node 2 (°F) 
Temperature at food node 3 (°F) 
Temperature at food node 4 (°F) 
Temperature at food node 5 (°F) 


I(R) if NSTR(1)=1 
otherwise, these i 
automatically set 
equal to R(71), 0 

76 

Thermal conductor from food node 
to 4 (Btu/hr °F) 

3 

0 

77 

Thermal conductor from food node 
to 3 (Btu/hr °F) 

2 

0 

78 

Thermal conductor from food node 
to 2 (Btu/hr °F) 

1 

0 

79 

Thermal conductor from food node 
to 5 (Btu'/hr °F) 

1 

0 

80 

Fraction of heater power on food 
node 4 (or 5) 


0 

81 

Fraction of heater power on food 
node 1 


0 

82 

Frozen food thermal conductivity 
(Btu/hr-ft-°F) 


1(0) 

83 

Unfrozen food thermal conductivity 
(Btu/hr-ft-°F) 

1(0) 
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3. 2. 2.6 Steady State K-Array Data Used for Transient Case 
Reference 

Location Description Data Type 

16 Food node number to be used for 1(0) 

automatic temperature control, 

1 through 5 

17 This parameter stores the heater ON/ 0 

OFF mode on previous time step, 

and is used for automatic temperature 
control. 0=0FF, 1=0N 

18 Ice/water transition indicator. 0 

0 = water in all food nodes is 

frozen. 1 = some food nodes are in 
ice/water transition. 2 = water in 
all food nodes is liquid. 


3. 2. 2.7 Transient V-Array Data 


Reference 

Location Description Data Type 


65 

66-68 

69 

70-79 

80 

81 

82-83 

84 

85 


Total average heater power for all 
heated food cavities during previous 
time step (Btu/hr) 

Same as for steady state 

Not used 

Same as for steady state 

Fraction of heater power and thermal 
mass assigned to food node 4 (or 5) 

Fraction of heater power and thermal 
mass assigned to food node 1 

Same as for steady state 

Heater control temperature (°F) 


Heater control temperature dead band 

(°F) fT ,=R(84) ± 

R(85)l control 


0 



1 ( 0 ) 
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3. 2. 2. 7 

(Continued) 

, 

Reference _ . _ 

Location 

Description 

Data lype 

86 

Specific heat of unfrozen food (Btu/lb 
wet food-°F) 

1(0) 

87 

Specific heat of bone-dry food (Btu/lb 
dry food-°F) 

0 

88 

Effective thermal mass of food heating 
tray connected to a single heating 
cavity (Btu/°F) 

1(0) 

89 

Moisture content of food (lb H20/lb 
total wet food weight) 

1(0) 

90 

Maximum heater power for a single 
heating cavity (watts) 

1(0) 

91 

Thermal mass of food node 1 (Btu/°F) 

0 

92 

Thermal mass of food node 2 (Btu/°F) 

0 

93 

Thermal mass of food node 3 (Btu/°F) 

0 

94 

Thermal mass of food node 4 (Btu/°F) 

0 

95 

Thermal mass of food node 5 (Btu/°F) 

0 

96 

Fraction of moisture in food node 1 \ 
which is frozen 


97 

Fraction of moisture in food node 2 
which is frozen 

^ I(R) only if initial 
food temp, is 32.0°F,0 

98 

Fraction of moisture in food node 3 

which is frozen 

99 

Fraction of moisture in food node 4 


which is frozen 


100 

Fraction of moisture in food node 5 
which is frozen 


101 

Food mass density (Ib/cu ft) 

1(0) 

102 

Dry food mass in a single food node 
(lbs) 

0 

103 

Mass of moisture in a single food 
node (lbs) 

0 
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Z.2.2,1 (Continued) 

Reference 

Location Description Data Type 

104 Food tray ON-OFF svn'tch. Set to -1.0 1(0) 

if tray is turned off, and to +1,0 if 

^.tray is in -a warming cycle. Note: Power 
to tray may be internally turned off, 
even when this switch is on, due to 
‘heater thermostat controller. Each 
time tray is turned on, reinitialize 
food temperature according to NSTR(l). 

105 Stores value of R(104) on previous 0 

iteration; used for ON/OFF initiali- 
zation logic. 


3.2.3 Analytical Model Description 
3. 2. 3.1 Geometric Considerations 

A single cylindrical food heating cavity is divided into five food nodes, 
as shown in Figure 3-5. Based on the assumption that all five nodes have 
the same volume, the following parameters based on geometry may be com- 
puted (see Table 3-7 for definition of symbols): 


Volume 


Height 


Vi = V2 = V3 = = V5 = 0.2 TF R^H 


hj = O'. 2 H 

h2 = hs = h^ = hg = 0.4 H 


(3.2.1) 


(3.2.2) 

(3.2.3) 
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3.2.3. 1 (Continued) 

Heater Areas 

Am = 7t R2 + 2 it RH 


Am = A =0 
^2 


Am “Am = 0.8 17 RH 

H4 Hg 


Am = 7T R2 + 0.4 ^ RH 
“1 


Outer Node Radii 


ri = ^4 = r-j = R 


rg = Tg = 0.707 R 


Radii of Node Centers (placed" at geometric center of node volume) 


r^ = 0.707 R 


i( ic 


r^* = rg* = 0.866 R 


Length Between Node Centers 


L 34 = l 2 g = 0.366 R 


(3.2.4) 

(3.2.5) 

(3.2.6) 

(3.2.7) 

(3.2.8) 

(3.2.9) 

(3.2.10) 

(3.2.11) 

(3.2.12) 

(3.2.13) 
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3. 2.3,1 (Continued) 


^23 " ^45 ■ ^ 


Lj 2 = ^ 0.09 H2 + 0.04285 R2 


Ljg = V 0.09 H2 + 0.02528 R2 


"Window" Area Between Nodes 


A 34 = A 25 = 0.5656 TT RH 


^23 “ ^45 “ ^12 ~ ^15 0.5 IT 


Heater Flux (?fe total heater power) 


Total Heater Area = 2 tt RH + ir R^ 


' ''Hj 


1 


2.5 + 1.25 R 
H 






0.4 + R 

H 

2 + R 
H 




Thermal Conductors (k = food thermal conductivity) 


G 34 - 622 = 1.548 TT Hk 


(3.2.14) 

(3.2.15) 

(3.2.16) 

(3.2.17) 

(3.2.18) 


(3.2.19) 


(3.2.20) 


(3.2.21) 



1.25 TT 



(3.2.22) 
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3.2.3. 1 (Continued) 

(3.2.23) 

(3.2.24) 

3. 2. 3. 2 Thermal Balance 

The thermal solution is divided into two parts: that within the food 

cavity, and thermal exchange between the food and ambient surroundings.. 

These two solutions are interdependent since: 

a. .The solution within the food cavity provides the average food temperature 

for computing ambient heat loss. 

b. The thermal balance with ambient surroundings provides the boundary 
conditions fo^ heat transfer within the food cavity. 

The heat loss from the food to ambient surroundings is found first using 
the G-189A subroutine QSURR. This routine determines the total conductive, 
radiative, and convective thermal loss from the food. The arithmetic 
average temperature of the five food nodes is used as the effective base 
temperature in the subroutine. The top food surface is assumed to be 
covered with a lid, as on Skylab, with negligible moisture evaporation. 

Thus, the boundary conditions for the food cavity are: 

a. The heat lost at the heater surface equals the total conductive heat 

(q ,) from the food to surroundings; 

^^cond' 

b. The heat lost at the top surface equals the total radiative (q^) and 

convective (l^onv^ surroundings. 


0.5 TT R^k 

yfo.09 H2 + 0.0429 R^ 

0.5 TT R^k 

Vo .09 H2 + 0.0253 R2 
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3. 2. 3. 2 (Continued) 

The thermal balance within the food cavity may then be made by applying 
the following standard equations at each food node: 





for a transient solution, and 


(3.2.25) 



(3.2.26) 


for a steady-state solution. The net heat into each food node is 
given in terms of the nodal surface areas (equations 3.2.17 and 3.2.18), 
heater areas (equations 3.2.5 and 3.2.6), heater flux (equations 3.2.18 
and 3.2.19) j and thermal conductors (equations 3.2.20 through 3.2.23), 

as follows: 


'^1 " \ 



(Tj - Tj) + (Tj 



(3.2.27) 


Z ^2 " ‘^12 ^"^1 ■ "^ 2 ^ ■^■'^23 "^ 2 ^ ®25 ^"^5 " "^ 2 ^ 


(3.2.28) 


' \onv^ ^ ^3 ("^2 ■ "^3^ *^34 ^''‘4 " ^3^ 



('^r’^^conv^ *^34 ^'^3"’’'4^ ^45 ^^5' 


"^cond ^15 '^25 ^’’’2"^5^ *^45 


(3.2.29) 

T^) (3.2.30) 

(3.2.31) 
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3. 2. 3.2 (Continued) 

Equations (3.2.25) through (3.2.31) are solved each time step using a 
standard forward-difference method as follows: 

T, = T, + 7 Z (3-2-32) 

^new ^old 

For steady-state solutions, the value used for y is given by 

0.8 ( 3 . 2 . 33 ) 

'^steady state maximum sum of conductors to ground 

Internal steady-state iterations are made using equations (3.2.32) and 
(3.2.33) until the old and new temperatures differ by 0.05°F or until 

25 iterations have been made. 


For transient solutions, the value for y in equation (3.2.32) is 


^transient 



(3.2.34) 


The thermal mass m c of the node will vary with the ice/water balance, 
as discussed further in the following section. The computing time incre- 
ment At is first computed as 0.'7 times the smallest of the time increments 
required for computational stability, where 

^ thermal mass (3.2.35) 

^^stability “ sum of conductors to ground 

The time increment is further limited to the system computing time step. 

In addition, it was found that large time steps could cause the automatic 
heater control to leave the heater on too long; thus , ‘-causing excessive 
temperature fluctuations in the food nodes adjacent to the heater. There- 
fore, as suggested in Reference 4, the internal time increment was further 
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3. 2.3. 2 (Continued) 

limited to 20 seconds. If necessary, internal iterations are performed 
until the required system time step is reached. 


3. 2. 3. 3 Thawing Effects 

Some thermal properties for typical foods, taken from Reference 4, are 
given in Table 3-6. It is apparent that some properties are highly dependent 
on the physical state (whether liquid or ice) of the food moisture. Thus, 
the thermal mass for each node is computed separately: 


'p, ' <Vy food "liquid water "ice 


(3.2.36) 


As a frozen food node reaches 32°F, its temperature is held constant and 
the heat input assumed to thaw its moisture. The water/ice balance during 
this transition is computed for evaluating thermal mass from equation 
(3.2.36). The food thermal conductivity is assumed to be that for frozen 
food at all nodes until all the water is melted. Then the conductivity 

for thawed food is used. 


3. 2.3.4 Heater Control 

For a transient solution, the heater power is controlled automatically to 
an input set point temperature. Any node may be used as a temperature 
sensor, but normally node 1 should be used since it will be hottest. A 
control deadband is also used. Since the internal computing time incre- 
ment may be less than the system time step, as discussed in Paragraph 3. 2. 3. 2, 
the heater may be cycled on and off several times during one transient time 
step. The integrated average heater power during one time step is computed 
and output in R-array location R(65). For a steady-state solution, only 
the average heater power should be input, and the heater assumed to be con- 
tinually energized. 
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3. 2. 3. 5 Pressure Drop Considerations 

The food warming/serving trays are normally exposed only to the cabin air 
convection currents and are not physically connected in an actual fluid 
flow loop. Thus, no pressure drop calculations would be required for this 
component. 
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TABLE 3-5 

DEFAULT VALUES FOR FTRAY SUBROUTINE INPUT DATA 


K(16)=l 


R-Array* 

Location 

54 

55 

57 

58 
60 
61 

.67 

68 

69 

71 


82 

83 

84 

85 

86 
88 

89 

90 
101 
104 


Default 

Value 


70 

0.148 

70 

0.0625 

70 

0.216 

1.875 

1.125 

7.9 
140 

- 10 
70 

. . 0.99 
0.25 
157 
2 

0.78 

1.9 
0.724 

52 

59 

1 


steady state 
transient, R(104) > 0 
transient, R(104) < 0 


*See Section 3.2.2 for definition of variables. 
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TABLE 3-6 

THERMAL PROPERTIES OF SELECTED SKYLAB FOODC 


Food 

Substance 

Water 

Content 

(2) 

Density 

lb 

7U 

Latent 
Heat of 
Fusion 
Btu 
lb 

a 

Specific 

Heat 

(Btu/lb -F) 

** 

Themal 

Conductivity • 
(Btu/hr-f t-F) 

** + 

. 

j 

Theraal | 

Diffusivitv [ 

(ft^/hr)' 1 

1 

|_ 

Pre-but- ■ 
tered Roll 

— 

23.8* 

15.* 

34. 

0.70 0.34 

0.19 0.45 

0.018 0.088 1 

Coffee 

Cake 

25.0* 

15.* 

36. 

0.70 0.34 

0.19 0.46 

' 

0.018 0.088 1 
» 

0.0063 0.042 i 
1 
t 

Filet 

Mlgnon 

63.2* 

58. 

91. 

0.71 0.40 

$ 

0.26 0.97 

* 

Chili with 
Heat 

66.9* 

58. 

96. 

0.74 

0.40 

0.26 0,99 

0.0061 0.043 j 
» 

Prlrac Rib 
of Beef , 

72.4* 

59, 

104. 

0.78 0.42 

0.27 1.04 

’ 

0.0059 0.042 1 
1 

Lobster 

Kewberg 

75.5* 

59. 

109. 

0.80 0.43 

0.28 1.09 

0.0059 0.043 1 
0.0054 0.044 1 

Stewed 
1 Tooatoas 

83.4 • 

61. 

129. 

0.91 

0.47 

0.30 1.25 


*Espericcutally detcnained for actual Skylab foods 
**Dn£ro:en 
tFrozea 


e POOR 



m OF TM 
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TABLE 3-7 

DEFINITION OF SYMBOLS FOR FTRAY SUBROUTINE DESCRIPTION 
' Am = Area of heater contact with food node "i" (sq ft) 

Am = Total heater surface area (sq ft) 

”t 

A. . = "Window" cross-sectional area between adjacent nodes "i" and "j" 

’J (sq ft) 

c = Thermal specific heat of node "i" (Btu/lb °F) 

Pi 

G. . = Thermal conductor between adjacent food nodes "i" and "j" (Btu/hr °F) 

• w 

H = Food cavity height (ft) 

= Total heater input electrical power (Btu/hr) 
h^. = Height of food node "i" (ft) 

k = Food thermal conductivity (Btu/hr-ft-°F) 

L. . = Length between centers of adjacent nodes "i" and "j" (ft) 

• J 

m = Total mass of food node (lb) 

*^cond " conduction from food to ambient surroundings (Btu/hr) 

qconv " convection from food to ambient surroundings (Btu/hr) 

qjj = Heater energy input into node "i" (Btu/hr) 

”i 

= Heat flow into node "i" (Btu/hr) 

= Heat lost by radiation from food to ambient surroundings (Btu/hr) 

R = Food cavity radius (ft) 

r^. = Outer radius of food node "i" (ft) 

r^ = Radius of center of food node "i" (ft) 
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TABLE 3-7 (concluded) 

DEFINITION OF SYMBOLS FOR FTRAY SUBROUTINE DESCRIPTION 

T^. = Temperature at node "i" (°F) 

At = Computing time increment (hrs) 

V^. = Volume of food node "i“ (cu ft) 

Y = Multiplier to compute new temperatures, defined by equation (3.2.32) 
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Figure 3-6. Logic Flow Chart for FTRAY Subroutine (concluded) 


OF 

PA&B is PWM 
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TABLE 3-8 

DEFINITION OF FTRAY SUBROUTINE VARIABLES 
ASP = Food cavity aspect ratio, radius/height (dimensionless) 

s 

COND = Food effective thermal conductivity (Btu/hr-ft-°F) 

DELTA = Maximum temperature change of all nodes during a single internal 
iteration (°F) 

DT = Temperature change computed for each food node during a single 
internal iteration (°F) 

DTH = System transient time step (hours) 

DTYME = Internal transient time step used (hours) 

GMAX = Maximum sum of thermal conductors to any food node (Btu/hr °F) 

H = Food cavity height (ft) 

JF = Number of food nodes with all water frozen 

JL = Number of food nodes with all water liquid 

LL = Maximum number of internal iterations of heat balance equation 

LPASS = Counter for number of iterations of thermal balance equations 

Q(l) >Q(2) ,Q(3) ,Q(4) ,Q(5) = Net heat transferred into food nodes 1, 2, 3, 

4, and 5 (Btu/hr) • 

QA3 = Heat transferred from ambient air into food node 3 (Btu/hr) 

QHTR = Net heat into food cavity from electrical heater (less the heat 

conducted into the tray); (Btu/hr) 

Q^.j = Heat transfer from food node "1" to "j" (Btu/hr) 

QMELT = Latent heat required to melt the ice in each food node (Btu) 

QW4 = Heat transferred from heater wall into food node 4 (Btu/hr) 

RAD = Food cavity radius (ft) 

SG1,SG2,SG5 = Sum of thermal conductors to food nodes 1, 2, and 5 (Btu/hr-°F) 
T1,T2J3,T4,T5 = Temperature of food nodes 1, 2 , 3 , 4, and 5 (°F) 

TAVG = Arithmetic average of the five food node temperatures (°F) 
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TABLE 3-8 (concluded) 

DEFINITION OF FTRAY SUBROUTINE VARIABLES. 

WDRY = Weight of dry food only in a single food node (lb) 

WH20 = Weight of water in a single food node (lb) 

X = Y = Multiplier in heat balance equation (3.2.32) 
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3.3 ROSHOS 

The RQSMOS subroutine is designated as G-189A No. 69. 

3.3.1 Subroutine Description 

Reverse osmosis is a process used for removing impurities from waste water. 
Fluid is forced by static pressures across a membrane in a direction opposite 
to which it would normally flow due to osmotic pressure alone. Most of the 
impurities are filtered out by the membrane, thus leaving a purified water 
flow through the unit. This technique is in use commercially in many water 
and waste treatment applications and is being developed to recycle spacecraft 
waste water. The subroutine will simulate the performance of an arbitrary 
type of reverse osmosis unit, shown in Figure 3-7, if test or design per- 
formance data for that unit are available. Only a single-stage unit may be 
simulated by a single component, but multistage operation may be achieved 
by connecting more than one unit together. 



3.3.2 Subroutine Data 
3. 3. 2.1 General Notes 

a. Feedwater must be inlet on the primary side. Treated product water 

will be outlet on the secondary side, and the brine outlet at the primary 
side. 


3-61 


D2-II8571-I 


3. 3. 2.1 (Continued) 

b. Only liquid flow codes (0 or 4) may be used. Results for flow code 4 
will only be correct if the recovery factors for each special flow 
type are the same. 

c. The reverse osnir is unit is unique among G-189A components in that 

the secondary (product) outlet pressure must be calculated every time to 
obtain a solution, even when normal pressure drop calculations are not 
being made. To do this, the user must input the value to be used for 
product water pressure in GPOLY logic in R(74). 


3. 3. 2. 2 Heat Loss V-Array Data 


Reference 

Location Description Data Type 

51 Average reverse osmosis unit core 1(0), 0 

temperature (°F) 


52 Summed thermal conductance between 0 

reverse osmosis unit and ambient - 
excluding conductance to liquid flow 
through unit (Btu/hr °F) 

53 Reverse osmosis unit net heat loss to 0 

ambient - excluding heat transfer to 
liquid flow through unit (Btu/hr) 

54 Ambient gas temperature (°F) 1(0) 


55 Thermal convection conductor from I(R) 

reverse osmosis unit outer surface to 
ambient gas (Btu/hr °F) 


56 Convective heat loss from reverse 0 

osmosis unit to ambient gas (Btu/hr) 

57 Ambient wall temperature (°F) 1(0) 

58 Thermal radiation conductor (A3-) I(R) 

between ambient walls and reverse 

osmosis unit surface (sq ft) 

59 Radiative heat loss from reverse 0 

osmosis unit to ambient walls (Btu/hr) 
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3. 3.2.2 (Continued) 



Reference 

Location 

Description 

Data Type 

60 

Temperature of ambient structure 
connected to reverse osmosis unit (°F) 

1(0) 


61 

Thermal conductor from reverse 
osmosis unit to attached ambient 
structure (Btu/hr °F) 

I(R) 


62 

Thermal heat loss from reverse 
osmosis unit to attached ambient 
structure (Btu/hr) 

0 


63 

Reverse osmosis unit external 
surface temperature (°F) 

0 


64 

• 

Thermal conductor throughinsulation 
around reverse osmosis unit (Btu/hr Fj; 
may be zero if unit is uninsulated 

l(R) 


3. 3. 2. 3 : 

Steady-State K-Array Data 



Reference 

Location 

Description 


Data Type 

16 

Number of individual types of s^.lts_ 
to be handled by reverse osmosis unit 
(extra V-array locations must be five 
times this value) 

1(0) 


3. 3. 2. 4 

Steady-State V-Array Data 



Reference ^ 

Location Descriptip_n_ 


Data Type 

65 

Overall design water recovery factor, 

m j u./mx j (dimensionless) 

‘"producr feed ^ 

1(0) 


66 

Overall design total salt rejection 
factor 

1(0) 


67 

Design product water pressure (psia) 

1(0) 


68 

Brine overall total salt concentration 
(lb salt/lb water) 

0 
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3. 3. 2. 4 (Continued) 

Reference 

Location Description Data Type 

69 Product water overall total salt 0 

concentration (lb salt/lb water) 

70 Design feedwater flow rate (Ib/hr) 1(0) 

71 Design feedwater overall total salt I(0) 

concentration (lb salt/lb water) 

72 Actual feedwater overall total salt 1(0) 

concentration (lb salt/lb water) 

73 . ‘Design^feedwater pressure (psia) 1(0) 

74 '- Actual product pressure (psia) input 1(0) 

in.-GPOLYl from downstream component 

75 Overall design average osmotic 1(0) 

pressure (psia) 


3. 3. 2. 5 Transient V-Array Data 


Reference 

Location Description Data Type 

76 Thermal capacitance -of reverse 1(0) 

osmosis unit-dry '(Btu/°F) 


3. 3. 2. 6 Extra V-Array Data [need 5J locations where J=K(16)] 


Reference 

Location Description Data Type 


77-*76 + 0 

Design concentration of individual 

1(0) 


feedwater salts (lb salt/lb water) 


77 J to 

Actual concentration of individual 

1(0) 

76 + 2J 

feedwater salts (lb salt/lb water) 


77 + 2J to 

Design individual salt rejection 

1(0) 

76 + 3J 

factor (dimensionless) 
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3.3. 2.6 (Continued) 

Description Data Type 

Actual concentration of individual 0 

salts in product water (lb salt/lb 
water) 

Actual concentration of individual 0 

salts in brine (lb salt/lb water) 

3.3.3 Analytical Model Description 

In a reverse osmosis unit, fluid is forced by static pressures across a 
membrane in a direction opposite to which it would normally flow due to 
osmotic pressure alone. This technique is in use commercially in many 
water and waste treatment applications and is being developed to recycle ^ 
spacecraft waste water. Since a particular spacecraft reverse osmosis unit 
design has not been determined, the subroutine was written in general terms 
to handle any unit for which some design or test data are available. The 
input design data are adjusted for the off-design conditions present during 
a run. The thermal balance of the reverse osmosis unit with ambient sur- 
roundings is handled using the standard G-189A subroutine QSURR. 


3.3. 3.1 Flow Balance 

The flow balance across a reverse osmosis module may be analyzed from two 
viewpoints of interest. First, one can consider the overall net effect of 
all the different types of solute impurities, lumped together and treated 
as a single homogeneous impurity. This approach is used in determining 
overall v^ater balances and initial sizing and design of water recovery com- 
ponents. Secondly, one can consider the effect of the reverse osmosis 
module on each different type of impurity present in the water with its own 
individual rejection factor. This subroutine uses the former "overall" 
approach for handling the total flow balance, and will also apply the second 
approach as an option for handling any number of individual Impurities 

desired. 


Reference 

Location 

77 + 3J to 
76 + 40 


77 + 40 to 
76 + 50 
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3.3.3. 1 (Continued) 

The overall water recovery factor for a reverse osmosis unit is defined 
s as 


Y 



(3.3.1) 


For conservation of mass across the unit, 

mp mg ~ mp 


(3.3.2) 

(3.3.3) 


where m refers to the total mass flow rate, and S is the effective 
overall mass flow rate of the solutes. The overall salt rejection factor 
for the unit is defined as 

§ 

R = (3.3.4) 

^F 


Combining equations (3.3.1) through (3.3.4), the overall total salt 
concentration in the brine and product water is given as a function of the 
inlet conditions as follows: 



R Cp 
(1 - Y) 


(3.3.5) 


Sp (1 - R) Cp 


(3.3.6) 
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3.3.3. 1 (Continued) 


Now the osmetic pressure of a given substance through a membrane has been 
^ found to vary with the concentration of that substance. It is assumed a 


proportional relationship exists such that 


^^osmosis '^^^osmosis 


design \^F 


design 


(3.3.7) 


The product water flow rate is assumed proportional to pressure drop across the 
membrane less the osmotic pressure; thus. 


mn P r 


aP. 


osmosis 


(3.3.8) 


The design pressures and flow rates of the reverse osmosis unit are input. 
Therefore, the flow rate of the product water may be found from equation (3.3.8) 
as follows: 



design 


(Pr - Pp - aP . \ 

Vf P osmosis; 

/^F ~ " ^^osmosisl . 

\ y design 


(3.3.?) 


To use equation (3.3.9) to compute the product water total flow rate ifip, the 
product pressure Pp must first be found. For an actual unit, this would 
be a function of the feedwater pressure, or it could be regulated by a down- 
stream valve or pump. Since these data are not known, the user must supply 
the product water pressure in a GPOLY statement based on data for the actual 
unit being simulated. The design product flow '^P. . in equation (3.3.9) 

is found from equation (3.3.1): 


mn - Y , . nV 
design esign f^^esign 


(3.3.10) 
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3.3.3. 1 (Continued) 

It is assumed that the flow rate of salts in the product will be proportional 
to (1) the feedwater total salt concentration , -and (2) the effective pressure 
drop across the membrane; that is, 

(3.3.11) 

Combining equations (3.3.5), (3.3.9), and (3.3.11), the overall concentration 
of salts in the product water may be found: 



(3.3.12) 


The design product concentration is found from equation (3.3.6) using input 
design values for R, Cp, and Y. Finally, the overall concentration of 
salts in the brine is found from equations (3.3.2) and (3.3.5): 




(3.3.13) 


3. 3. 3. 2 Individual Salt Concentration Balance 

The overall flow equations in Paragraph 3. 3.3.1 are sufficient 
to compute the outlet primary and secondary flows for the G-189A system. 
However, as an option, the subroutine will compute the product water and 
brine concentrations of any number of individual specific solutes. The 
same analysis in Paragraph 3. 3. 3.1 for the overall water/salt balance may 
be made for each individual salt with similar results. This would give, 
fO‘" the brine and product individual salt concentrations [see equations 

(3.3.12) and (3.3.13)] : 




design 



) ; i=l, total number 
of salts 


(3.3.14) 
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3. 3. 3. 2 (Continued) 


m - m 

- E E E — t i=i, total number 

io m of salts 


( 3 . 3 . 15 ) 


Equations (3.3.14) and (3.3.15) will be evaluated in the subroutine for 
each individual solute requested by the user. For each solute, the design- 
values (or test values if available) must be input for feedwater salt con- 
centration and salt rejection factor . Also, the actual feedwater 
concentration V each solute must be computed in GPOLY logic and input. 
Again, the design product concentration is found from equation (3.3.6): 


(1-Ri, . )S- 

desman 


(3.3.16) 


design 


'design 


3. 3. 3. 3 Thermal Balance 

Thermal exchange between the reverse osmosis unit and ambient surroundings 
is accounted for using the standard G-189A subroutine QSURR. This model 
was included to accommodate certain designs requiring elevated temperatures, 
and simulates conduction, radiation, and convection from the core and insu- 
lation effects. The net heat into the single core node is given by 


^in " "’f 


(Tp - TJ - q, 


For the steady-state case, the net heat input must be zero; so the average 
core temperature is given by 


steady state 


nip Cp 
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3. 3. 3. 3 (Continued) 

For transient solution, the core temperature is found as follows: 


°new %ld Cp 


At 


(3.3.19) 


The computing time increment At is computed as 0.5 times the maximum 
time allowed for computational stability, where 


At 


core thermal mass 


stability “ sum of conductors to ground 


(3.3.20) 


If the system time step is less than the stability time increment, then 
the system time step is used directly for At. If it is greater, then 
internal iterations of equation (3.3.19) are repeated until the system 
time step is reached. 


3. 3. 3. 4 Pressure Drop Considerations 

The reverse osmosis unit has one flow inlet (primary feedwater) and two 
outlets (primary brine and secondary product water). This component is 
still in the design stage, and pressure drop model data cannot yet be 
specified. Typical designs now being considered have an approximate 
pressure range of 1000 psia feedwater, 50 psig drop to the brine, and 
50 psia product water pressure. The feedwater flow rate may or may not 
be controlled by a downstream pressure regulator. If it were, the brine 
pressure should be set in GPOLY logic based on the regulator setting. This 
would not require use of the standard G-189A pressure drop options. If 
a regulator were not used, the G-189A option (1) would again probably be 
best to use for computing brine pressure (with the effective flow path 
area, length, and hydraulic diameter being input) unless actual experi- 
mental data were available. 
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3. 3. 3. 4 (Continued) 

The reverse osmosis unit is unique among G-189A components in that the 
^ product water (secondary outlet) pressure must be known in order to obtain 
a solution. Thus, it must be either input or calculated every time step 
whether a system pressure drop analysis is being made or not. Since the 
product is on the component secondary side, and there is no secondary 
inlet flow (feedwater is on the primary side), the standard G-189A pressure 
drop options may not be used to compute its pressure. Thus, the product 
pressure must be input or computed in GPOLY logic based on experimental 
data. In the Space Station system model described in Section 6, the product 
pressure is specified in R(74) as 45 psia. When the reverse osmosis unit 
design becomes finalized and pressure drop data are available, it is 
recommended they be incorporated into the ROSMOS subroutine; thus, elimi- 
nating this additional user input. 
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TABLE 3-9 

DEFAULT VALUES FOR ROSMOS SUBROUTINE INPUT DATA 


R-Array* 

Location 

Default 

Value 

51 

75 

54 

70 

57 

70 

60 

70 

65 

0.96 

66 

0.93 

67 

50 

70 

9 

71 

. 0.01 

72 

0.008 

73 

900 

74 

45 

75 

200 

76 

2.5 


*See Section 3.3.2 for definition of variables. 
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TABLE 3-10 

DEFINITION OF SYMBOLS FOR ROSMCS SUBROUTINE DESCRIPTION 

i 

C = Solute concentration §/iri (lb solute/lb total fluid) 

Cp = Thermal specific heat (Btu/lb °F) 

m = Fluid mass flow rate (Ib/hr) 

P = Fluid pressure, psia 

*^surr ~ reverse osmosis unit to ambient surroundings 

(excluding heat transfer to the flowing fluid^- (Btu/hr) 

R = Overall total solute rejection factor defined by equation (3.3.4) 

S = Effective overall solute mass flow rate (Ib/hr) 

T = Temperature (°F) 

Y = Overall water recovery factor defined in equation (3.3.1) 

At = Computing time increment (hours) 

Subscripts 

B = Brine 

F = Feedwater 

i = Individual solute 

0 = Reverse osmosis unit core 

P = Product water 
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■W INITIALIZE COMPONENT DATA 


COMPUTE TOTAL OUTLET 
FLOW RATES R(l), R(20) 


DETERMINE OUTLET 
product PRESSURE 


f 

DO 50 L 

= 15, 19 

SET OUTLET CONSTITUENT 
FLOW RATES 


COMPUTE OUTLET OVERALL BRINE 
AND PRODUCT SALT CONCENTRATION 



DO 70 L=l, K(l) 


COMPUTE' INDIVIDUAL OUTLET 
SALT CONCENTRATIONS 



RETURN 


3 



DETEPJ-IINE 
COMPUTING 
TIME INCREMENT 
DTYME 




i 


f 

0.5 

.A . . .... 


nTVMP 

^ rG 

••4 



c 


LPASS=0 
T 



CF Figure 3-8. Logic Flow Chart for ROSMOS Subroutine 
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Figure 3-8. Logic Flow Chart for ROSMOS Subroutine (concluded) 
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TABLE 3-n 

DEFINITION OF ROSMOS SUBROUTINE VARIABLES 

CPRD = Concentration of solutes in product water under design (or test) 
conditions (dimensionless) 

DTH = System transient time step (hours) 

DTYME = Internal transient time step used (hours) 

LL = Maximum number of internal iterations of heat balance equations 

LPASS = Counter for number of iterations of heat balance equations 

POSM = Effective osmotic pressure of salts across membrane (psi) 

R51 = Variable used for temporary storage of R(5l) during thermal balance 

WPRD = Product mass flow rate under design (or test) conditions (Ib/hr) 

X = Y =, Multiplier used in heat balance equation; net heat input times X 
equals required temperature change (°F/Btu/hr) 
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3.4 SHOWER 

The SHOWER subroutine is designated as G-189A No. 67. 


3.4.1 Subroutine Description 


This subroutine models the thermal and evaporative mass exchange in a shower 
stall, as shown in Figure 3-9. The G-189A subroutine QSURR is used to model 


Ai r Water 

inlet inlet 

(primary) (secondary) 



Total outlet 
(primary) 


Figure 3-9. Shower Model Flow 
Schematic 


the thermal exchange between the 
shower stall frame and ambient 
environment. Evaporation is modeled 
by mass transfer equations as a 
function of air inlet flow rate, 
humidity and properties. The effect 
of the shower occur ant is included 
by standard metabolic equations 
based on input metabolic rate and 
respiratory quotient. 

3.4.2 Subroutine Data 

3.4.2. 1 General Notes 

a. Air inlet is on the primary side; 
water inlet is on the secondary. 
Total outlet flow is on the pri- 
mary side. . 


b. All gaseous flow codes are allowed (1, 2, or 3) on the primary side, and 
all liquid flow codes (0 or 4) on the secondary. 

c. Oxygen and COg balance will only be made for primary flow codes 2 and 3. 

d. The shower stall is considered to be occupied if-and only if there is flow 
inlet to the primary (air) side. 
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e. The thermal, metabolic oxygen, and CO 2 effects of the crewman are 
accounted for within the shower stall when occupied. Therefore, the 
occupant should not also be included elsewhere in the G-189A system 
model during the time the shower is occupied. 

f. When the water is turned on, the total (primary and secondary) water 
inlet and outlet flowrates are equal. When the water is turned off 
with air flow only, the outlet water vapor flow is equal to the inlet 
flow plus the amount of water evaporated. To preserve the overall 
water balance, the moisture retained in the towel and lost by door 
opening should be subtracted from the total water usage in GPOLY logic 
as required. 

g. To satisfy the total cabin heat balance, the heat transferred through 
the shower stall frame, R(53), should be added to the associated cabin 
sensible heat addition [cabin component location R(66)] in GPOLY logic 
each iteration. 


3. 4,2. 2 Heat Loss V-Array Data 


Reference 

Location 

Description 

Data Type 

51 

Effective average shower stall frame 
temperature (°F) 

1(0), 0. 

52 

Effective summed conductance between 
shower stall and ambient surroundings 
(Btu/hr-°F) 

0 

53 

Shower stall frame heat loss to 
ambient, excluding heat transfer to 
air and water flowing through stall 
(Btu/hr) 

0 

54 

Ambient gas temperature (°F) 

1(0) 

55 

Thermal convection conductor between 
shower stall external surface and 
ambient surroundings (Btu/hr-‘^F) 

1(0) 
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3. 4. 2. 2 (Continued) 


Reference 

Location 

56 

Description 

Convective heat loss from shower stall 
outer surface to ambient gas (Btu/hr) 

57 

Ambient wall temperature (°F) 

58 

Thermal radiation conductor (A^) between 
ambient walls and shower stall (sq ft) 

59 

Radiative heat loss from shower stall 
to ambient walls (Btu/hr) 

60 

Temperature of ambient structure 
connected to shower stall (°F) 

61 

Thermal conductor from shower stall to 
attached ambient structure (Btu/hr °F) 

62 

Thermal loss from shower stall to 
attached ambient structure (Btu/hr) 

63 

Same as R(51) on previous iteration 

64 

Thermal conductor through insulation 
around shower stall. No insulation 
is assumed, and this value must be 
set to zero. 


3. 4. 2. 3 Steady-State K-Array Data 

Reference ^ . 

Location Descn^tiM 

2^5 This location is used to store flow 

condition on previous iteration. 

0 = shower in use with water on; 

1 = shower in use with air flow 
only; -1 = shower not in use 


Data Type 

0 

1(0) 

0 

1(0) 

1 ( 0 ) 

0 

0 

I(R) 

Data Type 

0 
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3. 4. 2. 4 Steady-State and Transient V-Array Data (locations used for 
transient runs only are marked with an asterisk) 


Reference ^ ^ , 

' Location Description Data Type 

65 Not used 

66 Crewman metabolic rate when shower is 1(0) 
occupied (Btu/hr) 


67 CO 2 generation rate of shower 0 

occupant (Ib/hr) 

68 Respiratory quotient of shower 1(0) 

occupant (moles CO 2 produced/moles 

O 2 used) 

69 Shower stall net internal void volume 1(0) 

when occupied (cu ft) 


70 

71 


72 

73 

74 


Oxygen usage rate of shower occupant 0 
(Ib/hr) 

Thermal radiation conductor when 1(0) 

shower is occupied, shower wall *to- 
occupant (sq ft) 

CO 2 partial pressure in shower stall 0 
(mm Hg) 

Oxygen partial pressure in shower 0 

stall (psia) 

Effective thermal conductor (conduc- 1(0) 

tion and convection) between shower 
occupant and walls when occupied (Btu/hr °F) 


75 Water evaporation correction factor: 1(0) 

ratio of actual evaporation rate to 
predicted value from equation (3.4.18); 
(dimensionless ) 


76 Total electrical power (for lights) in 1(0) 

shower stall when in use (watts) 
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Reference 

Location Description 


77 

Shower occupant actual skin temperature 
(“F) 

0 

78 

Effective surface area of shower 
occupant skin for air drying (sq ft) 

1(0) 

79 

Effective heat transfer coefficient 
within shower, air-to-occupant/ 
walls (Btu/hr sq ft °F) 

1(0) 

80 

Diffusion coefficient for water vapor 
in air at 77°F and 14.7 psia (sq ft/ 
hr) 

1(0) 

81 

t 

Shower stall total inside wall surface 
area for evaluating thermal convection 
(sq ft) 

1(0) 

82 

Fraction of shower stall inside wall 
area that is wetted (for evaluating 
evaporative drying) 

1(0) 

83 

Water evaporation rate inside shower 
stall during drying (Ib/hr) 

0 

*84 

Shower stall frame thermal capacitance 
(Btu/°F) 

1(0) 

*85 

Mass of shower occupant when occupied 
(lb) 

1(0) 

*86 

Maximum transient computing time 
increment (seconds) 

1(0) 

*87 

Shower occupant initial skin temp (°F) 

1(0) 

88 

This location is used to store the 
outlet water vapor partial pressure on 
the previous iteration (psia) 

0 
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3.4.3 Analytical Model Description 
3.4.3. 1 Thermal Balance 

^The thermal nodal network assumed for the shower component is shown in 

Figure 3-10 for the occupied and unoccupied cases. T|ie shower is assumed 

to be occupied if and only if there is flow inlet to the primary (air) 

Water 



sk 


(a) Unoccupied 


(b) Occupied 


Figure 3-10. Thermal Model of Shower Stall Component 


3. 4. 3. 1.1 Unoccupied Case 


Considering the unoccupied case first, an effective ambient surrounding 
temperature T^ is first found. To do this, the radiation conductor is 
first linearized as follows: 



cjA 3 - 


/t 2 + T 





(3.4.1) 
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Where temperatures are all in absolute degrees Rankine. The effective 
s ambient temperature is then given by: 


T = 


^rad "^rad ^ ^conv ^conv ^ccnd "^cond 


G j + G +6 
rad conv cond 


(3.4.2) 


For the steady-state case, the shower wall temperature will be equal to 
the effective ambient temperature T . 


For a transient solution, the effective thermal conductance to the ambient 
temperature is computed as 


Sff ' ®rad ‘^cond ^ ®conv (3,4.3) 

Note that equation (3.4.3) requires that all three thermal conductors be 
connected directly to the shower wall node. This is a slight deviation 
from the standard usage of the G-189A subroutine QSURR, and requires that 
the insulation conductance [R-array location R(64)] be set to zero. This 
restriction presents no problem since no insulation is planned for use on 
the zero-gravity whole body shower currently designed for space use. The 
differential equation for heat transfer from the unoccupied shower stall 
may then be expressed as: 


dQ = C dT = - 

0) 0) 

®eff (T, - T„) dt 

(3.4.4) 

dT /G _ 

to , ( eff 

\ j _ ®eff 

(3.4.5) 

dt M C 

\ to 

i w C 

/ to 


The solution to this equation is: 


-(H‘ 

e “ j (3,4.6) 

REPROD U ClRlLi 1 i ■ OF THE 
ORXONAL page JR POOR 
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3. 4. 3. 1.1 (Continued) 

For the transient case, the unoccupied shower stall temperature is computed 
.from equations (3.4.1), (3.4.2), and (3.4.6). 

3. 4. 3. 1.2 Occupied Case 

When the shower is occupied (inlet air flow is present), the thermal model 
in Figure 3-10(b) applies. The basic thermal balance for the shower wall 
and occupant skin temperatures is the same whether the water is On or off. 
New temperatures are computed as follows: 

^new ” ^old ^ ^^*^into node (3.4.7) 

” '^e’‘^evap "*^surr''’‘^sen ’^*^c^'*'sk”'*'a)' ^^^r^"^sk ”'*'03 ^ (3.4.8) 

' w u 

I^sk “ ‘^MET''^evapg|^’^'^sengl^'’’^c^^co”'''sk^''’'^*^r^'^w^’'*'sk'^^ (3.4.9) 


For steady-state cases, the value used for y in equation (3.4.7) is 


0.4 


'steady state sum of conductors to ground 


(3.4.10) 


Internal steady-state iterations are made until the nodal temperature change 
computed is 0.005°F or until 25 iterations have been performed. 


For a transient solution, the value of y is given by 


_ At 

^transient in C 

P 


(3.4.11) 


where m is the nodal thermal mass. The computing time increment At 
is taken to be 0.4 times tie stability time increment where 



At 


thermal mass 

stability ~ sum of conductors to ground 


(3.4.12) 


3-84 


D2-11857N 


3. 4. 3. 1.2 (Continued) 

Initial runs with a 90-second computing time step had computational 
instability problems. This was overcome by further limiting the internal 
' computing time increment to a user input value in R(86). A default value 
of 10 seconds is set by the subroutine for R(86) and was found satisfactory 
for the shower cases considered. If the system time step is less than the 
stability time increment, the system time step is used directly for At. 

If it is greater, then internal iterations of equations (3.4.7) through 
(3.4.9) are repeated until the system time step is reached. 

Water turned on 


When the water is turned on, it is assumed that the outlet air/water 
mixture is saturated. The fluid is assumed to have intimate contact with 
both the shower stall walls and occupant equally. Thus, its outlet 
temperature T^ will be the arithmetic average temperature of walls T^ \ 
and skin T^j^. Therefore, for this case the thermal resistance represented \ 
by conductors and in Figure 3-10(b) is neglected. The sensible 

heat input from the fluid to the shower wall and occupant in equations (3.4.8) 
and (3.4.9) is thus given by 


’sen 


= %m_c 


0) 


P P, 


(T^ -T )+hmc (T -T ) 
, Pin “ ^ Ps ^-n “ 


(3.4.13) 


q = ^m„c„ (T 

Mk P Pp 


-T I )+^am^c^ (T 
Pin = Ps 


s. -’sk> 
in 


(3.4.14) 


The mixing of the air and water mixture is done using the standard G-189A 
subroutine HBALNC. This subroutine evaluates the amount of water evaporated 
to give saturation conditions. The latent cooling effect from this amount 
of evaporation is placed 70?^ on the occupant skin and 30% on the walls. 

This 70/30 split was determined from the relative temperature and surface 
area of the occupant and walls. 
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3. 4. 3. 1.2 (Continued) 

Water turned off 

“ When air flow is present with the water (secondary flow) off, evaporative 
drying is assumed to occur. The latent heat of evaporation applied to 
the walls and occupant skin is evaluated in Paragraph 3. 4. 3. 2. With the 
convection heat transfer coefficient input in R(79), the sensible heat 
input from the fluid to the shower wall and occupant in equations (3.4.8) 
and (3.4.9) is given by 

qc^„ = h A (T -T ) (3.4.15) 

^sen c u 0 0 ) ' ' 

0 ) 


where the wetted wall and skin areas are input. The outlet air temperature 
is computed as follows: 



T„ + h^ A T + h A , T , 
P Pp p.j„ C gj 0 ) c sk sk 

m C + h^ (A + A^. ) 
p Pp c ' 03 sk' 


(3.4.17) 


3. 4. 3. 2 Evaporative Drying 

When air flow is on and the water is turned off, the evaporative mass 
transfer coefficient is evaluated according to the following relation: 



( 3 . 4 . 18 ) 


( 3 . 4 . 19 ) 
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3.4. 3. 2 (Continued) 

The water vapor diffusion coefficient D is input for 14.7 psia and 77°F. 
'To adjust for the actual conditions, the following equation from Reference 5 
is used: 



(3.4.20) 


where T and P are in absolute units. The constant 0.43 in equation (3.4.18) 
was determined by correlating actual shower runs with experimental data. 

If further correlation is required, the multiplying factor Cj^^ was added, 
which is input in R(75). If Cj^j is not input, the subroutine assumes a 
default value of 1.0. The evaporation rate is then given by 


m 


hpA 


evap R 


H20 


*^^*^surface 
(^surface ^ ^60) 


H20, 




■46ory 


(3.4.21) 


and the latent cooling is given as 


^evap evap 


(3.4.22) 


3. 4. 3. 3 Crewman Metabolic Simulation 

The thermal, metabolic oxygen, and CO2 effects of the crewman are accounted 
for within the shower stall when occupied. Therefore, the occupant should 
not also be included elsewhere in the G-189A model during the time the 
shower is occupied. The shower is considered to be occupied if and only 
if there is inlet flow on the primary (air) side. The occupant metabolic 
rate and respiratory quotient are input (default values of 650 Btu/hr and 
0.9 will be used if these data are not input). As recommended by Refer- 
ence 7, one third of the crewman thermal mass is used as a skin thermal 
node. Only the thermal effects on this node are included, as the central 
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3.4.3. 3 (Continued 

core mass is assumed to remain at constant temperature. The oxygen con- 
sumption rate and carbon dioxide generation rate are computed from the 
relation!: from Reference 8: 



X 10-3 


0.1708 - 0.0123 


/ rQ - 0.707 
\ 0.293 


(3.4.23) 



(3.4.24) 


The carbon dioxide balance within the shower stall is made assuming the flow 
model shown below: 



Complete and immediate mixing is assumed within the shower enclosure. 
Therefore, the total outlet flow is identical in composition to the gas 
within the shower, or 



(3.2.25) 


A mass balance within the shower stall over a differential time increment 
dt i s 


\ 02 .^ ' 


( 3 . 4 . 26 ) 
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3. 4. 3. 3 (Continued) 


Combining equations (3.4.25) and (3.4.26) and rearranging; 


dm, 


€02 


m. 


stall 


out 


dt 


m. 


m 


C02 


m, 


stall 


stall 


C02 


+ 0 ) 


in 


C02 


( 3 . 4 . 27 ) 


For continuous operating conditions, It is assumed the CO^ inlet flow rate 
and generation rate (and the total outlet flow) remain constant over a time 
step. Since the amount of COj in the stall is relatively small compared to 
the total fluid mass, the fractional variation in total mass is assumed 
negligible compared with the variation of CO 2 . Therefore, the solution 


to equation (3.4.27) is as follows: 


mcQo (t=t +At) 

‘^^^stall 


= e 


CO 2 


•[ 


”“Zstall'‘'"°’ 


m-T 


^C02j 


3co2 - 


\ 02 . ^ “CO 
in 


stall 


m. 


'out 


0 


out 




m^ 


1 At 


stall, 


V (3.4.28) 




For a transient solution, equation (3.4.28) is used to compute the total 
mass of CO, in the shower stall at the end of the system time step at 

For the stLdy-state case, the solution at at=- applies; whereby, C 

For the steacy sc , , j „te is obtained 

is equal to Brna- the outlet 

from equation (3.4.25). The COj partial pressure within the st 

in either case assuming a perfect gas; whereby. 


m„, Rco ^ 

“^stall 


C02 


2.785 V 


(3.4.29) 


The preceding analysis for CO^ levels may also be made for the “VS®" ^ 

balance. The final result, giving the total mass of Qp 

Tg POOR 
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3. 4. 3. 3 (Continued) 

after a time step At, is given by 


^09 “ ^02 * "’02 " '"02 

°2sta11 ° L stall ° 


(t=t +At) = 3, 


(t=t«) - 009 e 




mj 

out 
nu 

*\ stall 




At 


) (3.4.30) 


3q2 - 


^^ 02 . " “02 
in 


stall V ni. 


out 




Again, for the steady-state case, the mass of oxygen in the shower stall 
is simply 3 q 2- The outlet oxygen flow rate in either case is given by 


m, 


02 


m, 


stal 1 


02 


out 


m. 


m. 


(3.4.31) 


stall 


out 


and the oxygen partial pressure within the stall is given by 



mo 2 


stall 


^02 


(T + 460) 


(3.4.32) 


During transient operation, the input control variable K(16) is used to 
determine when the shower is being first turned on after being unoccupied. 
If it is, the oxygen and CO 2 partial pressures within the stall are 
reinitialized at their values for the inlet air. 
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3. 4. 3. 4 Pressure Drop Considerations 

The shower may be operated either with or without inlet water flow. Average 
^ shower use time is reported to be 10 minutes, with an average of 2.5 liters 
(0.66 gal) of water used. For a water supply flow rate of 75 kg/hr (165 
Ib/hr), the water will be turned on for a cumulative total of 2 minutes 
per shower. This will occur during several water applications. A 

vacuum pick-up device is located in the bottom of the shower stall to collect 
the water and air flow. Clearly, two-phase flow will sometimes be occurring, 
with the pressure drop across the vacuum collection device fluctuating widely 
as slugs of water are picked up with air flow. The G-189A program should ^ 
not be used to' predict such fluctuations due to possible computational 
instability problems. Therefore, when water flow is being inlet to the 
shower, the pressure drop of shower components should not be made. 

The shower subroutine does not account for the time lag of liquid water 
flowing through the stall. Thus, when the inlet water is turned off, air 
flow only is assumed to occur, and pressure drop may be computed. The 
pressure drop predictions of Reference 19 were used to determine data for 
the shower components. The values given for the various shower components 
are shown in Table 3-13. To derive G-189A pressure drop model data from 
these figures, the method of Reference 20 was used with the following 
relation: 


paP 

oc 


• 2-cc 

cc m 


From the data presented in Reference 20, and assuming turbulent flow, a value 
of 0.2 for was assumed. Thus, neglecting the air viscosity difference 
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3. 4. 3. 4 (Continued) 


between actual conditions and those estimated in Table 3-1 3 > the shower 
' component pressures may be predicted in G-189A runs from the following 
equation: 


ap = 





1.8 


T' 

y (3.4.33) 


where subscript t refers to the data in Table 3-13 (or any future test 
data which should become available), and the unsubscripted variables refer 
to the changing shower conditions computed during a G-189A run. Equation 
(3,,4.33) represents G-189A standard pressure drop option (3) (see Table 3-1), 
and the value of C' and the exponent 1.8 must be input. Assuming the 
air density to be 0.075 Ib/cu ft, the values of C' in equation (3.4.33) 
for the shov/er components described in Table 3-13 are as follows: 


Shower 

Component 

Shower stall (including 'vacuum 
pick-up device, line, fittings) 

Air heater (with lines, flowmeter 
and fittings) 

Water/air separator-air side 
(with attached fittings) 


Value of C for 
6-189A Pressure 
Drop Option 3 

4.74 X 10"5 
3.79 X 10-5 
8.70X10-5 


These values for C , together with the exponent of 1.8 in equation (3.4.33) , 
may be input as component pressure drop model data using the G-189A 
option (3) in Table 3-1. 
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TABLE 3-12 

DEFAULT VALUES FOR SHOWER SUBROUTINE INPUT DATA 

I 


R-Array* Default 

Location Value 


51 

70 

54 

70 

55 

14 

57 

70 

58 1 

14 

60 

70 

61 

30 

66 

650 

68 

0.9 

69 

39 

71 

13 

74 

4.5 

75 

1 

76 

24 

77 

90 

78 

16.575 

79 

0.8 

80 

1 

81 

46 

82 • 

1 

84 

10 

85 

170 

86 

10 

87 

80 


*See Section 3.4.2 for definition of variables. 
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TABLE 3-13 


ESTIMATED PRESSURE DROP (AND RISE) FOR VARIOUS 
ZERO-GRAVITY WHOLE BODY SHOWER COMPONENTS 

SHOWER 

COMPONENT 

PRESSURE 
DROP (OR RISE) 

Shower stall (including vacuum pick-up 
device, line, fittings) 

345 N/m^ 
(7.20 psf) 

Air heater (with lines, flowmeter 
and fittings) 

276 N/m2 
(5.76 psf) 

Water/air separator (with attached 
fittings) 

786 N/m2 
(16.42 psf) 

Blower - pressure rise indicated 
(with attached lines, fittings) 

1407 N/m2 
(29.38 psf) 


AIR 

FLOW 

81.6 kg/hr 
(180 Ib/hr) 

81.6 kg/hr 
(180 Ib/hr) 

91.9 kg/hr 
(202.5 Ib/hr) 

91.9 kg/hr 
(202.5 Ib/hr) 
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TABLE 3-14 

DEFINITION OF SYMBOLS FOR SHOWER SUBROUTINE DESCRIPTION 

A = Area (sq ft) 

C = Thermal mass (Btu/°F) 

Cj. = Multiplying factor to adjust analytical mass transfer coefficient 
(dimensionless) 

Cp = Thermal specific heat (Btu/lb°F) 

D = Diffusion coefficient for water vapor in air (sq ft/hr) 

e = Naperian base logarithm, 2.7183... 

3 = Thermal radiation gray-body shape factor (dimensionless) 

G = Thermal conductor (Btu/hr-°F) 

6 . = R(61) = Effective thermal conductor between shower stall frame 

" and attached ambient structure (Btu/hr °F) 

G = R(55) = Effective thermal conductor between shower stall frame 
and ambient gas (Btu/hr °F) 

G = Effective thermal conductor between shower stall frame and 
^ ambient surroundings (Btu/hr °F) 

G . = Effective thermal radiation conductor between shower stall and 

^ ambient wall (Btu/hr °F) 

h »= Thermal convection heat transfer coefficient in shower (Btu/hr- 
^ sq ft-°F) » 

hp = Evaporative mass transfer coefficient (ft/hr) 

= Fluid thermal conductivity in shower (Btu/hr-ft-°F) 

a 

m = mass (lbs) 

m = Mass flow rate (Ib/hr) 

P = Absolute pressure (Ib/sq ft) 

Pco 2 ■ R(72) = CO 2 partial pressure in shower stall (mm Hg) 


3-95 


D2-II8571-I 

TABLE 3-14 (continued) 

DEFINITION OF SYMBOLS FOR SHOWER SUBROUTINE DESGRIPTION 

Pq 2 = R(73) = Oxygen partial pressure in shower stall (psia) 
Pr = Air Prandtl number in shower (dimensionless) 


q = 3.413 X R(76) = Total electrical power (for lights) in shower 
® stall (Btu/hr) 

q = Latent heat of evaporation (Btu/hr) 

’GVap 

^MET ~ R(55) = Shower occupant metabolic rate (Btu/hr) 

q = R(53) = Shower stall frame heat loss to ambient, excluding heat 

transfer to air and water flowing through stall (Btu/hr) 

R = Ideal gas constant (ft-lb/lb °R) 

t 

RQ = R(68) = Metabolic respiratory quotient of shower occupant 
(moles CO 2 produced/moles Og used) 

Sc = Schmidt number for shower air (dimensionless) 

T = Temperature (°F, except where otherwise noted) 

t = Mission time (hours) 

T j = R(60) = Temperature of ambient structure connected to shower 
frame (°F) 


T = R(54) = Ambient gas temperature surrounding shower enclosure (°F) 
conv 

00 • 

t^ = Mission time at beginning of current time step (hours) 

T . = R(57) = Effective average ambient wall temperature surrounding 

’^^“00 shower enclosure (°F unless otherwise specified) 

U 02 = Metabolic oxygen consumption rate of shower occupant (Ib/hr) 

V = Shower stall internal void volume (cu ft) 

e = Constant in metabolic balance equation (3.4.28 , 30) 

a = Stefan-Boltzmann constant (0.1714x10"^ Btu/hr-sq ft- °R‘^) 

At = Transient solution time increment (hours) 


3-96 


'nii-'TTNAL FA6rS IS iOOS 


D2-U8571-I 

TABLE 3-14 (concluded) 

DEFINITION OF SYMBOLS FOR SHOWER SUBROUTINE DESCRIPTION 

= Air viscosity in shower (Ib/hr-ft) 

= Air density in shower (Ib/cu ft) , 

= Multiplier to compute new temperatures, defined by equation (3.4.7) 
= Metabolic carbon dioxide generation rate of shower occupant (Ib/hr) 

Subscripts 

a = Shower air 

c,cond = Thermal conduction 

conv = Therma-1 convection 

C02 = Carbon dioxide 

evap = Evaporation 

H20 = Water vapor 

in = Inlet condition 

new = Condition at end of current time step or iteration 
0 = Shower air outlet condition 

old = Condition at beginning of current time step or iteration 

02 = Oxygen 

p = Primary side 

r.rad = Thermal radiation 

s = Secondary side 

sen = Sensible heat 

sk = Shower occupant skin node 

stall - Condition inside shower stall ^ 

T = Total fluid property 

oo = Ambient condition surrounding shower enclosure 
0 ) = Shower stall walls/frame node 


y 

p 

Y 

u 


C02 
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METABOLIC BALANCE 

EVALUATE 0 ? CONSUMPTION AND CO? GENERATION 
PERFORM 0 ? AND CO? BALANCE IN STALL. 
COMPUTE 0 ? AND CO? OUTLET FLOW RATES. 


ADD WATER TO AIR FLOW AS 
ENTRAINED LIQUID. 


Figure 3~11. Logic Flow Chart for SHOWER Subroutine 
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1. Logic Flow Chart for SHOWER Subroutine (continued) 
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Figure 3-11. Logic Flow Chart for SHOWER Subroutine (concluded) 
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TABLE 3-15 

definition of SHOWER SUBROUTINE VARIABLES 

Variable used to store inlet air specific heat (Btu/lb °F) 

Effective thermal capacitance of shower occupant skin node 
(Btu/hr °F) 

Effective thermal capacitance of shower wall node (Btu/hr °F) 

Water vapor diffusion coefficient in air inside shower stall 
(sq ft/hr) 

: System transient time step (hours) 

■ Internal transient computing time step (hours) 

: Water' evaporation rate on shower occupant (Ib/hr) 

= Water evaporation rate on shower walls (Ib/hr) 

= Logical variable set to TRUE for first internal iteration with 
water either on or off; otherwise FALSE 

= Linearized radiation thermal conductor from shower stall walls 
to ambient wall (when shower is unoccupied) or to occupant 
shower is occupied); (Btu/hr-°F) 

= Sum of thermal conductors to shower occupant (Btu/hr °F) 

= Thermal conductor from inlet air/water to shower occupant (Btu/hr F) 

= Sum of thermal conductors to shower stall (Btu/hr °F) 

= Thermal conductor from. inlet air/water to shower walls (Btu/hr F) 

= Net heat input from fluid to wall and skin nodes respectively 
(Btu/hr) 

~ Heat transfer coefficient; air to shower occupant and walls 
(Btu/hr-sq ft-°F) 

= Mass transfer coefficient for water evaporation inside shower 
(ft/hr) 

= Total enthalpy of outlet air/water mixture (Btu/.hr) 

= Counter for number of iterations to determine saturated conditions 
= Internal heat and mass balance iteration counter 
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Revised 10-16-75 


‘TABLE 3-15 (concluded)’ 

DEFINITION OF SHOWER SUBROUTINE VARIABLES 

NTYME = Maximum number of internal heat and mass balance iterations 
required 

PR = Air Prandtl number inside shower (dimensionless) 

PVAP - Water vapor partial pressure resulting from mass transfer equation 

(psia) 

PWIN = Water vapor partial pressure of inlet air (psia) 

PWINF = Effective water vapor pressure within shower stall (psia) 

PWOUT = Outlet water vapor pressure computed on previous iteration (psia) 

QSKIN = Latent heat of evaporation on shower occupant (Btu/hr) 

QWALL = Latent heat of evaporation on shower walls (Btu/hr) 

R1,R2,R6,R10,R12 = Variables used to store inlet air total flow, 'temperature, 

humidity, O 2 and CO 2 flows during current iteration 

R51,R77 = Variables used to store wall and skin node temperatures during 
heat balance 

RHOAA = Variable used to store inlet air density (Ib/cu ft) 

RR6,RR7 = Variables used to store component outlet water vapor and liquid 
flows from previous system iteration 

SC = Air Schmidt number inside shower (dimensionless) 

IWIN = Inlet air temperature (°F) 

TWINF = Effective air temperature within shower stall (°F) 

TWOUT = Air/water outlet temperature (°F) 

VISCAA = Variable used to store inlet air viscosity (Ib/hr ft) 

WC02 = Weight of CO 2 in shower stall (lb) 

WOXY = Weight of oxygen in shower stall (lb) 

X51,X77 = Heat balance multiplier for computing new wall and skin node 
temperatures, respectively, defined by equation (3.4.7) 

XKAA = Variable used to store inlet air thermal conductivity (Btu/hr ft °F) 
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3.5 WASDRY 

The WASDRY subroutine is designated as G-189A No. 70. 

s 

3.5.1 Subroutine Description . - 

This subroutine will simulate operation of the following appliances: 

0 Clothes washer 

0 Clothes dryer 

0 Clothes washer/dryer combination 

0 Dishwasher 
0 Dish dryer 

0 Dishwasher/dryer combination ' - 

0 Towel /cloth drying rack 

A flow. schematic of the WASDRY component is shown in Figure 3-12. Seven 
operational usage phases may be simulated: 

Phase 0 - Unit off 

Phase 1 - Wash water fill 

Phase 2 - Wash (circulate) 

Phase 3 - Spin dry - wash water out 

Phase 4 - Rinse water fill 
Phase 5 - Rinse (circulate) 

Phase 6 - Spin dry - rinse water out 
Phase 7 - Dry 

The routine will control the switching between phases based on input cycle 
schedules if requested by the user. Thermal exchange between the tub and 
frame and the ambient environment is modeled using the standard G-189A 
library routine QSURR. During the drying phase, the evaporation process 
is modeled in detail as a function of air inlet flow rate, humidity and 
properties, velocity within the tub, and water retention in the load. 

The subroutine models the thermal /mass exchange in the agitator or tub only, 
with the peripheral pumps, valves, accumulator, etc. to be simulated by 
standard G-189A component routines. 


D2-118571-1 


3.5.1 (Continued) 


Air in 




(a) Washer/dryer (b) Towel/cloth rack dryer 

Figure 3-12. WASDRY Component Flow Schematic 

3.5.2 Subroutine Data 


3. 5, 2.1 General Notes 

a. Water for the wash cycle is inlet and outlet on the component primary 
side. All liquid flow codes (0 and 4) may be used. Air for drying 
is inlet and outlet on the component secondary side. All gaseous 
flow codes (1, 2, or 3) may be used. 

b. If a towel/cloth rack dryer is being simulated, only phase 0 and 7 
may be used. 

c. To satisfy the total cabin heat balance, the heat lost through the 
washer structure, R(53), should be added to the associated cabin 


3-105 



1 


D2-IIS57N 


3.5.2. 1 (Continued) 

component sensible heat addition [cabin component location R(66)] 
in GPOLY logic each iteration. 

3.5. 2. 2 Instruction Options 

NSTR(l): Type of unit to be simulated 

= 0 Washer and dryer combination 

* 1 Washer only 

= 2 Dryer only 

= 3 Towel /washcloth drying rack 

NSTR(2): Steady-state criterion to be used for dryer 

0 Constant evaporation rate assumed equal to total water evaporated 
[derived from R(77), R(94) and R{78)] divided by total time allowed 
in R(79) for drying. Note: This value is limited internally to 

a drying rate corresponding to saturated outlet conditions. 

= 1 Assume saturated dryer outlet conditions 

=2 Use normal transient dryer evaporation equations, while holding 
the residual water loading constant 

NSTR(3): Maximum air absolute velocity in dryer tub 

= 0 Use input value in R(80) 

= 1 Use air inlet velocity based on flow rate and cross-sectional 
area in R(81) 

= 2 Compute velocity based on tub outer diameter and rotational 
speed in R(82) and R(83) 

= 3 Same as 2 above, except limited to value in 1 above 
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3. 5. 2. 2 (Continued) 


NSTR(4): 

Type of load 


= o' 

Clothes 


= 1 

Dishes 


NSTR(5): 

Cycle phase logic, request for transient run 


= 0 

Do not perform cycle phase logic. User will control the 
cycle phase in R(95) 

= 1 

Perform cycle phase logic. Cycle phase to be used each time 


step will be determined and output in R(95) 


3. 5. 2. 3 

Heat Loss V-Array Data 


Reference 

Location 

Description 

Data Type 

51 

Effective average unit frame 
temperature [or drying rack fixture 
temperature if NSTR(1)=3]; (°F) 

1(0), 0 

52 

Effective summed conductance from unit 

0 

frame (or drying rack fixture) to 
ambient surroundings (Btu/hr °F) 


53 

Total heat loss from unit frame (or 
drying rack fixture) to ambient sur- 
roundings (excluding heat transfer to 
tub and clothes/dishes load); (Btu/hr) 

0 

54 

Ambient gas temperature (°F) 

1(0) 

55 

Thermal convection conductor between 
unit frame external surface (or 
drying rack fixture) and ambient 
gas (Btu/hr °F) 

1(0) 

56 

Convective heat loss from unit frame 
(or drying rack fixture) to ambient 
gas (Btu/hr) 

0 

57 

Ambient wall temperature (°F) 

1(0) 


3-107 


3. 5. 2. 3 (Continued) 


D2-II8571-I 


Reference 


Location 

Description 

Data Type 

58 

Thermal radiation conductor between 

unit frame (or drying rack fixture) 
and ambient walls (sq ft) 

1(0) 

59 

Radiative thermal loss from unit 
frame (or drying rack fixture) to 
ambient walls (Btu/hr) 

0 

60 

Temperature of ambient structure 
attached to unit frame (or drying 
rack fixture); (°F) 

1(0) 

61 

Conductive thermal conductor between 
unit frame (or drying rack fixture) 
and attached ambient structure 
(Btu/hr °F) 

1(0) 

62 

Conductive thermal loss from unit 
frame (or drying rack fixture) to 
attached ambient structure (Btu/hr) 

0 

63 

External surface temperature of 
insulation (if any) around unit 
frame (°F). Same as R(51) if no 
insulation is used. 

0 

64 

Thermal conductor, through insulation 
around unit frame (Btu/hr °F). May 
be zero if no insulation is used, and 
should be zero for towel /washcloth 
drying rack fixture. 

I(R) 

3. 5. 2.4 

Steady-State and Transient K-Array Data 


Reference 

Location 

Description 

Data Type 

16 

Steady-state operational cycle phase 
to be used 

I(R) 

17 

Specifies next mission time the unit 
is to be turned on; may be either 1, 

1(0) if NSTR(5)=1, 0 


2, or 3, corresponding to being first 
turned on at mission time in R(99)j 
R(IOO), or R(lOl), respecti vely . 
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3. 5.2. 4 (Continued) 


R0^6 ^0 nc0 

Location Description Data Type 

18 Cycle phase used for previous iteration 0 

* 

3. 5. 2.5 Steady-State and Transient V-Array Data 

Note: (1) Locations used for transient runs only are marked 
with an asterisk. 

(2) All locations used for towel /cloth drying rack 
[NSTR(1)=3] are marked 


Location Description Data Type 


65— Heat input to component from drum 0 

motor (Btu/hr) 

66“ Clothes or dishes load, dry (or 1(0) 

weight of dry clothes on drying 
rack); Obs) 

67 — Weight of water in tub (or in clothes I(R) 

for drying rack); (lbs). This value 
used for steady-state analysis only. 

68 Drum motor input electrical power for i(0) 

continuous wash/rinse/dry (watts) 


69 Not used 

70 Drum motor input electrical power for 1(0) 

high-speed spin dry (watts) 

71 Not used 


72-- Tub and contents temperature (°F); 1(0), 0 

(or temperature of clothes on drying 
rack) 

73 Weight of detergent used per wash (lb I(R) 

detergent per 1b water). Mote: This 

item may only be used if the primary 
flow code is 4 since the detergent 
in the outlet wash water is put in 
the Special Flow No. 3. 

74 Not used 
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3. 5. 2. 5 

(Continued) 

Reference 
Locati on 

Description 

s 

*75 — 

Effective average specific heat of 
clothes or dishes load-dry (Btu/lb °F) 

76— 

Thermal conductor between tub and unit 
frame (or between towels/cloths and 
drying rack fixture); (Btu/hr °F) 

77 — 

Residual water retention after spin- 
dry (or initial water loading for 
towels/cloths on drying rack), used 
for clothes load only (lb water/lb 
dry clothes) 

78— 

Residual water retention required to 
terminate the drying process (lb 
water/lb dry clothes or dishes) 

79- 

Maximum time to be allowed for drying 
process (minutes). Note: For auto- 

matic control [NSTR(5)=1], the drying 
phase may be terminated sooner, 
depending on the amount of water 
retained in the load. 

80— 

Effective maximum air absolute 
velocity in drum during drying 
(ft/hr); or air velocity over towel/ 
cloth drying rack. 

81 

Total cross-sectional area of air 
inlet to dryer drum (sq in) 

82 

Diameter of dryer drum (inches) 

83 

Rotational speed of dryer drum (rpm) 

84 

Not used . 

85 

Not- used . . 
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Data Type 

1 ( 0 ) 

1 ( 0 ) 

1(0) if NSTR(4)=0 


1(0) if NSTR(1)?^1 


1(0). if NSTR(1)?«1 


1(0) if NSTR(1)5^1 
and NSTR(3)=0, 0 


1(0) if NSTR(1)5^1 
and NSTR(3)=1 or 3 

1(0) if NSTR(1)^1 
and NSTR(3)=2 or 3 

1(G) if NSTR(1)?^1 
and NSTR(3)=2 or 3 
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3. 5. 2. 5 (Continued) 


Reference 

Location Description 

86 Ratio of effective average air velocity 

in dryer drum relative to clothes or 
dishes load to maximum absolute air 
velocity in R(80)(dimensionless) 

87 — Effective exposed surface area of 

clothes or dishes load for drying 
(sq ft/ lb dry load) 


88- - Diffusion coefficient for water vapor 

in air at 77°F and 14.7 psia (sq ft/ 
hr) 

89- - Experimental correlation factor to 

adjust analytical evaporative drying 
rate (dimensionless) 


90 

91 

92- 

93- 
94 


Not used 

This location used to store outlet 
water vapor partial pressure on 
previous iteration (psia) 

Water evaporation rate during drying 
before adjusting for capillary action/ 
surface tension effects (Ib/hr) 

Actual water evaporation rate during 
drying (Ib/hr) 

Residual water loading after spin- 
dry, used for dishes .load only 
(lb water/lb dry dishes) 


*95-- Current cycle phase to be simulated 

during transient run 

*96-- Weight of water in tub and clothes/ 

dishes load during transient run 
(lbs) ; or water retained in towels/ 
cloth for rack dryer 

97 Time to be allowed for wash cycle 

(mi nutes ) 


Data Type 
1(0) if NSTR(1)?«1 

1(0) if NSTR(1)?«1 
1(0) if NSTR(1)?«1 
1 ( 0 ) 

0 

0 

0 

1(0) if NSTR(4)=1 

I (R) . Also output if 
NSTR(5)=1 

I(R). 0 

1(0) if NSTR(5)=1 
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3. 5. 2. 5 (Continued) 


Reference ^ 

Location Description 

98 Time to be allowed for rinse cycle 

(minutes) 

Note: The following three cycle start times 

are automatically repeated every 24 
hours. Zero start times for R(IOO) 
and R(lOl) are allowable and will 
be ignored. 


*99 First mission time (hrs) for unit to 

turn on 

*100 Second mission time (hrs) for unit to 

turn on 


*101 

102 

103 

*104- 


Third mission time (hrs) for unit to 
turn on 

Total quantity of water to be used 
during wash cycle (lb) 

Total quantity of water to be used 
during rinse cycle (lbs) 

Internal computing time increment used 
(hrs) 


*105— 

*106— 

*107- 

*108 


Maximum computing time increment 
allowed for computing evaporative 
drying (sec) 

Thermal mass of emipty tub only 
(Btu/°F); input zero for rack dryer 

Thermal mass of unit frame (Btu/°F); 
or of towel/cloth drying rack fixture 

This location is used to store the 
mission start time of the current 
cycle phase (sec), used for cycle 
phase determination logic 


3-112 


Data Type 
1(0) if NSTR(5)=1 

I(R) if NSTR(5)=1 
I(R) if NSTR(3)=1 
I(R) if NSTR(5)=1 
1(0) if NSTR(5)=1 
1(0) if NSTR(5)=1 
0 , 

1 ( 0 ) 

1 ( 0 ) 

1 ( 0 ) 

0 
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3. 5. 2. 5 (Continued) 


Reference 
, Location 


Description 


Data Type 


109— Critical moisture content for water 1(0) if NSTR(l)?* 

evaporation from clothing load (1b 
water/lb dry clothes). For clothes, 
this value represents the point where 
capillary action/surface tension 
effects begin to retard the evaporation 
from a free liquid surface. For 
dishes, this value represents the 
point where the initially complete 
wetted area begins to decrease. 


110 


Time to be allowed for spin-dry (min) 1(0) if NSTR(5)=1 


3.5.3 Analytical Model Description 


3. 5. 3.1 Thermal Balance 
Nodes having thermal mass 


The thermal model assumed for the WASDRY component is shown in Figure 3-13. 
The temperatures of the washer/dryer tub and frame (or towel rack and 
cloths) are computed as follows: 



Y m 


into node 


(3.5.1) 


For a washer/dryer component, the net heat into the tub and frame is: 


“ "'s S, 




<Tf - ft) 


° ®tf ^''"t ■ “’motor ■ PsutT 


For a towel /cloth drying rack, the net heat into the cloth and fixture is: 


H '^a \ ‘a^^^ 't j ^evap '"tf ^^f " ^t^ 


"" ®tf ^"^t “ ' “^surr 
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Mater Ai’ r 

inlet inlet 




(a) Washer/dryer 


(b) Towel/cloth rack dryer 


Figure 3-13. Thermal Model for WASDRY Component 
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3. 5.3.1 (Continued) 

The latent heat of evaporation Pgygp is zero except for the drying 
•> phase 7. It is evaluated during drying from convective mass transfer 
considerations discussed in Paragraph 3. 5. 3. 2. The total motor heat 
generated is assumed to be simply the electrical power for the current 
cycle phase which is input. The ambient heat loss to the surroundings 
*^surr computed using the standard G-189A subroutine QSURR which 
models the conduction, convection, and radiation thermal exchange with 
the surrounding environment. 

For the steady-state case, the value used for y in equation (3.5.1) is 

_ 0^5 

"'^steady state ” sum of conductors to ground 

Internal steady-state iterations are made until the nodal temperature change 
computed is less than O.Ol^F or until 20 iterations have been made. For 
a transient solution, y is computed as 

_ At 

''^transient ~ Cp 

where m Cp is simply the thermal mass of the node. The computing time 
increment At is taken to be 0.5 times the smaller time increment required 
for computational stability for the two nodes having thermal mass, where 

. _ thermal mass 

^^stability sum of conductors to ground 

Initial runs vyith a 90-second computing time step had computational 
instability problems during evaporative drying. This problem was overcome 
by further limiting the internal computing time increment to a user input 
value in R(105) during the drying phase. A default value of 10 seconds 
is set by the subroutine for R(105) and was found satisfactory for the 
clothes and dishes dryer cases considered. If the system time step is 
less than the stability time increment, then the system time step is 
used directly for At. If it is greater, then internal iterations of 
equation (3.5.1) are repeated until the system time step is reached. 
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3. 5. 3.1 (Continued) 

Fluid flow 

In the case of a washer/dryer, the fluid within the tub is assumed to be 
at the same temperature as the tub and contents. Thus, the air and water 
inlet flow conductors in Figure 3-13 are evaluated as Jti Cp. For a towel 
rack dryer, a convective coefficient is computed, The Reynolds number is 

evaluated as 

Re = (3.5.2) 

y 

The air velocity V is input. From the rack surface area A^, also input, 
an effective average flow length L is computed assuming a square surface: 

L = (3.5.3) 


Thus, the Reynolds number for flow over a plane surface is computed, with 
the fluid properties y and p evaluated from the inlet conditions. The 
convection heat transfer coefficient, assuming laminar flow (Re<5xl05), 
is then found from the relation: 


^3 K 33 

h = 0.664 Re Pr-'*'' 
C L 


(3.5.4) 


The air outlet temperature is then computed as 


m C T + 6. T 
s p_ a 


s in 


a ‘t 


(3.5.5) 


out 


m^ + G. 
s p_ 


where the convection conductor G is the coefficient h^, times the cloth 


surface area which is input. 
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3. 5. 3. 2 Evaporative Drying 

During the drying phase 7, the evaporative cooling is determined from mass 
transfer relations for water vapor in air. For a towel rack dryer, since 
the convection heat transfer coefficient is known from equation (3.5.4), 
the mass transfer coefficient may be given by 




(3.5.7) 


For a dryer tub, since the convection coefficient was not used directly, 
the form of the equation for mass transfer over wetted wall columns was 
used (References 5, 6): 


hp = 0,023 Cj^ (3.5.8) 

The Schmidt number Sc is defined in equation (3.5,7). An effective 
Reynolds number within the dryer tub is found as follows. First, the 
effective average air velocity V relative to the clothes or dishes load 
may be input, or it may be found from the drum inlet air flowrate and duct 
size or the drum rotational speed, according to input instruction NSTR(3). 
Then, for flow continuity, 



(3.5.9) 
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3. 5. 3. 2 (Continued) 


An effective hydraulic diameter of the air flov/ oath within the drum is 


s given by 



(3.5.10) 


The effective Reynolds number within the dryer drum may then be found from 
the relation 


Re 



(3.5.11) 


Note that in both the mass transfer equations (3.5.6) and (3.5.8), it is 
primarily the form of the equation that is most important. It relates the 
relative effects of variations in inlet flow properties, flow rates, and 
dryer design conditions. The actual magnitudes of the mass transfer co- 
efficients thus predicted must be correlated with experimental data for 
reliable accuracy.and the correlation factor was included for this 
purpose. In the absence of experimental data with which to correlate the 
analysis, the subroutine will assume a default value of 0.4 for which 

was determined from limited analyses of Space Station clothes and dishes 
dryer concepts. 

The water vapor diffusion coefficient D is input for 14.7 psia and 77°F. 
To adjust for actual dryer conditions, the following equation from 
Reference 5 is used (where T and P are in absolute units): 


D 


.^i.75 

~P~ 


(3.5.12) 


With the mass transfer coefficient hj^ known, the evaporation rate during 
the drying process is computed from the relation: 


m 


ho A 


evap R 


H20 


Tr 


“surface 

77^ + 460) 

surface ' 


(Tq -t- 460)' 


(3.5.13) 
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3. 5. 3. 2 (Continued) 

and the latent cooling is given by 

'’evap ^ "^evap 

Fv. poration rate adjustment for " dry" co.n^tioni 

When the clothes or dishes load is very "«et." evaporation ^ ^ 

free-standing water surface, and the drying rate may be accurately predicted 
from equation (3.5.13). However, as drying continues, a point ^ 
beyond which the actual evaporation rate is less than the computed value. 

For dish drying, the reason for this is that the exposed dish area become 

partially dried, and the effective wetted area is less ‘^^" 'ITirus’ed 
Lh the dishes completely wet. A critical moisture content (W^) is used 

to determine the point where some of the dishes surface first begins to ^ 
become dried. The value ot must be determined experimentally and 

subroutine input, although a reasonable default value 1s assumed if it is 
not. After the critical moisture content is reached, the wetted surface 
area, and hence the drying rate, is assumed to vary in direct proportion 
the amount of water left on the dishes according to the relation. 

(3.5.14) 


• — rn 

evap(yj<|j ^ (equation 3.5.13) 


(y 


For a clothes load, as the material becomes dried, capillary action 
surface tension effects retard the evaporation process and must be accoun e 

for. In early runs, a relatively recent theory assuming a submerge 

evaporative interface in a porous material was used (References 9. 10, ). 

However, it was found this theory did not apply to such thin materials a 
cloth fabric, consequently, a simplified approach was taken. Reference Id, 
in which the drying phase is divided into two distinct zones. With very 
wet material, a free-standing water surface is assumed and the evaporation 
rate computed from equation (3.5.13) is assumed valid. As drying continue 
a critical moisture content is reached and a "falling-rate peno 

begins. During this latter period, the evaporation rate predicted by 
equation (3.5.13) is adjusted by the same equation (3.5.14) given a ove 

for a (dishes loa(d. 
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3 . 5.3.3 Washer/Dryer Phase Switching Logic 

The user may choose any of seven phases in which the washer dryer will 
s operate, as described in Paragraph 3.5.1, by setting the input value of 
K(16) or R(95). However, if input instruction NSTR(5) is given a value 
of 1, the subroutine will perform its own phase switching logic for trans- 
ient runs. The program checks each time step to see if the current phase 

is complete. If it is, the cycle phase R(95) is set to the next phase. 

It is up to the user, then, to check the cycle phase R(95) in 6P0LY logic 
each time step and set the proper inlet flow conditions. For example, if 
it is in phase 1 (wash water fill), an inlet water flow should be set. The 
criterion used to terminate each phase is listed below: 

PHASE 0 (unit off): The mission time desired to turn unit on is input. 

Up to three times per day may be input in R(99, 100, 101). When the turn- 
on time is reached, either phase 2 (for washers) or phase 7 (for dryers 

only or towel /cloth drying rack) is selected. Each turn-on time is auto- 

matically repeated every 24 mission hours. 

PHASE 1 (wash water fill): The weight of water in the washer drum is 

compared with -Jhe maximum weight to be used for the wash cycle in R(102). 
When that point is reached, the cycle phase is set to 2. 

PHASE 2 (wash-circulate) : The beginning time of this cycle is stored. 

When the time allowed for washing [input in R(97)] is reached, cycle 
phase 3 is selected. ^ 

PHASE 3 (spin d rv-vjash''v)at’er out): The residual water left after spin-dry 

is input and also the total time allowed for spin-dry. These inputs deter 
mine the water outlet flow rate. When the water loading in the drum has 
decreased to the input residual level, cycle phase 4 is selected. 

PHASE 4 (rinse water fill): The weight of the water in the washer drum is 

compared with the maximum weight to be used for the rinse cycle in R(103). 
When that point is reached, the cycle phase is set to 5. 
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3. 5. 3. 3 (Continued) 

PHASE 5 (rinse-circulate): The beginning time of this cycle is stored. 

, When the time allowed for rinsing [input in R(98)] is reached, cycle phase 6 
is selected. 

PHASE 6 (spin dry-rinse water out): The residual water to be left after 
spin-dry is input and also the total time allowed for spin-dry. These 
inputs determine the water outlet flow rate. When the water loading in 
the drum has decreased to the input residual level, either cycle phase 0 
(for washer only) or 7 (for dryer only or washer/dryer) is selected'. 


PHASE 7 (dry): ‘Two criteria are used to terminate this cycle phase. First, 
the maximum time allov/ed is input in R(79). Secondly, the drying phase 
will be automatically terminated sooner if the residual water loading is 
less than the input value in R(78). When this point is reached, phase 0 
is first selected, and then the next turn-on time is checked as discussed 
previously. 
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3. 5. 3. 4 Pressure Drop Considerations 

A clothes or dishes washer or dryer have not yet been designed or built 
, for space application. Thus, component pressure drop data, which will 
depend on the specific hardware design, cannot yet be specified. The 
clothes dryer design could be quite different from a conventional house- 
hold dryer due to an apparent optimum drying time of 4 hours (from 
Reference 21). Thus, much lower, air flow rates may be required. The 
clothes and dishes washers do no't. actually impose a pressure drop on the 
spacecraft water system. This is tecause there is never a water flow from 
the water system directly through the washer subsystem. Considering the 
washer cycle phases described in Paragraph 3.5.1, water enters the washer 
from the water system during the wash and rinse water fill cycles only. 

At this time, there is no water flow out of the washer and hence no 
pressure drop through the washer subsystem. Water from the washer enters 
the water system only during the spin-dry phase. At that time, the washer 
water pump imposes only an outlet pressure on the water system. 
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TABLE 3-16 

DEFAULT VALUES FOR WASDRY SUBROUTINE INPUT DATA 

K(17)=l 


Default 

Value 


R-Array* 

Location 

Default 

Value 

R-Array 

Location 

Clothes washer 
NSTR(4)=0 

Dishwasher 

NSTR(4)=1 

. 51 

70 

66 

4 

15.2 

54 

70 

75 

0.2 

0.25 

55 

3 

78 

. 0 05 

0.0001 

57 

70 

79 . 

240 - 

90 

58 

1 

80 

3343 

10587 

60 

70 

87 

10 

0.45 

61 

4 

89 

0.4 

0.4 

68 

50 

96(Phases 

R(77)xR(66) 

R(94)xR(66) 



4 or 7) 


70 

150 

102 

55 

15 

72 

70 

103 

55 

15 

76 

12 

106 

4 

4 

77 

0.25 

107 

10 

10 

81 

8 

109 

0.2 

0.01 

82 

16 




83 

120 




86 

0.8 




88 

1 




94 

0.01974 




97 

30 




■ 98 

30 




105 

10 




no 

6 





*See Section 3.5.2 for definition of variables. 
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TABLE 3-17 

DEFINITION OF SYMBOLS FOR WASDRY SUBROUTINE DESCRIPTION 

•> A = Surface area (sq ft) 

C = Thermal mass (Btu/°F) 

C = Experimental correlation factor to adjust analytical mass transfer 
^ coefficient (dimensionless) 

Cp = Thermal specific heat (Btu/lb °F) 

D = Diffusion coefficient for water vapor in air (sq ft/hr) 

d = Effective dryer air flow path hydraulic diameter (ft) 

H 

G = Thermal conductor (Btu/hr °F) 

G = Thermal conductor between washer/dryer tub and frame nodes (or 
between towel/cloths and drying rack fixture); (Btu/hr-°F) 

h - Thermal convection heat transfer coefficient of air over dryer 

^ contents (Btu/hr-sq ft-°F) 

hp = Evaporative mass transfer coefficient (ft/hr) 

k = Air thermal conductivity (Btu/hr-ft-°F) 
a 

L = Effective average air flow path length through dryer (ft) 

m = Mass flow rate (Ib/hr) 

Pr = Prandtl number for dryer air (dimensionless) 
q = Heat transfer (Btu/hr) 

a = Latent heat of evaporation in dryer (Btu/hr) 

^evap 

q^OTOR “ Washer/dryer total motor heat dissipated (Btu/hr) 

a = Washer/dryer frame node net heat loss to ambient (Btu/hr) 

^surr 

= Ideal gas constant (ft-lb/lb-°R) 

? Effective air Reynolds number (dimensionless) 

= Schmidt number for dryer air (dimensionless) 

= Pressure (psi a) 


R 

Re 

Sc 
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TABLE 3-17 (concluded) 

DEFINITION OF SYMBOLS FOR WASDRY SUBROUTINE DESCRIPTION 


'T 

At 

V 

P 

u 

Y 

w 

w_ 


= Temperature (®F, except where otherwise noted) 

= Transient solution time increment (hours) 

= Effective air velocity relative to dryer contents (ft/hr) 

= Dryer air density (Ib/cu ft) 

= Dryer air viscosity (Ib/hr-ft) 

= Multiplier to compute new temperatures, defined by equation (3.5.1); 
(dimensionless) 

= Moisture content in washer/dryer load (lb v;ater/lb dry clothes or 
dishes) 

= Critical moisture content below which evaporation rate is less than 
predicted rate for free-standing water surface 


Subscripts 
a = Dryer air 

f = Washer/dryer frame node (or towel drying rack node) 

H20 = Water 

in = Inlet condition 

new = Condition at end of current time step or iteration 

0 = Averaqe dryer air condition away from cloth or dishes surface 

old = Condition at beginning of current time step or iteration 


p = Primary side (water) 

r - Drying rack 

s = Secondary side (air) 

t “ Washer/dryer tub and contents node (or towels/cloths node 

for rack dryer) 

0 ) = Water 
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‘STORE TEMPERATURES AND DRYER AIR FLOW RATES. 

R20*R(20) evaluate* 

R25»R(25) EVALUATE. 

“IN 

T. 



DETERMINE TRANSIENT COMPUTING TIME 

STEP (DTYME) AND NUMBER OF ITERATIONS (NTYME) 


LPASSnO 


LPASS=LPASS+1 


NPHASE:7 


BEGIN EVAPORATIVE DRYING CALCULATIONS. 
EVALUATE P , T 
Wq ° 


SET constant” STEADY STATE 
EVAPORATION RAuE TO INPUT VALUE r 


W 



EVALUATE DRYER 
AIR PROPERTIES 
-VELOCITY 

-DIFFUSION COEFFICIENT 
-SCHMIDT NUMBER 


Figure 3-14. Logic Flow Chart for WASDRY Subroutine (Continued) 












STEADY 


R(96)=R(96)-m^At 


R(96):0 


m^=R(96)/At 

R(96)=0 


REPRODUCIBILITY OP THE 
ORIG3LJAL PAGE IB POOR 


EVALUATE FLOW AND THERMAL INPUTS 
FROM EVAPORATION. PROCESS. 
*^EVAP lily 
^MOTOR"^'^^^ R(68) 

R{25) =R25+nij^ , 

R(20) =R(24HR ( 25 ) 


Fiqure 3-14. Logic Flovi Chart for WASDRY Subroutine (Continued) 
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Figure 3-14. Logic Flow Chart for WASDRY Subroutine (Concluded) 
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TABLE 3-18 

DEFINITION OF WASDRY SUBROUTINE VARIABLES 

C = Combined thermal mass of washer drum, water and load (Btu/®F) 

CSTART = Start time of current cycle phase (seconds) 

DH = Effective hydraulic diameter of air flow path inside dryer drum (ft) 
DIFF = Water vapor diffusion coefficient in air (sq ft/hr) 

DTH = System transient time step (hours) 

DTYME = Internal transient computing time step used (hours) 

EVAP = Water evaporation rate (Ib/hr) 

GAIR = Thermal conductor from washer load to air inlet (Btu/hr-°F) 

HC = Heat transfer coefficient of air over clothes on drying rack 
(Btu/hr-sq ft-°F) 

HD = Mass transfer coefficient for water evaporation (ft/hr) 

HMIX = Total enthalpy of dryer inlet air (Btu/hr) 

KPASS = Iteration counter for determining saturated conditions 

LPASS = Internal heat and mass balance iteration counter 

LPHASE = Variable used to select washer phase for next iteration 

NPHASE = Operational cycle phase assumed (0 through 7) 

NTYME = Number of internal transient computing iterations 

PR = Air Prandtl number (dimensionless) 

PVAP = Computed outlet water vapor pressure (psia) 

PWIN = Inlet air water vapor pressure (psia) 

PWIMF = Effective water vapor pressure in drum (or over drying rack); (psia) 

PWOUT = Outlet water vapor pressure on previous iteration (psia) 
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QEVAP = 
QMOTOR = 
R20,R25 

RE 

RESID = 
SC 

TTUB = 

TWIN = 
TWINE = 
TWOUT = 
VMAX = 
VREL = 
W 


•TABLE 3-18 (concluded) 

DEFINITION OF WASDRY SUBROUTINE VARIABLES 

Latent heat of evaporation (Btu/hr) 

Heat added to structure from drum motor (Btu/hr) 

= Variables used to store value of R(20) and R(25) at beginning 
of system iteration 

Air Reynolds number (dimensionless) 

Residual water retention required to terminate drying process 
(lb water/lb dry clothes or dishes) 

Air Schmidt number (dimensionless) 

Variable used to store washer load temperature prior to making 
thermal balance (°F) 

Inlet air temperature (°F) 

Effective air temperature in drum (or over drying rack); (°F) 
Outlet air temperature on previous iteration (°F) 

Maximum absolute velocity of air over clothes or dishes (ft/hr) 
Air velocity relative to clothes or dishes (ft/hr) 

Water retention in load (lb water/lb dry clothes or dishes) 
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3.6 WASTEC 

The WASTEC subroutine is designated as 6-189A No. 68. 

3.6.1 Subroutine Description 

This subroutine will simulate a urine/fecal waste collector applicable to 
space use, such as a urinal or dryjohn. A flow schematic of the WASTEC 
component is shown in Figure 3-15. Three operational usage phases may be 
simulated: 

Phase 0 - Unit off 

Phase 1 - Urine collection 

Phase 2 - Fecal collection 

Phase 3 - Combined urine/fecal collection 


Air/ 

Air flush water 



Collector Urinal 

outlet outlet 


Figure 3-15. WASTEC Component Flow Schematic 

Commode contents may be under vacuum, if requested, during phases 0 and 
1. Thermal exchange between the collector and ambient environment is 
modeled using the standard G-189A library routine QSURR. The subroutine 
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3.6.1 (Continued) 

models the thermal/mass exchange in the collector and urinal only, with 

the peripht-al pumps, valves, etc. to be simulated by standard G-189A 

component routines. 

3.6.2 Subroutine Data 

3. 6. 2.1 General Notes 

a. The flow^ paths for the collector and urinal are isolated. The collector 
is on the component primary side, and the urinal is on the secondary. 

b. All gaseous flow codes (1,2, and 3) are permitted for both the collector 
and urinal.. For a urinal flush, the water must be inlet as entrained 
liquid in a gaseous fluid. 

c. Vacuum exhaust and drying of the collector are optional according to 
input instruction NSTR(l). If NSTR(1)=0, a collector vacuum is auto- 
matically assumed during phases 0 and 1. Each time vacuum is 
initiated following operation of a fecal collection phase (2 and 3), 

a pumpdown or blowdown is automatically performed. The time required for 
this is input. The vacuum pump, if used, is not a part of the WASTEC component. 

d. To satisfy the total cabin heat balance, the heat lost through the 
waste collector structure, R(53), should be added to the associated 
cabin component sensible heat addition [cabin component location R(66)] 
in GPOLY logic each iteration. 

3. 6. 2. 2 Instruction Options 

NSTR(l): Controls vacuum drying of commode contents 

= 0 Commode contents are assumed to be under vacuum continually 

except during fecal collection (phases 2 and 3). Vacuum pumpdown 
is activated automatically (during transient run only) immedi- 
ately after fecal collection or at beginning of run. 

= 1 Commode contents are assumed at ambient pressure. 
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J.6.2.3 Heat Loss V-Array Data 

Reference 

Location Description 

51 Collector contents temperature (°F) 

52 Summed thermal conductance between 
collector contents and ambient, 
excluding conductance to air flow 
through collector or urinal 
(Btu/hr °F) 

53 Collector contents net heat loss to 
ambient, excluding heat transfer to 
air flow through collector and urinal 
(Btu/hr) 

54 Ambient gas temperature (°F) 

55 • Thermal convection conductor from 

collector outer surface to ambient 
gas (Btu/hr °F) 

56 Convective heat loss from collector 
outer surface to ambient gas (Btu/hr) 

57 Ambient wall temperature (®F) 

58 Thermal radiation conductor(A3-) between 

ambient v;alls and collector outer 
surface (sq ft) 

59 Radiative heat loss from collector 
outer surface to ambient walls (Btu/hr) 

60 Temperature of ambient structure 
connected to collector (°F) 

61 Thermal conductor from collector node 
to attached ambient structure 
(Btu/hr °F) 

62 Thermal loss from collector to 
attached ambient structure (Btu/hr) 

63 Collector outer surface temperature 
(°F) 

64 Thermal conductor between collector 
contents and collector outer surface 
(Btu/hr °F) 
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Data Type 
1 ( 0 ). 0 
0 


0 

1(0) 

1(0) 

0 

1 ( 0 ) 

1 ( 0 ) 

0 

1(0) 

1 ( 0 ) 

0 

0 

1(0) 
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3. 6. 2. 4 

Steady-State V-Array Data 


Reference 


Location 

Description 


65 

Net heat generated in unit by.slinger 
motor (Btu/hr) 

0 

66 

Operational usage phase (0 through 
4 as described previously) 

I(R) 

67 

Slinger motor rated electrical power 
(watts) 

1(0) 

68 

Not used 


69 

Urine flow rate during micturition 
(Ib/hr) 

1(0) 

.70. ■ ' • 

t 

Not used 


71 

Not used 


72 • 

Not used 


73 

Not used 


74 

Exposed surface area of collector 
contents, used for evaporative drying 

1(0) 


calculations (sq ft) 


75 

Evaporation coefficient f in equation 
(3.6.6) for vacuum drying (dimension- 
less) 

I(R) 

76 

Evaporation or sublimation rate of 
collector contents computed from 
equation (3.6.S); (Ib/hr) 

0 

77 

Effective thermal conductor between 
air flow through collector and 
collector contents (Btu/hr °F) 

1(0) 

78 

Effective thermal conductor between 

1(0) 

flow through urinal and collector 
contents (Btu/hr °F) 


79 

Urine and feces inlet temperature 
(°F) 

1(0) 


Data Type 


if NSTR(1)=0 
if NSTR(1)=0 
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3.6. 2. 5 

Transient K-Array Data 



Reference 

Location 

Description 


Data Type 

16 

Operational cycle phase used on 
previous time step 

0 


17 

This variable is used in certain 
program logic to indicate pressure 
inside commode during the previous 
time step, as follows: 

0 



- 1 ambient pressure 
0 ■> hard vacuum 
5 -> vacuum pumpdown or 
blowdown applied 



3. 6.2. 6 

Transient V-Array Data 



Reference 

Location 

Description 


Data Type 

80 

Time required for vacuum pumpdown 
or blowdown (min) 

1(0) 

if NSTR(1)=0 

81 

Internal empty volume of fecal 
collector (cu ft) 

1(0) 

if NSTR(1)=0 

82 

Total cumulative mass of urine 
collected during current usage 
phase (lb) 

0 


83 

Total urine generated per crewman 
Ob/use) 

1(0) 


84 

Total mass of free gas in collector 
to be evacuated during vacuum pump- 
down (lb) 

0 


85 

Effective thermal mass of empty 
collector node (Btu/°F) 

1(0) 


86 

Total fecal mass added to collector 
per usage (lb) 

1(0) 


87 

Inlet fecal moisture content (lb H 2 O/ 
lb wet fecal material) 

1(0) 


88 

Specific heat of dried feces (Btu/lb ° 

F)I(O) 
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3. 6. 2. 6 

(Continued) 


Reference 

Location 

Description 

Data Type 

89 

Moisture present in collector in 
liquid phase (lb) 

1(0), 0 

90 

Moisture present in collector in 
solid phase (lb) 

1(0). 0 

91 

Total dry mass of collector contents 
(lb) 

1(0), 0 

92 

Not used 


93 

Start time of current drying phase 
(seconds) 

0 

94-98 

Not used 


3. 6. 2. 7 

Extra V-Array Data (must be specified for 

a transient case) 

Reference 

Location Description 

Data Type 

99 

Reserved to store density (Ib/cu ft) 
of upstream primary source. Used for 
initial condition for vacuum pumpdown. 

0 

lOO^x 

Reserved to store' R-array flow data 
for inlet primary flow on previous 
iteration. Used to store initial flow 
and constituents for vacuum pumpdown. 
Reserve 4+5*NPF locations. 

0 

8.A.3 Analytical Model Description 



3. 6. 3.1 Thermal Balance 

The thermal model assumed for the WASTEC component is shown in Figure 3-16 
The collector contents temperature is computed by a standard forward- 
difference method: 


= T. 


'new 


'old 


^'^into node 


(3.6.1 
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3.6.3. 1 (Continued) 

^^into node " (^^out’^^) 


The total slinger motor heat 
generated is assumed to be simply 
the electrical input power for 
the phase being simulated, The 
latent heat of evaporation or 

sublimation ^evap 
except when the collector is 
being vacuum dried. During 
vacuum drying, the latent heat is 
computed as discussed in Para- 
graph 3. 6. 3. 2. The ambient heat 
loss to the surroundings 
is computed using the standard 
G-189A subroutine QSURR which 
models the conduction, convection 
and radiation thermal exchange 
with the surrounding environment. 

For the steady-state case, the 
value used for r in equation 
(3.6.1) is 


■^^M0T0R“^surr"'’evap (3.6.2) 


Collector Urinal 

inlet inlet 



Figure 3-16. Thermal Model for 
WASTEC Component 


_ 0.5 

^steady state ~ sum of conductors to ground 

Internal steady-state iterations are made until the nodal temperature change 
computed is less than 0.05°F or until 20 iterations have been made. For a 
transient solution, y is computed as 

At 

^transient C 

c 
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3. 6. 3.1 (Continued) 

The computing time increment At is taken to be 0.4 times the time increment 
s required for computational stability, where 


^^stability ~ sum of conductors to ground 


If thq system time step is less than the stability time increment, the 
system time step is used directly for At. If it is greater, then internal 
iterations of equation (3.6.1) are repeated until the system, time -step is 
reached. After computing the collector contents- temperature - Jj,, the fluid 
outlet temperatures are calculated as follows: 


out 


%'p Tp. 

^ *^p ^in 


(3.6.3) 





(3.6.4) 


3. 6. 3. 2 Vacuum Drying 

During vacuum drying of the collector, the evaporation or sublimation rate 
is assumed to be governed by the Hertz-Knudsen relation (References 13, 14): 



2.784 a 


( 


surface 



RT. 


(3.6,5) 


where pressures are in mm Hg and the temperature is in °R. The accommodation 
coefficient a has a theoretical maximum value of 1, but varies widely with 
the purity of the liquid being evaporated and the substrate material (Ref- 
erences 15, 16 5 17). A discussion of the applicauion of equation (3.6.5) 
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3. 6. 3. 2 (Continued) 

to freeze-drying of foods is given in Reference 18. In standard freeze- 
, drying practice (Reference 13), the gas pressure above the evaporating 
or sublimating surface is assumed negligible compared with the water vapor 
partial pressure ^“surface’ water loss rate is given by: 


m 

(ji 


f A, P 

“ ^ “surface 


V 


M 

ijii 

2 IT RTq 


Since essentially only water is being evaporated in the collector, the 
constants M , tt, R, and f were combined into one constant f, and 
the following equation 


m 

U3 


^ “surface 


(3.6.6) 


was used in the subroutine to compute the evaporation or sublimation rate 
of the moisture in the collector contents, where 

The evaporation coefficient f is a subroutine input and may be varied to 
correlate with experimental data. A default value for f of 3.0 is used 
in the subroutine, which was found to agree best with limited test results. 
The outlet flow rate is assumed to be saturated water vapor, and the G-189A 
subroutine PSAT is used to compute the outlet gas pressure. The latent 
cooling is given by: 


^evap 


1042 m (evaporation, T ^ 32°F) 

U} 

1225 m (sublimation, T_ ^ 32°F)‘ 
0) c 


(3.6.7) 


During a transient solution, the water/ice balance is kept track of. When 
passing through the freezing point, the collector contents temperature is 
held constant at 32°F until all the water has either melted or frozen. 
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3. 6. 3. 3 Mass Flow Additions 

During phases 1 (urine collection), 2 (fecal collection), and 3 (urine/ 
fecal collection), mass is added to the dryjohn and/or air flow. During 
phases 1 and 3, the urine is added as entrained water in the inlet air flow 
according to the micturition flow rate which is input. The total void mass 
per man is also input, and when the total urine flow has reached this amount 
it is terminated for the current phase. It will automatically begin again 
when phase 1 or 3 is reinitiated. Each time phase 2 or 3 is initiated, the 
collector contents mass and water contents are increased by the total defeca 
tion mass per crewman, which is input. 

3. 6. 3.4 Pressure Drop Considerations 

The current dryjohn subsystem, which will be flown on Shuttle, is still in 
the design stage; and pressure drop data for the various components are not 
yet available. Consequently, the G-189A pressure drop model input data 
cannot be specified at this time. However, when the data become available, 
the method described previously for the shower data [depicted by equation 
(3.4.33)] may be used. One possible exception to this is during urine 
collection, when there can be two-phase flow of urine mixed with air. For 
this condition the pressure between the urinal and water separator may 
fluctuate, and an attempt to include this in the model could cause compu- 
tational instability problems. 
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TABLE 3-19 

DEFAULT VALUES FOR WASTEC SUBROUTINE INPUT DATA 


R-Array* 
Locati on 


Default 
Va1 ue 


51 

60 

54 

70 

55 

2 

57 

70 

58 

0.8 

60 

70 

61 

1 

64 

2.5 

67 

20 

69 

108 

74 

1 

75 

3 

77 

1.5 

78 

0.15 

79 

97 

80 

3 

81 

1.2 

83 

0.6 

85 ‘ 

0.5 

86 

0.3 

87 

0.75 

88 

0.25 

89 + 90 

.08 

91 

2.5 


*See Section 3.6.2 for definition of variables. 
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TABLE 3-20 

DEFINITION OF SYMBOLS FOR WASTED SUBROUTINE DESCRIPTION 
A = Exposed surface area of collector contents (sq ft) 

^ ' 

C = Nodal thermal mass (Btu/°F) 

C = Thermal specific heat (Btu/lb°F) 

f = R(75) = Modified vacuum evaporation coefficient used in subroutine, 
equation (3.6.6) 

f^ = Evaporation/sublimation coefficient for vacuum drying (dimensionless) 

G = R(77) = Effective thermal conductor between air flow through collector 
^ and collector contents (Btu/hr °F) 

G = R(78) = Effective thermal conductor between flow through urinal and 
^ collector contents (Btu/hr °F) 

m • = Mass flow rate (Ib/hr) 

M = Molecular weight 

= Water vapor partial pressure (psia unless otherwise noted) 

P = Gas pressure above collector surface (mm Hg) 

oo 

q = Heat transfer rate (Btu/hr) 

^evap ~ Latent heat of evaporation or sublimation in collector (Btu/hr) 

^MOTOR ” 'input to collector from motor (Btu/hr) 

q „ = Collector node net heatj, loss to ambient (Btu/hr) 

^surr 

R = Universal gas constant (-ft lb/mole-°R) 

T = Temperature (°F unless otherwise noted) 

At = Transient solution time increment (hours) 

Y = Multiplier to compute new temperatures, defined by equation (3.6.1); 
(dimensionless) 

a = Accommodation coefficient for vacuum evaporation/sublimation rate, 
equation (3.6.5); (dimensionless) 

1 / 


3-144 


D2-II8571-I 


TABLE 3-20 (concluded) 

DEFINITION OF SYMBOLS FOR WASTEC SUBROUTINE DESCRIPTION 

s Subscripts 

c = Collector 

in = Inlet flow condition 

new = Condition at end of current time step or iteration 
old = Condition at beginning of current time step or iteration 
out = Outlet flow condition 
p = Primary side (collector) 

s = Secondary side (urinal) 

u = Urinal 

to = water 
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CALL QSURR 


EVALUATE 
MASS OF AIR 
IN COLLECTOR 


LPASS=1 


: STEDDY : 


DETERMINE STABILITY TIME 
INCREMENT (DTH) AND 
NUMBER OF THERMAL PASSES 
(LPASS) REQUIRED 


p, L1.L=LLL+1 


TIME SPENT 
CURRENT 


PHASE: VACUUM 
PUMP-DOWN TIME 


PERFORM VACUUM 


PUMPDOWN. 


COMPUTE OUTLET 

, ■ 

PRESSURE, FLOW RATE 


NPHASE:2 


EVALUATE VACUUM DRYING 

m =fA,P , /ypn 

^ “surface ^ 


NPHASE:3 

NSTR(1):1 

STEDDY 

K(2):0 


STEDDY 


EVALUATE q^vap'^R 
STEADY-STATE CASE. 
COLLECTOR LIQUID IS 
ALL EITHER WATER 
OR ICE. 


EVALUATE^^ FOR 

TRANSIENT CASE 
BASED ON WATER/ ICE 
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PRESSURE, TEMPERATURE 
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Figure 3-17. Logic Flow Chart for WASTEC Subroutine (Continued) 
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TABLE 3-21 

DEFINITION OF WASTEC SUBROUTINE VARIABLES 

CAP = Thermal capacitance of collector and contents (Btu/°F) 

CPBB = Variable used to store value of inlet air property CPB 

DTMIN = Time elapsed from start of vacuum pumpdown (minutes) 

DTYME = System transient time step (hours) 

HMIX T Total enthalpy of inlet urinal air mixture (Btu/hr) 

LLL = Internal heat balance iteration counter 

LPASS = Number of heat balance iterations required 

NPHASE = Operational usage phase assumed (0 through 3) 

QEVAP = Latent heat of evaporation or sublimation of collector contents 
moisture (Btu/hr) 

Q.IN- = Net heat transfer into collector and contents (Btu/hr) 

QMOTOR = Heat added to collector from slinger motor (Btu/hr) 

RHOBB = Variable used to store value of inlet air property RHOB 
SUBEVA = Vacuum evaporation or sublimation rate (Ib/hr) 

UFLOW = Urine flow rate addition (Ib/hr) 
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4.0 SUBROUTINE MODEL VERIFICATION RESULTS 

Each; of the appliance subroutines was initially checked for accuracy by 
modeling a single G-189A component and running it separately (i.e., not 

\ 

connected in an all-up system model) with dummy ambient and inlet flow 
conditions. These conditions were held constant during a run and a steady- 
state and transient solution obtained. The results of these checkout runs 
are described in the following paragraphs for each appliance subroutine. 
Calculations are presented to verify conservation of energy within the 
components, and the solutions are compared with test data where available. 
Many of these appliances have not yet been built for spacecraft application, 
and available test data required for correlation of the models are limited. 
For such cases, additional testing and model correlation are recommended. 

The. inlet flow conditions tabulated in the following paragraphs for each 
component consist of the standard G-189A A- and B-arrays. The A-array 
contains the primary side inlet flow rates, temperatures, and pressures; 
the B-array contains the same data for the secondary side. The definition 
of these arrays is given in Table 4-1, which is taken from the G-189A 
Program Manual (Reference 2). In addition, the R-array locations 1 through 
38, which have the same definition for all components, are described in 
Table 4-1. These R-array locations comprise the outlet flow conditions 
for the component solutions given in the following paragraphs. For a more 
complete description of these arrays, the G-189A Program Manual (Reference 
2) should be consulted. The rest of the component solutions, as well as 
their model input data, are stored in the remainder of the R-array; and 
these R-array locations for each component are described in the subroutine 
descriptions in Section 3. 




TABLE 4-1 
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r-xj 
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O 
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-1 

o 

9 
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G-189A COMPONENT INLET AND OUTLET FLOW DATA DEFINITION 




SaeoiidasT 

Prinary ' 

Secondary, 



Source 

Source 

Side 

Side 

. 


Flov 

Flov 

Coaponeat 

Component 



Data 

Ite.ta 

Data 

Data 


Ie 

A(l) 

B(l) 

R(l) 

R(20) 

Total Plow (1-b/hr)' 


A(2) 

E(2) • 

R(2) 

R(21) 

Fltd.d Teaiperat’jre (®P) 

? ! 

A(3) 

B(3) 

R(3) 

R(22) ■ 

Upstream Duct Outlet Pressure (psla) 

4-2 

A(i4) 

B(L) 

. R(4) 

R(23) 

Component Outlet Pressure (psla) : 

i' 
f ^ 

A(5) 

B(5) ' 

R(5) 

R(2U) 

Non-CondensaBles Flow (it /hr) 


A(6) 

B(6) • 

R(6) 

R(25) 

Condensable Vapor Flow (ih/hr) f 

1 i 

A(7) 

B(7) 

R(7) 

R(26) 

Condensable Entrained Liquid Flov (ITj/hr) 

i i 

A(8) 

B(8) 

R(8) 

R{27) 

Hon-Condensahles Specific Heat (Btu/lb-®F) 

II 

A(9) 

B(9) 

R(9) 

R(28) 

Ron-Condensables Molecular Weight (ib/nol) 

ll 

A(IO) 

B(10) 

R(10) 

R(29) 

Oxygen Flow (ib/hr) 

r 

A(ll) 

B(ll) • 

R(ll) 

R(30) 

Diluent Flow (lb /hr) 

ll 

A{12) 

3(12) 

R(12) '• 

R(31) 

■ Carbon Dioxide Flov (ib/hr) 

■. 5 

A(13) 

B(13) 

R(13) • 

R(32) • 

Trace Contaminants Flov (ib/hr) 

' i 

A(lL) 

B(lL) 

• r(iA) 

R(33) 

Special Flov No. 1 (ib/hr) 

, a 

A(15) 

B(15) 

R(15) 

R(3J+) 

, Special Flov No. 2 (lb /hr) 


A{16) 

B(l6) 

R(l6) 

R(35) 

Special Flov No. 3 ( Ib/hr) • 

\ 

A(17) 

B(17) 

R(17) 

R(36) 

Special Flow No. A ( Ib/hr) 

1 

A(18) 

B(l8) 

. R(18) 

R(37) . 

Special Flow No. 5 (ib/hr) . . ' 

f 

A(19) 

E(19) 

R(19) 

R(38) " 
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4.1 CHILLR CHECKOUT RUN 

The CHILLR subroutine was used to model a refrigerator and freezer of the 
.approximate size of those used on Skylab. For the refrigerator simulation, 
an externally chilled Freon-21 coolant was routed through cooling tubes 
embedded in the refrigerator walls. For the freezer simulation, a self- 
contained cooling unit was assumed with a constant coefficient of per- 
formance of 1.07. This value represents a feasible vapor-compression 
refrigeration unit once considered for use on Shuttle (Reference 22). 

For both cases a steady-state and transient solution were obtained, and 
the results presented in Sections 4.1.1 through 4.1.4. In addition, a 
CHILLR run simulating the Stirling cycle freezer designed for Shuttle 
(Reference 32) is described in Section 4.1.5. 

4.1'. 1 CHILLR Steady-State Refrigerator Case 

The 6-189A component input data for the refrigerator case are listed in 
Table 4-2. The final steady-state solution is tabulated in Table 4-3. The 
results for the run are plotted in Figures 4-1 and 4-2. As most Skylab 
refrigerator test data were not readily available, the model performance 
was compared with independent analysis. A steady-state condition was 
achieved within 10 iterations, as seen in the figures. Subsequent changes 
have been made within the subroutine to allow additional internal iterations 
thus achieving steady state conditions within only one or two system 
iterations. The final total energy input to the locker by conduction and 
door opening is seen in R(53) and Figure 4-2 to be 6.47 watts (22.06 Btu/hr) 
This is also equal to the heat transferred to the cooling fluid, as required 
For this condition, the coolant fluid temperature rise should be 

.-r ' ^ d _ 22.06 Btu/hr _ ^ rror 

^'coolant fluid • (25 lb/hr)(0.248 Btu/lb-°F) ^ 


For the inlet temperature of 37°F, the outlet coolant should therefore be 
40.56°F, which is the temperature shown in Figure 4-1. 
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TABLE 4-2 

G-189A INPUT DATA FOR REFRIGERATOR CHECKOUT RUN 


Id** 

97 

1 

REFr I GE 

KB AS 

97 

u 

7 1 

N S T R ' ■ 

9 7 

03001 D ' 

Vary 

97 

S 1 

70. 

Vary 

■97“ 

5 7 

70. 

Vary 

97 

5 5 

7 5 . 

V aR Y ■ 

97 

s; 

7 0. 

VaR y 

9 7 

58 

20 . 

vary 

97 

" 0 0 

70 . 

Vary 

97 

6 1 

,08 

V aRy ■■■ 

■97' 

• '6 J- 

7 D . 

Vary 

97 

6 '4 

. 7 

vary 

97 

6 5 

2 . 6 

V aR y 

9 7 

' Gy 

2.09 

Vary 

97 

7 C 

1 » 0 

Vary " 

97 ■ 

■ 7 I 

2 7 . “■“ 

VARY 

97 

72 

2 8,35 

vary 

97 

-73 

, 2 ■ ' “ 

Vary 

9 7 

7 7 

1 0 . 

Vary 

97 

' 75 

70. “ ■ 

Vary 

97 

78 

7 , 

Vary 

97 

78 

6.3 

Vary 

97 

8 2 

I 3 . 

Vary 

97 

9,!. 

70 .■ 

Vary 

97 

97 

.075 

Vary 

97 

9a ■ 

0 . 

Vary 

97 

9 V 

.27 

Vary- 

97 

103 

1 . -- 

Vary 

97 

1 CH 

37 . 

Vary 

97 

1 1 7 

1 77Qi 

Vary 

97 

J 1 b 

. 5 

Vary 

97 

1 1 8 

.5 

V aRy 

9 7 

1 19 

. 2 
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thermal conductivity Ratio, fRuzenvunfrozen 
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ambient aIR temp (F) 
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aMbieHt air humidity “ “ - 

ambient aIR specific HeaT 

COOl, INg coils thermal mass dry ■ - 
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DRY food specific HEaT 


D2.118571-I 






TABLE 4-3 

FINAL STEADY-STATE SOLUTION FOR REFRIGERATOR CHECKOUT RUN 


COMP' NO' "~9T~“CHaiR "SOaR NO' 71' PRr“sor 

failure flags -- COMP»_ 


SEC*“SOR- 

LOOP" 0 


-COMf-PASS-NO 10" 




.0 sec 


All)* 


Blip 


_ 25.000 

”'.'6oooo 

CPA* 

.00000 

.00000 

CpB* 


37.000 

".OOOOO' 

.zhboo 
'.00000 
.00000 
. DOOOQ 


.00000 

■'.00000 

wTHAS 

.00000 

.00000 

rtTHSu" 


ZH0»00_ 
t) 0 0 0 0~ 
IOZ.90 
'.ooooo' 
.00000 
'.ooooo ' 


.00000 . 

.00300 

khoa* 

.00000 

.00000 

RHOB»" 


•OOOOO 
“iODOOO 
66.800 
■'.OOOOO 
.00000 
■ . 0:0000 


•ooooo 

-.00000 

VI 5 CA* 

'.00000 

•ooooo 

vises*' 


•ooooo 

-.00000* 

.82200 

*.00000 

•ooooo 

-.00000 


•ooooo 
“•ooooo 

AKA* 

-.00000 ■ ■ 
•ooooo 

AKB* 


•ooooo 


. 61200-01 

-.00000 

• oootio 


I 

cn 


COM 

VR I 
VR 1 
VR 1 
vR 1 
VR I 
vR I 
vR 1 
VR I 
VR t 
vR 1 


NO 


1 ) " 
22 ) "■ 
B5 ) * 
ill* 
67 ) • 

7 3 1 * 

7 ? ) = 

BS i * 

9 1 ) = 

9 7 1 * 
VR I 1 03 ) •_ 
VRI 109 )* 
VRIIIS )* 

vR I 1 2 I ) * 


- 

— 

f'ai'lO're 

FLAGS -- 

C'd'hP 

■ 0 

LOOP* 

0 

elap 

25.000 

VRl 

2 ) “ 

H 0 .S 5 H 

VR 1 

3 )- 

2 H 0.00 

VRl 

H»* 

.00000 

VR 1 

231 * 

.00000 

VRl 

51 »- 

H 1 .779 

VRl 

521 * 

H 5 . 000 

VRl 

56 ) * 

“ 13 .H 60 

VR 1 

571 * 

70*000 

,VRl 

SSI* 

a snnr)Q*o 1 

V H 1 

62 ) * 

- 2.2578 

VRl 

631 * 

69.701 

VRl 

6 HI* 

.03000 . 

VR 1 

68 ) = 

2.6000 

' VKl 

69 ) ■ 

2.0900 

• VRl 

701 * 

. 2 0 C 0 0 

VR < 

7 H ) * 

1 0 • 000 

VR 1 

751 * 

H 1 . / 7 9 

VRl 

76 1 ■ 

H 0 . 55 H 

VR I 

OO ) = 

.00000 

VR 1 

8 1 ) * 

.00000 

VRl 

32 I ■ 

.00000 

VR I 

86 ) * 

.00000 

VR 1 

B 7 1 ■ 

.00000 

VRl 

58 1 * 

t OOOOO 

VRl 

92 ) « 

.00000 

VR 1 

93 ) ■ 

.00000 

VRl 

9 HI* 

. 75 OQO -01 

VR 1 

9 B ) » 

. ooooo 

VRl 

99 ) ■ 

. 2 HOOO 

VRl lOOl* 

1 .0300 

VR 1 1 OH ) ■ 

H 0 . 55 H 

VRl 1051 * 

.25908 

VRl 106 >*_ 

. OOOOO 

VR I 1 1 0 ) * 

.00000 

VRll 1 1 )- 

•OOOOO 

VR 1 11 2 1 ■ 

.50300 

VR 1 1 1 6 ) * 

.50000 

VRU 17 1 * 

•OOOOO 

VR 1 1 1 8 » ■ 

• 00000 

VR 1 1 22 ) ■ 

•OOGDO 

VR 1 1 23 ) » 

.00000 

VR 1 1 2 H ) * 



CpP- 

. 2 H 80 O 


ATMP* 

102.90 





CPS* 

•OOOOO 


«TMS* 

•00000 

• 

RHQS* 


COMP 

I ME' '1 N 


2H0.00 

.77258 

20.000 

• 70000 

I .0000 

"h.ooco " 
1 u.ooo 
.00000 
.00000 

• ooooo 

• ooooo 
".00000 

• ooooo 

22.Q32 

B6 • aoo 

. DODDO 


PASS NO 
"MlUtrSET' 
VRl 20 »■ 
VRi 53)*- 
VRI 5»)*_ 
VRl 65 ) •" 

VRl 71 )• 

VRl 77)* 
VRl 83)* 
VRl 89)*'" 
VRl 951*^ 
VRl IQI I* 
VRl 107>*_ 
VRl I 131* 
VRl I 1 9 I • 
VRl 125)* 

VlSCP* _ 

Vises* 


10 

»i-CP* 

.00000 

-22.063 

>6.0869 

22.063 

29.000 

■ HI .779* 
• OOOOO 
.00000 
•OOOOO 
.ooooo' 

_ .00000 
• OOOOQ ' 
• 2Q000 
.00000 
.62200 

■ .ooooo' 



• 0 see 


o 

ro 

00 

-oi 


pp* — 
VRl 2l»* 
VRl 5H)*“ 
VRl 60l* 
VRl 66»* 
VKl 721* 
VRr78»*'" 
VRl 8H 1 ■ 
VKl 901*' 
VRl 961* 
VR I I 02 I ■” 
VRl 10B»* 
VHll IH I •" 
VRIUQI* 
VKl 126)*' 
AKP* 
XKS* 


0 

tOOOOO 
70.000 " 

70.000 

HO.OOO 

26.350 

6.3000 
.00000 
.00000 
70.000 _ 

• OOOOO 
.00000 
IHHD.O ■ 

• OOOOO 
■ .OOOOO 

•6|200-0I 

. .00000 
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4.1.1 (Continued) 

The food inner and outer node temperatures are equal at steady ” 

they Should be, at 5.43«C (41.780F). This value results from equating heat 

in by conduction and door opening to heat out to the fluid. The beat input 

from door opening 's given by 

Odoor ' '=air '' S '^ambient ‘ '^food' 


= 0.075 


^ (0.51 cu ft) ^0.24 (70»F-41.78»F) 


=0.259 Btu 

This agrees with the value of R(105) in Table 4-3. Since 24 door openings 
per day were assumed in R(71). the steady-state average door opening heat 
added is 0.076 watts (0.259 Btu/hr). The effective thermal conductance 
between the food and ambient environment is seen in R(52), Table 4-3 o 
be 0 773 Btu/hr-°F. This value was determined from Skylab freezer an 
refrigerator d.esign data. The heat input by conduction to the food is 

then given by: 


^conduction in 


0.773 (70°F-41.7SeF) 


21.80 Btu/hr 


Thus the total heat into the refrigerator by conduction and door opening 
should be 0.259 + 21.80 = 22.06 Btu/hr = 6.47 watts. As shown previously, 
this is the same energy input computed by the program; thus, verifying the 
energy balance made in the subroutine is valid. 
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4.1,2 CHILLR Transient Refrigerator Case 

The G-189A input data for the transient refrigerator run were identical to 
^ the steady-state case shown in Table 4-2. A 4-hour transient solution was 
obtained. The final solution at the end of 4 hours is tabulated in 
Table 4-4. The transient history is plotted in Figures 4-3 and 4-4. Time 
did not permit running a longer solution; however, the temperatures in 
Figure 4-3 are approaching the required steady-state temperatures in 
Figure 4-1. (Note that a 10-hour refrigerator solution is shown in the 
Space Station model. Appendix G, for similar input conditions, showing the 
equilibrium food temperature is indeed approaching the steady-state value 
in Figure 4-1.) 

The heat balance may be verified by checking the energy conservation 
equation: 

' ""=p (ir) 


Checkinq this equation at a time of 4 hours, the heat input and output are 
seen in Figure 4-4 to be 22.5 and 17.9 Btu/hr, respectively. The net energy 
i.nput is therefore 4.6 Btu/hr. The heat storage must be computed as the 
sum of (m Cp AT/At) for the (1) inner food node, (2) outer food/walls/ 
shelves node, and (3) cooling coils. Taking the thermal mass from the 
solution in Table 4-4 and estimating the heauing rates from Figure 4-3 the 
total heat storage rate at 4 hours time is given by: 

Heat storage (inner food) =. (25,3 Btu/°F) (0.123°F/hr) 

= 3.11 Btu/hr 

Heat storage (outer food/walls) = (16.32 Btu/°F) (0.088°F/hr) 

= 1.43 Btu/hr 
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TABLE 4-4 


FINAL TRANSIENT SOLUTION FOR REFRIGERATOR CHECKOUT RUN 
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•00000 
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APPLIANCE CHECKOUT TRANSIENT 
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Figure 4-3. Temperatures for Refrigerator Transient 
Checkout Run 
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Figure 4-4. Heat Flow for Refrigerator Transient 
Checkout Run 
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4.1.2 (Continued) 

Heat storage (cooling coils) = (1.0 Btu/°F)(0.063°F/hr) 

= 0.06 Btu/hr 


Thus, the combined heat storage rate at 4 hours time is 4.60 Btu/hr. As 
shown previously, this is also equal to the net heat input to the refrig- 
erator, thus verifying the accuracy of the overall transient energy balance. 


m.L.R Steady-State Freezer Case (Vapor Compression) 

The G-189A component input data for the freezer case are listed in Table 4-5. 
The basic freezer locker design is identical to the refrigerator locker 
discussed previously, except that the cooling fluid is supplied by a self- 
contained 1 ef rigerati on unit with a coefficient of performance of 1.07. For 
this steady-st'^Le case, the refrigeration unit was assumed to be on con- 
tinuously at normal operating capacity. The final steady-state solution is 

tabulated in Table 4_6. The results for the run are plotted in Figures 4-5 
through 4-7. 

The net heat conducted into the food locker from the ambient environment 
is seen in Figure 4-6 to be 74.9 Btu/hr. At steady state, all this heat 
must be transferred directly to the cooling coils. The coefficient of 
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TABLE 4-5 

G-189A INPUT DATA FOR FREEZER CHECKOUT RUN 




ca 


'1 D» • 

-42 

- -f FREEZER" -“-WITH-SI 

K B A S 

42 

0 71 7 

NSTR 

42 

0101 1 1 \ 

KART 

42 

17 0 

vary 

4 2 

51 i-3' - 

Vary 

42 

Sh 70. 

Vary ■■ 

42 ■ 

-■SS MS. 

Vary 

42 

57 70. 

vary 

4 2 

5H 20. - ■“ . 

vary 

42 

60 70. 

VARY 

42 

6 1 . OB ' ■ 

vary 

42 

4J 70. 

vary 

42 

6M .7 . 

vary 

62 

6 fa -10. ' 

VARY 

62 

67 1.0 

Vary 

42 

6 2. a 6 

Vary 

42 

67 2.09 

Vary ' 

42 

7 0 I . 0 ' 

Vary 

62 

72 26.35 

Vary 

62 

7 3 a 2 

V^RY 

62 

7 m 10. 

V jRY 

62 

7 Ij -IQ. 

vaRt 

62 

74 M . 

Vary 

42' 

- 77-10. 

Vary 

42 

7 U 6 a 3 

Vary 

67 

77 -10. 

Vary 

42 

6 2 6 a 

Vary 

62 

8J -15. " 

vary 

62 

Bm -5. 

Vary 

42 
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vary 

42 

9/ .075 

Vary 

62 

’■•i 0. 

Vary 

62 

9v . 2M 

vary 

42 

1 on 1 . 07 

Vary 

42 

I 0 2 -28. 

vary 

42 

■ 1 0 J 1 . 

Vary 

42 

1 OM - to. 

■ Vary 

42 

1 15 .5 

Vary 

42 

1 1 a .5 

■ Vary 

62 
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Vary 

‘ 62 

l20 .08 

Vary 

42 
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vary 

62 

l2S .25 

Vary 

62 
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62 
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62 
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Vary 

42 

1 30 2. ' 
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TABLE 4-6 


COHP ’NO"' 


A (1 I » 


B ( I 1 ■ 


VR I 1 > " 
VRt 71- 
VR ( 131* 

VR ( S 1 ) • 
VR ! 57 1- 
VRt 63)* 
VR ( 6V 1 - 
VRt 7bl ■ 

V R I 6 1 ) ■ 
VRt 871* 
vR t 93 )■ 

VR I 99)» 
VR n Ob ) « 
vRtini- 
VRdl71- 

V R U 2 3 1 - 

V R t I 2 9 ) - 


FINAL STEADY-STATE SOLUTION FOR FREEZER CHECKOUT RUN 


-6Z CHILLR 
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4.1.3 (Continued) 

‘performance for the coolinq unit was input’ in R(IOO) 

, as 1.07 according to Reference 22. Thus, from equation (3.1.15), the 
total electrical input power required should be 

~ iToT " Btu/hr 


This value is correctly reflected in Figure 4-6. From equation (3.1.16), 
the total heat rejected by the cooling unit condenser to ambient should be 

qn == \ Btu/hr 

This value is also correctly reflected in Figure 4-6, thus verifying the 
accuracy of the thermal energy balance. 

4.1.4 CHILLR Transient Freezer Case 

The G-189A input data for the transient freezer run were identical to the 
steady-state case shown in Table 4-5. The basic freezer locker design is 
identical to the refrigerator locker discussed in Paraqraohs 4. 1 . 1 . and 4.1.2. 
except that the cooling fluid is supplied by a self-contained refrigeration 
unit with a coefficient of performance of 1.07. A 4-hour transient solution 
was obtained. Initially, the refrigeration unit was turned off. The locker 
gradually warmed until the food temperature reached -5°F (after about 2 
hours), when the refrigeration unit was turned on. The cooling unit remained 
on for the remainder of the run. The final results at the end of 4 hours 
are tabulated in Table 4-7. The transient history is plotted in Figures 4-8 
through 4-10. The solution during the'second 2 hours, when the cooling 
unit is turned on, is approaching the steady-state case discussed in Para- 
graph 4.1.3 and plotted in Figures 4-5 through 4 _y^ The temperatures and 
heat balance were checked against the analytical model in the same manner 
as was done for the previous cases and found to be valid. 
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Figure 4-8. Temperatures for Freezer Transient 
Checkout Run 
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4.1.5 CHILLR Simulation of Shuttle Stirling Cycle Freezer 

A preliminary design of a food and medical sample freezer using a Stirling 
' cycle refrigeration unit was developed for use on Shuttle (Reference 32). 

The CHILLR subroutine has been used to simulate that freezer, and excellent 
agreement has been obtained. The model input data for this simulation have 
been built into the subroutine as default data. These data are listed in 
Table 3-2. 

A detailed three-dimensional thermal model of the freezer is described in 
Reference 32. In the absence of hardware test data, the analytical model 
results were used to correlate the G-189A subroutine model. Both the G-189A 
subroutine and the independent three-dimensional model described above 
were run with identical operating conditions, and their results compared in 
Figures 4-11 through 4-16. The transient run, Figures 4-14 through 4-16, 
included five door openings near time 3 hours, and a fresh medical sample 
(with 2.15 lbs. water) initially at 80°F inserted at 6 hours. The results 
from the two models agree very closely— especially the overall heat 
balance, which is the most important G-189A consideration. For steady 
state operating conditions in a 70°F ambient environment with no door 
openings or medical sample insertion, the refrigeration unit duty cycle 
(fraction of total time on) for both models was 62 percent. 
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Figure 4-11. Correlation of G-189A Shuttle Freezer 
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Figure 4-12. Correlation of G-18SA Shuttle Freezer 
Steady State Thermal Solution 
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Figure 4-13. Correlation of G-1S9A Siiuttle Freezer 
Steady State Heat Balance 
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4.2 FTRAY CHECKOUT RUN 

The FTRAY subroutine v/as used to model a standard Skylab-type food heating 
V tray. Two cases were rim:' one steady state and one transient. For both 
cases the subroutine input default data, listed in Table 3-5, were used to 
simulate heating of a frozen meat product stored at -25°F to approximately 
150°F. 

4.2.1 FTRAY Steady-State Case 

The final component data for the steady-state food tray case.' are listed 
in Table 4-8. The results for the run are plotted in Figures 4-17 and 
4-18. The temperature data are compared with actual test data from Refer- 
ence 23, both for an average heater power of 7.9 watts. The test data were 
taken with a single thermocouple buried roughly near the center of the 
heated food. Since such a large temperature gradient exists within the 
food (43°F), it is not possible to compare temperatures at specific nodes 
with the test data point. However, the test data fall generally within 
the computed temperature range. In Figure 4-18, the heat dissipated to 
ambient surroundings is equal to the input heater power of 7.9 watts 
(26.95 Btu/hr) at steady state, as required by an energy balance on the 
component. 

4.2.2 FTRAY Transient Case 

The final component data for the. transient food tray run are listed in 
Table 4-9. The results for 3 hours of heating time are plotted in Figures 
4-19 and 4-20. This time was long enough to reach steady-state conditions, 
and the final temperatures are seen to reach the same values shown in 
Figure 4-17 for the steady-state case. Transient test data from Reference 
23 are plotted in Figure 4-19. These data were taken with a single ther- 
mocouple without a clearly defined location, but which was roughly near 
the center of the food. Since the temperature gradients within the food 
are so large, the test temperatures cannot be compared with any specific 
food node. However, they do fall roughly within the computed temperature 
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Figure 4-20. Heat Flow for Food Tray Transient 
Checkout Run 
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4.2.2 (Continued) 

range. The heat flow into and out of the food is shown in Figure 4-20. 
After 3 hours warming time, the total heat dissipated to ambient was 
7.9 watts (26.95 Btu/hr), which was also equal to the average heater 
input power. The average equilibrium heater power required in the tests 
of Reference 23 was also 79 watts, thus indicating good correlation 
between the G-189A model and Sky lab food tray performance. 


4.3 ROSMOS CHECKOUT RUN 


The final design of the reverse osmosis unit to be selected for future space 
missions has not yet been determined, and consequently experimental data 
could not be used for comparison. However, a conceptual design was assumed 
for analytical simulation with currently feasible model data used as input. 
Two cases were run: one steady state and one transient. The component 

input data for the two cases v/ere identical and are listed in Table 4-10. 
These data correspond to the subroutine input default data listed in 
Table 3-9 except that a thermal conduction path between the unit and ambient 
surroundings has been added in R-array locations 54 through 64. In its 
planned mode of operation, the reverse osmosis unit will operate near 
ambient temperature, so no thermal path was included in the default data. 

For this case, the fluid outlet temperature is equal to its inlet tempera- 
ture. However, some unit designs require elevated temperatures, so this 
case was run to check out the subroutine thermal balance. 


4.3.1 ROSMOS Steady-State Case 


The final component data for the steady-state reverse osmosis unit simulation 
are listed in Table 4-11. The results for the run are plotted in Figures 4-21 
through 4-23. The outlet brine flow rate should be computed as described in 
Paragraph 3.3, equations (3.3.7), (3.3.9), and (3.3.10). For the input data 
assumed for the component, these calculations should yield: 


aP 


osmosis 


= R(75) 


\ini 

= 200 

’ 0.008* 

LR(71)J 

C_W\J 

0.01 


= 160 psia 


111 


productd^si 


R(70) 


gn 


= (9)(0.96) = 8.64 
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TABLE 4-10 

Grl89A INPUT DATA FOR REVERSE OSMOSIS UNIT CHECKOUT RUN 
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VARY 

7 

SR 

75, 

VARY 

7 

SS 

3 . 

VARY 

7 

67, 

75, 

VARY 

7 

58* 
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thermal capacitance of unit - DRY 


D2-II857I-I 




TABLE 4-11 

FINAL STEADY-STATE SOLUTION FOR REVERSE OSMOSIS UNIT CHECKOUT RUN 
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Figure 4-21. Flow Rates for Reverse Osmosis Unit 
Steady-State Checkout Run 
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4.3.1 (Continued) 

[ 4(4) - R(74) - 

"productj^jjij, ■ '"ptoductj^j^.g^ L '^(”) - R(67) - R(75) J 
= 8.64 = 9.770 Ib/hr 

The proiduct flow rate thus calculated agrees exactly with the flow plotted 
in Figure 4-21 and the actual value [R(20)] shown in Table 4-11. The brine 
flow is seen in Figure 4-21 and Table 4-11 to be the difference between the 
feed (inlet) and product flows as required. 

The thermal data in Figure 4-23 show that, in steady-state conditions, the 
heat into the reverse osmosis unit (from the warm fluid flowing through it) 
is equal to the heat dissipated through the unit structure to ambient 
surroundings. The energy input from the fluid is 

q^, . _ = m c„ AT = A(l) * CPA [A(2) - R(2)l 

^fluid in p 

= (10.2 lb/hr)(l Btu/lb °F)(165°F - 160.67°F) = 44.2 Btu/hr 


The energy out to ambient surroupdings by conduction, convection, and 
radiation is given by 

W “ Wtive ‘T = 9(52) (9(61) - R(54)l 

= (0.516 Btu/hr °F)(160.67°F - 75°F) = 44.2 Btu/hr 

Thus, the predicted reverse osmosis unit temperature agrees with the steady 
state heat balance as required. 
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4.3.2 RQSMOS Transient Case 

The final component data for the transient case are listed in Table 4-12. 

The results for 60 minutes of operation are plotted in Figures 4-24 through 
4-26. The conditions assumed for the transient case were identical to those 
used in the steady-state case discussed in Paragraph 4.3.1. The 60 minutes 
operating time was sufficiently long to achieve near equilibrium conditions, 
and the final results may be compared directly with the steady-state solution 
in Table 4.11 and Figures 4-21 through 4-23. It is apparent in the figures 
and tabular data that the final transient and steady-state solutions are 
identical, as they should be. 

4.4, SHOWER CHECKOUT RUM 

The SHOWER subroutine was used to model the zero gravity whole body shower 
described in References 19 and 24. Only the conditions in the new shower 
stall component are simulated in the checkout run since the peripheral 
equipment (pump, fan, heater, etc.) can be modeled using standard G-189A 
subroutines already operational. Note that the complete shower subsystem 
is included in the Space Station model described in Section 6. A steady- 
state and transient case were run, as described in the following paragraphs. 
For both cases the subroutine input default data listed in Table 3-12 

were used. 

4 ^ 4_1 Steady-State Shower Case 

The final solution for the steady-state shower run is given in Table 4-13. 
The results of the run are plotted in Figures 4-27 through 4-30. The 
inlet air and water flow rates and temperatures were set to actual test 
conditions. For the first 10 iterations, the shower water was turned on. 
Thus, the relative humidity in the stall was assumed in the model to be 
100 percent as Y'eflected in Figure 4-29. For the second 10 iterations, 
the inlet water flow was turned off, and evaporative drying within the 
stall was computed. In Figure 4-28 is shown the total heat balance for 
tho shower stall. At steadv-state. both for water-on and water-off con- 
ditions, it is apparent that the heat input is equal to the total heat 
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TABLE 4-12 


FINAL TRANSIENT SOLUTION FOR REVERSE OSMOSIS UNIT CHECKOUT RUN 
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TABLE 4-13(a) 

FINAL STEADY-STATE SOLUTION FOR SHOWER CHECKOUT RUN WITH WATER TURNED ON 
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TABLE 4-13(b) • 

FINAL STEADY-STATE SOLUTION FOR SHOWER CHECKOUT RUN WITH WATER TURNED OFF 
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Figure 4-27. Temperatures for Shower Steady-State 
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Figure 4-29. Evaporation Parameters for Shower 
Steady-State Checkout Run 
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4.4.1 (Continued) 

out as required. The oxygen and CO 2 pressures within the stall, in 
' Figure 4-30, agree -with the computed values from equations (3.4.29) 
through (3.4,32). The shower temperatures, Figure 4-27, were checked 
by hand calculations and found to agree with the thermal balance given 
in Paragraph 3.4.3 for the walls and occupant nodes. During the air 
drying period, with water turned off, the results in Figure 4-29 show 
an evaporation rate slightly over 0.5 Ib/hr. This value agreed favorably 
with limited test results in Reference 19. However, additional test 
data would be desirable to further refine and verify this result. The 
subroutine may be readily modified to fit other data using the correlation 
factor Cf^ discussed in Section 3. 4. 3. 2. 

4.4.2 Transient Shower Case 

The G-189A input data for the transient shower run were identical to the 
steady-state case. For the first 1% hours, the water was turned on and 
the remainder of the time it was off. Final solutions for the v/ater-on 
and water-off conditions are given in Table 4-14. The results for the run 
are plotted in Figures 4-31 through 4-34. It is evident in the plots that 
both steady-state and transient runs reach the same final solution as 
required. 


4.5 WASDRY CHECKOUT RUN 

The WASDRY subroutine was used to model a clothes and dishes washer/dryer 
drum. The peripheral equipment (pump, fan, accumulator, etc.) were not 
included in the checkout runs described here since they can be modeled 
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TABLE 4-14(b) 

FINAL TRANSIENT SOLUTION FOR SHOWER CHECKOUT RUN WITH WATER TURNED OFF 
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Figure 4-31. Temperatures for Shower Transient 
Checkout Run 
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4.5 (Continued) 

using standard G-189A subroutines already operational. The complete 
V washer/dryer subsystem is included in the Space Station model described 
in Section 6. 

To minimize the high power consumption of an electrical clothes dryer, it 
is recommended in Reference 21 that a drying time of 4 hours be allowed. 

This will permit much smaller air flow rates, and possibly lower inlet 
temperatures, than are used in commercial dryers. Also, a much faster 
spin-dry speed could be used to remove more of the water before drying. 

Typical household clothes washers achieve a 60 to 100 percent water reten- 
tion (lb water/lb dry clothes) in the clothes (Reference 25). However, a 
high-speed spin is reported in Reference 26 to achieve an 11 percent water 
retention. Therefore, commercial dryer data could not be used to correlate 
the dryer model. The model input data were made to correspond to the selected 
Space Station washer/dryer in Reference 1, and they also agree with the data 
from References 21 and 27. Since these represent conceptual designs only, 
and not current hardware, it was necessary to estimate reasonable values 
for some input data (such as thermal conductors within the washer frame) 
which will depend on the final design details. 

Steady-state and transient checkout runs were made. Results were obtained 
for four cycle phases: (1) wash water fillj (2) wash circulate; (3) spin- 

dry, wash water out; and (7) dry. The rinse cycles 4, 5, and 6 are compu- 
tationally equivalent to the wash cycles and hence were not considered 
here. They are included, however, in the all-up Space Station model in 
Section 6. 

4.5.1 Clothes Washer/Dryer Steady-State Run 

Both the steady-state and transient clothes washer/dryer runs used the sub- 
routine input default data listed in Table 3-16. The final steady-state solution 
is tabulated in Table 4-15. The steady-state results are plotted in Figures 
4-35 through 4-39. These figures show the solution with the unit turned 
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TABLE 4-15 

FINAL STEADY-STATE SOLUTION FOR CLOTHES WASHER/DRYER CHECKOUT RUN 
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Figure 4-35. Operational Cycle Phase for Clothes Washer/Dryer 
Steady-State Checkout Run 
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Figure 4-37. Teinperatures for Clothes Washer/Dryer 
Steady-State Checkout Run 
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Figure 4-38. Heat Flow for Clothes Washer/Dryer 
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4.5.1 (Continued) 

off and for cycle phase 1 (wash water fill), 2 (wash-circulate) , 3 (spin- 
^ dry), and 7 (dry). Since the solution was steady state, rather than transient, 
the amount of water in the drum was held constant, even in cases' (such as 
wash water fill) where there was a net flow of water into or out of the drum. 

During the initial phase with the unit turned off, all fluid flow rates and 
heat flow in and out of the drum are seen in Figures 4-36 and 4-38 to be 
zero; and the washer temperature was equal to the assumed ambient tempera- 
ture of 70°F as required. During phase 1, wash water fill, the heat input 
from the assumed 50 watt drum motor was 170.6 Btu/hr. The net heat conducted 
through the structure to ambient is shown in Figure 4-38 to be 450 Btu/hr. 

Thus, the total heat input from the inlet fluid must have been 450-170.6 = 
279*.4 Btu/hr. For the inlet water flow rate of 300 Ib/hr, its energy con- 
tribution may also be expressed as: 

^input from water ” 

279.4 Btu/hr = (300 lb/hr)(l Btu/lb °F)(150°F - (4.5.1) 

Solving the above equation results in a required final water (and tub) 
temperature of 149°F. This is the actual predicted temperature shown in 
Figure 4-37, thus verifying the accuracy of the energy balance. 

For the steady-state run, the next cycle phase 2 (wash-circulate) differs 
from phase 1 only in the inlet water flow rate which was 4000 Ib/hr as 
shown in Figure 4-36. This flow rate was recommended in Reference 21 to 
provide adequate agitation for washing. Substituting this flow rate into 
equation (4.5.1) results in a final water and tub temperature of approxi- 
mately 150°F. This value is reflected in Figure 4-37. 
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4.5.1 (Continued) 

During phase 3 (spin-dry), the only heat input to the washer drum is from 
, the 150 watt motor or 511.8 Btu/hr. Since all this heat is applied to the 
frame node in the thermal model, the tub and frame should be at the same 
temperature recording to the following relation: 


^MOTOR 


- r (T 

frame-to-ambient ' 


tub ' "'^ambient^ 


511.8 Btu/hr = 8.216 Btu/hr-°F - 70®F) 


Solving this equation results in a washer tub and water temperature of 
132. 2°F, which is also the value reflected in Figure 4-37. 

The final steady-state solution was for the drying phase 7. There was 
assumed to be 0.4 lb' of water [R(67)l in the clothes load of 4 lb [R(66)]. 
This gave a water loading of 10 percent as shown in Figure 4-39. For this 
condition, and with 12 Ib/hr inlet air flow rate, the predicted evaporation 
rate was 0.185 Ib/hr in Figure 4-39. Since typical household clothes dryers 
use much higher air flow rates, these data cannot be correlated with con- 
ventional dryer tests. If test data become available in the future for a 
spacecraft dryer, correlation of the model with the data may be made easily 
using the evaporation rate correlation multiplier, R(89), described in 
Paragraph 3. 5, 3. 2. For the 0.185 Ib/hr evaporation rate predicted here, 
the latent cooling effect would be 193 Btu/hr. This is nearly balanced by 
the drum motor input power of 50 v^atts or 171 Btu/hr. In addition, the heat 
conducted through the dryer structure and dissipated to ambient is given by 


*^structure loss 


= G at = (8.076 Btu/hr °F)(39.G°F " 70°F) 


= 158.3 Btu/hr 
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4.5.1 (Continued) 

Thus, the net heat input required from the air flow should be (193-171+158)= 
180 Btu/hr. This heat should be governed by the equation 

« 

*^out of air ■ ^"^in ’ ^out^ 

180 Btu/hr = (12 Ib/hr) (0.243 Btu/lb °F)(150°F - T^^^^^) 

and the final air outlet temperature should thus be 88.3°F. This is indeed 
the predicted outlet ai r temperature in Figure 4-37, thus verifying the 
accuracy of the energy balance in the model. 

4.5.2 Clothes Washer/Dryer Transient Run 

The G-189A input data for the transient clothes washer/dryer run were 
identical to the steady-state case. The final 
transient solution is tabulated in Table 4-16. The transient plots are 
shown in Figures 4-40 through 4-44. The washer was initially off, then 
cycled through phase 1 (wash water fill); 2 (wash-circulate) ; 3 (spin-dry, 
wash water out); and 7 (dry). The rinse phases were omitted here since 
they are equivalent to the wash phases. (Note that they are included in 
the all-up Space Station run in Section 6.) 

During the initial phase with the unit turned off, all flow rates and heat 
flow are zero; and the washer temperatures are equal to the assumed ambient 
temperature of 70° F. The solution during phase 1 (wash water fill) can be 
seen to be generally approaching the steady-state solution discussed pre- 
viously. The solution during phase 2 (wash water circulate) does reach the 
steady-state solution described previously. The drying phase 7 can be com- 
pared at time 2.5 hours with the steady-state case discussed previously 
since the water loading in the clothes is the same for both. The tv^o 
solutions are very nearly equal, with the transient temperatures being only 
slightly higher due to the stored energy in the dryer drum from the previous 
wash cycle. 
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TABLE 4-16 

FINAL TRANSIENT SOLUTION FOR CLOTHES WASHER/DRYER CHECKOUT RUN 
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4.5.3 Dishes Washer/Dryer 

Both the steady-state and transient dishwasher/dryer runs used the sub- 
routine input default data listed in Table 3-16. The final steady-state 
. solution is tabulated in Table 4-17. The steady-state results are plotted 
in Figures 4-45 through 4-49. The transient solution is shown in Table 
4-18 and Figures 4-50 through 4-54. These results were checked in the same 
manner as was done for the clothes washer/dryer case and found to be valid. 

4.6 WASTEC CHECKOUT RUN 

The WASTEC subroutine was checked out by modeling the dryjohn planned for 
use on Shuttle. Only the collector and urinal were modeled for the checkout 
runs since the peripheral equipment (pump, valves, lines, etc.) can be 
modeled using standard G-189A subroutines already operational. The complete 
dryjohn subsystem is included in the Shuttle and Space Station models in 
Sections 5 and 6. 

A steady-state and transient run were made as described in the following 
paragraphs. The subroutine input default data, listed in Table 3-19, were 
used for both cases, except for two minor exceptions. First, the initial 
collector temperature in R(51) was set at 70°F (instead of the 6 Qof 
default value). Secondly, the vacuum evaporation correlation multiplier, 
R(75), was input at 2.5 (instead of the 3.0 default value). The reason 
for this slight difference is due to subsequent data analysis which showed 
the default values to be slightly more realistic. The default values 
v/ere used in the Space Station models in Section 6. 

The three possible usage phases described in Paragraph 3.6 were simulated: 
(1) urine collection; (2) fecal collection; and (3) combined urine/fecal 
collection, as well as phase 0 (unit off). During phases 0 and 1, the 
collector contents were assumed to be under vacuum. 

4.6.1 Dryjohn Steady-State Run 

The final steady-state solution is tabulated in Table 4-19. The steady- 
state results for the four usage phases are plotted in Figures 4-55 through 
4-58. In each case, steady-state conditions were achieved in about five 
iterations or less. 
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TABLE 4-17 

FINAL STEADY-STATE SOLUTION FOR DISHWASHER/DRYER CHECKOUT RUN 
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TABLE 4-18 

FINAL TRANSIENT SOLUTION FOR DISHWASHER/DRYER CHECKOUT RUN 



1 1 ” 

•COOOO 

VR { 

2 J" 

109^99 

VR< 

3)- 

•00000 

VR( 

9l- 

• OOOOO 

VRl 

20 )• 

•OOOOO 

vr« 

*!*• 

IOV*RR 

VH ( 

22)- 

•00000 

VR< 

23)- 

.00.000 

VR ( 

29 I-'. 

.00000 

VRl 

251- 

• OOOOO 

VRI 

26 ) • 

•OOOOO 

VRI 

2T|» 

,00000 

V(? ( 

23 1 » 

•00000 

VR t 

29 > • 

.00300 

VR 1 

30)- 

•00000 

VRl 

3D- 

• OOOOO 

vrc 

32)» 

•OOOOO 

VrI 

331- 

•ooooo 

VR t 

5 I ) " 

95.57*) 

VR < 

52l- 

8.0953 

VR 1 

53 ) - 

209.09 

VR( 

59)- 

70.000 

VRl 

55) ■ 

3.0000 

VrI 

56 I • 

77.985 

VR ( 

S7 ) « 

70.000 

VR ( 

58 ) - 

I *0000 

VR ( 

59 ) - 

28,291 . 

VR ( 

6o ) - 

70.000 

VRi 

61)- 

9.0000 

VrI 

62)- 

103,31 

VR ( 

4 3 ) • 

95.828 

VR ( 

69 ) - 

•ccooo 

VR 1 

65)- 

•ooooo 

VR 1 

66 ) - 

16.200 

VRl 

67)- 

1 5 . 000 

VrI 

68 1 ■ 

50.000 • 

VR ( 

49 J • 

•OOGOO 

VR ( 

7o>* 

1 5c .00 

VR ( 

71 )- 

•00000 

VR( 

72)- 

109.99 

VRl 

73)- 

• OOOOO 

VRl 

79 1- 

• OOOOO 

VR ( 

75 ) » 

• 2SC00 

VR ( 

74 ) ■ 

12.000 

VR < 

771- 

•25000 

VR 1 

78)- 

. lOoCO-03 

VRl 

79 ) - 

290.00 

VRl 

BO ) • 

I0587t 

VR ( 

8 1)- 

8 . 0000 

VR ( 

32 ) - 

16. 000 

VR ( 

83 ) - 

120.00 

VR t 

89 ) - 

• 95oCO 

VRl 

65 ) - 

• OOOOO 

vrI 

B6I- 

.80000 

VR { 

37 ) • 

• 9 5000 

VR 1 

88 ) - 

1 .0000 

VR 1 

(19)- 

•90000 

VR« 

9q)- 

• ooooo 

VRl 

91)- 

♦ 1 8799 

VRl 

921- 

.00000 

VR t 

93 ) ■ 

• 00000 

VR < 

?9 ) » 

. 19790-01 

VR ( 

95 )- 

•OOOOO 

vrc 

96)- 

• 13560-02 

VRl 

97 ) - 

30.000 

VR 1 

98 ) ■ 

30.000 

VR ( 

99 1 - 

29.250 

VR 1 

too’* 

•GCDOO 

VR(lOl)- 

•OOOOO 

VR( 102)- 

15.000 

VRl 103)“ 

is.ooa 

VrI 109 1 • 

.25000-01 

VR{ 105)« 

1 0 . 000 

VR t 

)04l- 

9.CC00 

VR<107)- 

10.000 

VR( 103)- 

- 1 .OOOO 

VHl 109)- 

• 1 0000-01 

VRl 1 101- 

6,0000 

V R ( i 1 1 1 • 

•00000 

VR 1 

112)“ 

•ocooo 

VR(1131- 

•OOOOO 

VR(I19)- 

•OOOOO 

VR 1 1 1 5 )■ 

•ooooo 

VrI||6)- 

.00000 

VR 1 1 1 7 I • 

•00000 

VR t 

1 1 8 ) ■ 

•00000 

VR ( 1 1 9 1 - 

•ooooo 

VR 1 1 20) - 

• OOOOO 

VRl 1211- 

•OOOOO 

Vr 1 1 22 ) ■ 

•OOOOO 


D2-I18571-1 


OISHVIASS'IER CYCLE PHASE 


s.Mc.tt , APPLIANCE CHECKOUT Tl 

7 . I 1 1 — lOOOOQOOOOiO OO - O - 0 1 » 


.0 1,2 1,6 

HISSION TINE 


Figure 4-50. Operational Cycle Phase fo 
Transient Checkout Run 



DISHES WASHER/DRYER FLOWRATES CLB/HR) 



MISSION TIME - HOURS 


0 INLET WATER I OUTLET WATER 2 INLET AIR 

3 OUTLET AIR 

Figure 4-51. Flow Rates for Dishwasher/Dryer 
Transient Checkout Run 


4-88 


Ikl 


DISHES VJASHER/DRYER TEMP »»* DEG F 


D2-1I857N 



0 inlet 1 OUTLET CTU3) 


Figure 4-52. Temperatures for Dishwasher/Dryer 
Transient Checkout Run 
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Figure 4-54. Evaporation Parameters for Dishwasher/Dryer 
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FINAL STEADY-STATE SOLUTION FOR DRYJOHN CHECKOUT RUN 
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4.6.1 (Continued) 

PHASE 1 - Urine Conection • ' 

’“The first steady-state case run was phase 1 (urine collection). The 
subroutine assumes -a urinal flow addition equal to the micturition flow 
rate which was input in R(69) as 108 Ib/hr. The default value of 97°F 
was used for the urine inlet temperature. The inlet air flow rate was 
85.2 Ib/hr total at 70°F, of which 0.665 Ib/hr was water vapor. The outlet 
water vapor flow rate was 2.007 Ib/hr. Therefore, 1.342 Ib/hr of water were 
evaporated in the ur1r,c/air mixing process. The latent cooling from this 
source would be approximately 1.342 x 1060 = 1422.5 Btu/hr. The heat con- 
ducted from the fluid to the dryjohn (through conductor in Figure 3-16) 
is given by 

%rinal fluid to dryjohn ^ ^^u ^^urinal out ‘ "'^dryjohn^ 

= R(78) [R(21) - R(51)l 
= 0.15 Btu/hr °F (81.7°F - 58.4°F) 

=3.5 Btu/hr ’ 


Therefore, the net heat out of the urinal fluid mixture was 1422.5 + 3.17 = 
1426 Btu/hr. This should correspond to the predicted sensible heat loss 
and corresponding temperature drop of the air and urine. This is given as 
follows: 


*^air sensible loss 


= m Cp aT = (85.2 1b/hr)(0.243 Btu/lb°F) (70^F-81 .7°F) 


= -242.2 Btu/hr 


^urine sensible loss 


= m c. 


at = (108 lb/hr)(l Btu/lb°F)(97 - 81.7) 


= 1652.4 Btu/hr 
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4.6.1 (Continued) 

Thus, the net sensible heat loss of the fluid is about 1652.4 - 242.2 = 1410 
Btu/hr. This is approximately equal to the required heat loss shown ^revi- 
' ously to be 1426 Btu/hr. The reason for the difference is that a more 
accurate and slightly different method is used in the subroutine than the 
simplified equations given previously to check the results. However, the 
agreement is within the limits of normal uncertainties in the model input 
data, and the urinal heat balance and outlet temperature are valid. 

The collector contents were assumed to be under vacuum during this phase. 
Their temperature is given in Figure 4-56, and it is only slightly warmer 
(0.2°F) than the solution shown for phase 0. The results for that phase 
are discussed later. The 0.2°F temperature increase during phase 1 is 
caused by the heat transferred in from the urinal which is in use. 

PHASE 2 - Fecal Collection 

During this phase, 68.2 Btu/hr are input to the collector from the 20 watt 
slinger motor as reflected in Figure 4-57. in steady state, all this heat 
should be dissipated to the air and structure. The heat input to the 63.9 
Ib/hr collector air flow is given by 

q . = m c„ AT = (63.9 lb/hr)(0.243 Btu/lb°F) (71 .63°F-70°F) 

SI y' p 

= 25.0 Btu/hr 

The heat loss through the dryjohn structure to ambient is given by 

^structure = ^ = (2.33 Btu/hr °F)(88.47°F - 70°F) 

= 43.0 Btu/hr 
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4.6.1 (Continued) 

Thus, the total heat loss to the air and structure is 25 + 43 = 68 Btu/hr 
'Which is equal, as shown previously, to the total heat input as required. 
No evaporative cooling is involved during this phase since the collector- 
is not under vacuum. 

PHASE 3 - Combined Urine/Fecal Collection 

The urinal temperatures for this case are shown in Figure 4-56 to be 
approximately equal to those described previously for phase 1 (urine 
collection) as required. The collector temperatures are shown in the same 
figure to La equal approximately to those described previously for phase 2 
(fecal col iction) as expected. The heat flow in Figure 4-57 is also 
approximately equal to that for phase 2 as expected. 

PHASE 0 - Unit Off 

The results for phase 0 (unit off) with collector being vacuum dried were 
included at the end of the run as shown in Figures 4-55 through 4-58. The 
collector contents are shov/n in Figure 4-56 to be at 58.1°F. As seen in 
Figure 4-58, the collector outlet is saturated-; that is, only water vapor 
is being outgassed from the fecal material. Water vapor pressure at 
58.1°F is approximately 0.24 psia. The evaporation rate, from equation 
(3.6.6), should therefore be 


m 


evap 



P 

0 ) 



= (2.5)(1.0 sq ft)(0.24 psia) 
^460 + 58.1°F 

= 0.0264 Ib/hr 


KEPEOBUCBlLl'rY OP THE 
IS POOB 
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4.6.1 (Continued) 

As seen in Figure 4-58 and Table 4-19, this rate is evaluated correctly. 

^ From actual dryjohn data. Reference 28, it is reported that 4 hours are 
required for new fecal material to lose one-half of their initial water 
content. Assuming 0.3 lb fecal mass per defecation with 75 percent water 
content, this would require an average evaporation rate of 0.028 Ib/hr 
to lose one-half its moisture in 4 hours. Of course, this rate will vary 
widely with time and with different fecal properties, and the rate pre- 
dicted in the subroutine is within a typical range of operating conditions. 
For the evaporation rate of 0.0264 Ib/hr, the cooling effect on the dryjohn 
should be (0.0264 Ib/hr) x (1042 Btu/lb) = 27.5 Btu/hr. This value is 
correctly reflected in Figure 4-58. This should also be equal to the heat 
conducted into the collector structure, as given by 


0 . = 6 .. X. aT = (2.32 Btu/hr°F)(70°F-58.13°F) 

^structure effective “ ' 

= 27.5 Btu/hr 

Thus, the energy into and out of the collector is equal, as required, and 
the computed collector contents temperature is valid. 

4.6.2 Dryjohn Transient Run 

The input model data for this case were identical to the data used for the 
steady-state case. The final solution is tabulated in 

Table 4-20. Four usage phases were run: (1) urine collection, (2) fecal 

collection, (3) combined urine/fecal collection, and (0) unit off with 
vacuum dry. The transient results are plotted in Figures 4-59 through 4-62. 
These results may be compared with the steady-state solution described pre- 
viously. The primary difference between the two solutions is the urinal 
outlet temperature during phases 1 and 3. For steady state, a continual 
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Figure 4-60. Temperatures for Dryjohn Transient Checkout Run 
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Figure 4-61. Heat Flow for Dryjohn Transient Checkout Run 


4-104 


1 

D2-118571-I 


■ill'* i»'r • 

24 . 


(/) 

UJ 

K 

LU 

r 

< 

(r 

< 

Q. 

(J) 

cn 

LU 

o 

o 

or 

a. 

o 

z 

*— I 
>- 
or 

Q 

<r 

o 

»- 

o 

LU 

_J 

o 

o 


20 . 


16. 


12 . 


8 . 


4 . 


APPLIANCE CHECKOUT TRANSIENT 


1^ 


< 1 1 ' 


.0 


.5 


3.5 


l.O . 1.5 2.0 2.5 3.0 

MISSION TIME - HOURS 

WATER LOAD (LB/LB) EVAP RATE CLB/HR) COMMODE REL. HUMDTY 


4.0 


Figure 4-62. Collector Evaporation Parameters for 
Dryjohn Transient Checkout Run 
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4.6.2 (Continued). 

urine outlet flow equal to the input micturition flow rate was assumed. 

For the transient run, however, this flow was limited to the average urine 
generation of 0.6 Ib/use which was input, and then it was stopped. Con- 
sequently, the urinal outlet temperature did not rise appreciably above 
ambient. 
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5.0 SHUTTLE ORBITER APPLIANCES SIMULATION 


To demonstrate that the appliance subroutines are operational in an all-up 
^ G-189A system simulation, they were run in a Space Shuttle Orbiter and 
Modular Space Station G-189A ECLSS model. The appliances included in the 
models were taken from trade studies of many candidate concepts (Reference 
1)' and represent the optimum set of appliances for each vehicle. The 
results of these trade studies for a four-man, 20. 5-day Shuttle Orbiter are 
shown in Table 2-1. Due to the short mission duration, the most simple 
concepts were selected in most cases. The only Shuttle appliances selected 
which required the new G-189A appliance subroutines were a space radiator 
refrigerator, food heating/serving trays, and a dryjohn. Since a refriger- 
ator is currently not included in the Shuttle design and since the Shuttle 
cabin coolant loop would be only marginally effective for the space radiator 
concept, the refrigerator was not included in the Shuttle model. Both 
refrigerators and freezers are included in the Space Station model in 

Section 6. A food and medical sample freezer has been developed for use on 
Shuttle, as described in Reference 32. This freezer is not included in the 
Shuttle model. However, the CHILLR subroutine has been adapted to simulate 
the freezer, as described in Sections 3.1.1 and 4.1.5. 

5.1 UNMODIFIED SHUTTLE CASE 

An available steady-state G-189A model of the Shuttle Orbiter ECLSS 
(described in Reference 3)', together with several model modifications supplied 
by Me Donnell Douglas Corporation (MDAC), was used for the basic Shuttle 
system model. This model obtains a single steady-state solution for up to 
23 different mission phases. For the purpose of verifying the appliance 
subroutines, this model was run for mission phases 12, 13, 14, and 15 which 
represent four different days of normal orbital operations. Two cases were 
run. First, an unmodified case was run with all program inputs exactly as 
supplied by MDAC. No appliances or modifications were included. The input 
model data and GPOLY subroutines for this run are included in Appendix B. 

The output tables and plots from this run are listed in Appendix C. 
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5.2 SHUTTLE CASE WITH APPLIANCES ADDED 

A dryjohn and food heating trays were then added to the basic Shuttle model 
described previously, and the same case was rerun. The changes made in 
the basic model and 6P0LY subroutines are listed in Appendix D. A flow 
schematic of the G-189A components added to the basic Shuttle case is 
shown in Figure 5-1. The dryjohn air inlet and outlet were connected 
-directly to the secondary side of component 2, which was the Orbiter cabin 
in the basic model. The dryjohn collector and urinal outlets (component 121) 
were directed to a water separator (component 125). This water separator 
also had a fluid inlet from the airlock to simulate the actual Shuttle 
configuration. However, since the basic Shuttle model used did not have 
the airlock included, a dummy input (component 130) was used for this con- 
nection with a flow rate of zero to simulate normal crew work activities. 
Another dummy input from the air revitalization system waste water (com- 
ponent 130) to the waste water storage tank (component 129) was also 
included to make the model complete, although it was also given a flow rate 
of zero. The heat conducted through the dryjohn structure to the cabin 
(stored in the dryjohn R-array location 53) was added in GPOLYl to the total 
cabin input heat load in cabin component R-array location 66. The waste water 
storage tank (component 129) was included in the model. However, the heat 
dissipated through the tank structure was not added to the cabin heat load 
since the tank model, not being part of the dryjohn system, was not suffi- 
ciently correlated to accurately simulate the actual Shuttle tanks. 

The food heating trays were not connected directly in any fluid flow loop; 
thus, no components are shown attached in Figure 5-1. The heat dissipated 
from the food trays to the cabin (stored in the food tray R-array location 
53) was added in GPOLYl to the total cabin heat load in cabin component 
R-array location 66. 

To exercise various appliance subroutine options, different conditions were 
assumed during each of the four mission phases. These are shown in Table 5-1. 
The output tables and plots from this run are listed in Appendix E. The 
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Figure 5-1. Flow Schematic of G-189A 'Appl i ance Components 
Added to Basic Shuttle Orbiter Model 
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5.2 (Continued) 

identical parameters were plotted here which were plotted previously in 
^Appendix C for the unmodified case. All differences between the two solu- 
tions are attributable to the addition of the dryjohn and food trays. In 
addition, plots were made of the dryjohn and food trays performance during 
the four mission phases in Figures 5-2 to 5-10. In all cases, they are 
seen to reach approximately the same solutions which were presented for the 
individual subroutines in Paragraphs 4.2 and 4.6, thus demonstrating that 
the new subroutines are working properly in an all-up 6-189A system model. 


TABLE 5-1 

APPLIANCE COMPONENT OPTIONS ASSUMED 
IN SHUTTLE ORBITER SIMULATION 




D R Y J ( 

) H N 


FOOD TRAYS 

MISSION 

USAGE 


CABIN AIR FL0W,CFM 

TOTAL NUMBER 

PHASE 

PHASE 

DESCRIPTION 

URINAL 

COLLECTOR 

BEING HEATED 

12 

0 

Not in use 
Vacuum dry 

0 

0 

0 

13 

1 

Urine collection 
Vacuum dry 

20 

0 

2 

14 

2 

Fecal collection 

20 

15 

3 

15 

3 

Combined urine/ 
fecal collection 

20 

15 

4 


5-4 















140 . 


160 . 


100 . 120 . 
TIME - HOUP<^ 


M. 


;ed Food Trays in Shuttle Appliances 






FOOD HEATING TRAYS HEAT FLOW BTU/HR 


D2-I1S57N 



40 . 60 . 80 . 100 . 120 . 140 . 160 . 


MISSION TIME - HOURS 

HEATERS Q-IN o LOSS TO AMBIENT 

U<IVAC not SC-A020 »*CE 1}, 

Figure 5-5. Food Trays Heat Loss in Shuttle Appliances Simulation 




1 


5-8 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 


DRY JOHN USAGE PHASE 









DRY JOHN DRYING PROCESS PARAMETERS 











D2-I18571-I 


6.0 SPACE STATION APPLIANCES SIMULATION 

To demonstrate that the new appliance subroutines are operational in an 
all-up G-189A system simulation, they were run in a Space Shuttle Orbiter 
and Modular Space Station G-189A ECLSS model. The appliances included in 
the models were taken from trade studies of many candidate concepts (Ref- 
erence l) and represent the optimum set of appliances for each vehicle. 

The results of these trade studies for a 180-day, six-man Modular Space 
Station are shown in Table 2-2. 

6.1 GENERAL SPACE STATION MODEL DESCRIPTION 

Since no operational model of the Space Station ECLSS was found prior to 
this appliance simulation effort, a simplified model of the pertitjent sub- 
systems was developed in which to check the appliance subroutines. The 
Modular Space Station concept (Reference 29) involves eight separate com- 
partments for crew habitability, work, experiments, etc. Only two of these 
compartments were included in the model. The first includes most of the 
personal hygiene equipment, and the other the crew eating and sleeping 
quarters. The G-189A flow schematic for these cabins and associated appli- 
ances is shown in Figure 6-1. The water loop used to supply these appliances 
is shown in Figure 6-2, Flow schematics for the other appliances not in- 
cluded in these two flow loops are shown in Figure 6-3. The following 
appliances were included in the Space Station model: 

- Refrigerators 

- Freezers 

- Food trays 

- Reverse osmosis unit 

- Shower 

- Clothes washer/dryer 

- Dishwasher/dryer 

- Dryjohn 

- Wet wipe wetting unit 


( 


J 



Figure 6-1. G-189A Flow Schematic of Space Station Cabin Gas Loop 
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Figure 6-2. G-189A Flow Schematic of Space Station Water Loop 
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6.1 (Continued) 

A description of each appliance is given in . the following paragraphs. The 
appliances were run in a 10-hour transient simulation according to a typi- 
cal daily schedule shown in Figure 6-4. The G-189A input model data and 
GPOLY subroutines are listed in Appendix F. 

6.2 REFRIGERATORS 

The CHILLR subroutine v/as used to model the Space Station refrigerator , 
component number 97, The input model data used for this component are 
included in Appendix F. It is identical to the refrigerator model dis- 
cussed in Paragraph 4.1.2. It was assumed the refrigerator could be 
connected into the Freon 21 radiator cooling loop. Since this loop was 
not required anywhere else in the model, a dummy input was made (as shown 
in Figure 6-3) to define constant Freon inlet flow and properties to the 
refrigerator. The appliance concept trade study (Reference 1) specified 
three refrigerators of equal size for Space Station. This was accomplished 
in a simplified manner for demonstration purposes by assuming all three 
refrigerator units were identical. The inlet Freon flow to the three 
refrigerators was multiplied by 1/3 immediately before solution of the 
component. Thus, only one unit was included directly in the flow loop. 

After component solution of the single unit, the outlet flow rate was then 
multiplied by 3 to include the total effect of all three refrigerator units. 
The heat loss from the refrigerators by door opening and conduction through 
the walls [in R(53)] was also multiplied by 3 and added to the total 
galley cabin sensible heat load (in R-array location 66) in the GPOLYl sub- 
routine. The locker door was assumed to be opened once every hour, as seen 
in the schedule of Figure 6-4. 

6.3 freezers 

The CHILLR subroutine was used to model the Space Station freezer, component 
number 62. The input model data used for this component are included in 
Appendix F. It is identical to the freezer model discussed in Paragraph 4. 1.4. 
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6.3 (Continued) 

A self-contained vapor compression refrigeration unit once under 
consideration for use on Shuttle (Reference 22) was assumed to provide the 
cooling. The freezer was included in the cabin gas loop, Figure 6-1, to. 
provide cabin air to cool the condenser. 

The appliance concept trade study (Reference 1) specified 13 freezers of 
equal size for Space Station. This was accomplished in the model in the 
same manner as was described in Paragraph 6*2 for the refrigerators; that 
is, all 13 units were assumed to have the same thermal performance, so a 
solution was actually obtained for only one. The air flow to the 13 freezer 
units was therefore divided by 13 in GPOLYl immediately before component 
solution and then multiplied again by 13 immediately after in GP0LY2. The 
heat loss from the freezer by door opening and conduction through the walls 
was also multiplied by 13 and added to the total galley cabin sensible 
heat load (in R-array location 66) in GPOLYl. Thus, the full effect of 
all 13 freezer units was included in the Space Station model. The freezer 
locker doors were assumed to be opened at random intervals during the day 
according to the schedule in Figure 6-4. 

6.4 FOOD HEATING/SERVING TRAYS 

The FTRAY subroutine was used to model the Space Station food heating trays, 
component number 95. The input model data for this concept are included in 
Appendix F. It is identical to the food tray model discussed in Paragraph 
4.2. The heating trays were not connected directly into any fluid flow 
loop, as indicated in Figure 6-3. The heat dissipated to the cabin from 
the trays (stored in the food tray R-array location 53) was added in GPOLYl 
to the total cabin heat load in galley cabin component R-array location 66. 
The travs are warmed for 2 hours just before the noon and evening meals, 
as seen in the schedule in Figure 6-4. It was assumed the trays for all 
six crewmembers were heated at the same time; thus, only a single component 
was used to simulate all six trays. It was assumed each tray had three 
heating cavities being warmed. Thus, the subroutine internally multiplies 
the heat loss for one heated cavity by 18 to account for all six trays. 
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6.5 REVERSE OSMOSIS UNIT 

The ROSMOS subroutine was used to model the Space Station reverse osmosis 
unit, component number 90, as seen in Figure 6-2. The input model data 
used for the component are included in Appendix F. The model design is 
identical to the unit described in Paragraph 4.3 except that it operated near 
ambient temperature, so no thermal conduction path was included between the 
unit and its surroundings. The product water was used as wash water for the 
clothes washer, dishwasher, shower, and wipe wetting unit. 


6.6 SHOWER 

The Space Station whole body shower interfaces with the cabin gas loop. 

Figure 6-1, and the water loop, Figure 6-2, for the required air and water 
flows. The SHOWER subroutine was used to model the shower stall, with 
the associated peripheral equipment (heaters, fan, pump, etc.) simulated 
with standard 6-189A component subroutines. The shower stall model data, 
component 3, are identical to the unit described in Paragraph 4.4. The 
flow schematic for the complete shower unit is shown in Figure 6-5. A 
description of the shower subsystem is given in Reference 19 and 24. The 
shower is operated six times in the morning for 10 minutes each, as shown 
in the schedule in Figure 6-4. The input data used to model each component 

are included in Appendix F. 

An electrical air heater v^as assumed (component number 19), with a 
thermostatically controlled temperature of 105°F. A heat exchanger was 
used to warm the inlet water. However, rather than model the heat exchanger 
a dummy component (number 22) was used with a fixed outlet water temperature 
of 105®F. A vortex water separator was used which does not require a drive 
motor. This component is treated using GP0LY2 logic (shown in Appendix F, 
GP0LY2 statement number 2000) which separates the inlet air and entrained 
water into two outlet flow streams: one gas, the other liquid. Since the 

air and water were mixed in the shower stall, no further mixing is required 
in this component. A fan and pump component (numbers 16 and 23) were used 
to circulate the air and water. 
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The normal order of solution for the shower loop components in Figure 6-5 is as 

i 

follows; 

22 15 16 -> 17 18 IP 3 24 ->• 20 -> 21 14 ->■ 23 

When the shower is not being used but the stall is at an average temperature 
above ambient due to previous usage, the normal solution path shown above is 
followed to permit thermal exchange between the stall and cabin environment. 

When the shower is not in use and the stall is in thermal equilibrium with 
the cabin, the solution of the loop is skipped over (solution path 22 23) 

to conserve computer time. 

When the shower subsystem is being used, the 90-second system computing time 
step was too large to achieve satisfactory accuracy of the shower loop mass 
balance. Instead, a 30-second interval was found to be adequate for solving 
the shower loop, and three successive solutions were obtained each system 
time step to achieve a total of 90 seconds. For each 30-second time step, 
an iterative solution was used to achieve a mass balance within the -loop. 

To begin the iterations, an estimate of the outlet flow conditions from 
component 20 is first made, after v/hich the loop is solved in the normal 
sequence given previously through component 20. The computed outlet condi- 
tions from component 20 are then compared with their estimated Values. If 
different, an updated estimate of the component 20 outlet flow is made, and the 
solution of the shower loop repeated, beginning with component 15. These 
iterations are continued until the estimate of the component 20 outlet flow 
conditions agrees with the actual computed values within 0.1 percent. The 
logic required to accomplish this iterative procedure is included in the 
GP0LY2 listing. Appendix F, under statement numbers 2000, 2200, and 2300 
(corresponding to components 20, 22 and 23 respectively). Usually only 3 or 4 
iterations are required to achieve convergence. These iterations assure that 
the overall mass balance within the shower subsystem (net mass inlet equal to 
net mass outlet) is valid. Since the shower is usually activated for only 
relatively short periods compared with other events (e.g., a maximum of six 
10-minute showers per day for Space Station), the required computer execution 
time for these added iterations is felt to be negligible. 
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6.7 CLOTHES WASHER/DRYER 

The Space Station clothes washer/dryer interfaces with the cabin gas loop, 
Figure 6-1, and the water loop. Figure 6-2, for the required air and water 
flows. The WASDRY subroutine was used to model the washer/dryer drum. 

The input model data for this component are identical to the washer/dryer 
model described in Paragraph 4.5. The associated peripheral equipment 
(heaters, fan, pump, etc.) were simulated with standard G-189A component 
subroutines. No washer/dryer unit has currently been built for space use. , 
However, preliminary design effort has been made, as described in Refer- 
ence 21, from which much of the Space Station washer/dryer design details 
have been taken. The flow schematic for the complete clothes washer/dryer 
unit is shown in Figure 6-6. The input data used to model each component 
are included in-Appendix F. The washer/dryer is operated twice per day, 
as shown in the schedule in Figure 6-4. 

Electrical heaters were used to warm the inlet air and water. The water 
heater (component number 29) was thermostatically controlled to 155°F. The 
air heater (component number 26) has one setting of 70.34 watts (240 Btu/ 
hr); it turns on when the dryer outlet air is below 135°F, and off when the 
air is above 155°F. A fan and pump (components 25 and 32) were used to 
circulate the air and water. Since the washer drum rotates, it was assumed 
the drum can be designed to act as a water separator: removing water at 

the outer periphery of the housing, and removing air near the center of 
tub rotation. Thus, no water separator was included in the flow schematic. 
An accumulator tank was used (component number 35) to store rinse water to 
subsequently be used as v/ash water. The tank dissipates heat to the cabin 
air through the tank walls. This heat is stored in R-array location 53 and 
is added in GPOLYl to the total cabin heat load in cabin component R-array 
location 66. Likewise, the clothes washer heat lost through its structure 
(in R-array location 53) is also added to the total cabin heat in R(66). 

The washer is operated in all seven cycle phases described in Paragraph 3.5, 
with 30 minutes allowed for washing and 30 minutes for rinsing. The auto- 
matic cycle phase switching logic, described in Paragraph 3.5, was used to 
control the washer/dryer operation. The drying phase was assumed to termi- 
nate when the residual water retention in the clothes material was 0.05 lb 
water/ lb clothes. 
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6.8 DISHWASHER/DRYER 

The Space Station dishwasher/dryer interfaces with the cabin gas loop. 

Figure 6~1, and the water loop. Figure 6-2, for the required air and water 
flows. The WASDRY subroutine was used to model the washer/dryer drum. The 
input model data for this component are identical to the washer/dryer 
model described in Paragraph 4.5. The associated peripheral equipment 
(heaters, fan, pump, etc.) were simulated with standard G-189A component 
subroutines. No washer/dryer unit has currently been built for space use. 
However, studies of various concepts have been made. Reference 27, from 
which some of the Space Station washer/dryer design details have been 
taken. The flow schematic for the complete dishwasher/dryer unit is shown 
in Figure 6-7. The input data used to model each component are included 
in Appendix F. The washer/dryer is operated three times per day. How- 
ever, only two of these times occur within the 10 hours of mission time 
being solved, as shown in the schedule in Figure 6-4. 

Electrical heaters were used to warm the inlet air and water. The water 
heater (component number 47) was thermostatically controlled to 155°F. 

The air heater (component number 40) has one setting of 70.34 watts (240 
Btu/hr); it turns on v/hen the dryer outlet air is below 135°F, and off 
when the air is above 155°F. A fan and pump (components 43 and 51) were 
used to circulate the air and water. Since the washer drum rotates, it 
was assumed the drum can be designed to act as a water separator; removing 
water at the outer periphery of the housing, and removing air near the 
center of tub rotation. Thus, no water separator was included in the flow 
schematic. An accumulator tank was used (component number 41) to store 
rinse water to subsequently be used as wash water. The tank dissipates 
heat to the cabin air through the tank walls. This heat is stored in 
R-array location 53 and is added in GPOLYl to the total cabin heat load 
in cabin component R-array location 66. Likewise, the dishwasher heat 
lost through its structure (in R-array location 53) is also added to the 
total cabin heat load in R(66). The washer is operated in all seven cycle 
phases described in Paragraph 3.5, with 30 minutes allowed for washing and 
30 minutes for rinsing. The automatic cycle phase switching logic, described 
in Paragraph 3.5, was used to control the dishwasher/dryer operation. 
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6.9 DRYJOHN 

The Space Station dryjohn interfaces with the cabin gas loop. Figure 6-1, 
and the water loop. Figure 6-2. The dryjohn flow schematic is shown in 
Figure 6-8. The WASTEC subroutine was used to model the urinal and commode 
collector, component number 121. The input model data for this componerat 
are identical to the dryjohn model described in Paragraph 4.6 -except for > the 
two minor exceptions noted in that section. The peripheral dryjohn 
equipment, shown in Figure 6-8, were simulated using standard 6-189A 
components. The G-189A input data used to model each component are 
included in Appendix F. 

The dryjohn was operated as a urinal and commode according to the schedule 
given in Figure 6-4. When not in operation, or when used as a urinal only, 
the’collector contents were assumed to be evacuated to space vacuum. A 
motor-driven centrifugal water separator, component 125, was operated at 
all times when the dryjohn was in use. The logic used to model this 
separator was included in GP0LY2, statement 12500, and is shown in Appendix 
F. This involved (1) diverting the entrained liquid to the component pri- 
mary side, and (2) adding the motor heat to the gas and liquid flows. The 
water from the separator was directed to the vapor compression-distillation 
unit, component 93, as shown in Figure 6-2. 

6.10 WET WIPE WETTING UNIT 

A wipe wetting unit was included ’in the Space Station model to provide 
hand washing and partial body washing. This concept, described in Refer- 
ence 30, dispenses water and soap into a washcloth, then rinses and wrings 
out the cloth after use. A motor-driven water separator is used with the 
component to provide inlet air flow and to separate the outlet gas and 
liquid. This device, component 79, is included in the cabin gas and water 
loops. Figures 6-1 and 6-2. 

The wipe wetting unit is operated intermittently during the day according 
to the schedule in Figure 6-4 taken from Reference 31. A single use of 
the unit is assumed to require one computing time step, which was 90 seconds 
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6.10 (Continued) 

During this time, the latent and sensible heat specified for this component 
in Reference 30 are added to the air and water flows. There is 0.00635 gm 
(0.014 lb) moisture added to the cabin air flow through the device and 
84 watts (286.6 Btu/hr) heat added from the separator motor. The air and 
water are assumed to leave the unit at the same temperature. The calcu- 
lations required to simulate this operation were included in GP0LY2, 
statement number 7900, and are given in Appendix F. 

6.11 SPACE STATION/APPLIANCES MODEL RESULTS 

The Modular Space Station G-189A mode! described in the preceding sections 
was run for a 10-hour transient simulation. The appliances included in the 
model ViSre operated according to the usage schedule shown in Figure 6-4. 

The results are plotted in Appendix G. These plots show cabin temperatures, 
humidity, CO^ levels, and the thermodynamic performance of the appliances. 

Ihe final solution for the run is also listed in Appendix G in the standard 
G-1S9A format. 

The total net effect of each appliance subsystem on the Space Station cabin 
air, coolant system, and water supply is shown in Figures 6-9 through 6-29. 

For the appliance subsystems requiring more than one G-189A component (clothes 
washer/dryer, dishwasher/dryer, shower, and dryjohn), the values plotted 
represent the combined effect from each component. The plotted sensible 
heat to cabin represents the heat transferred to the cabin air flowing 
through the unit as well as the structural heat leak from each component 
in the appliance subsystem. Some of the figures show heating or flow rates 
as short single-line "spikes" (e.g., Figures 6-9, 6-24, and 6-27). These 
values represent momentary component operation for one system time step, 
which was 90 seconds. 

The total appliance thermodynamic results shown in Figures 6-9 through 6-29 
are shown to illustrate typical appliance operation, and must be interpreted 
in the context of the specific assumptions used for the Space Station 
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6.11 (Continued) 

demonstration run. For example, the shower water usage and humidifying effect 
will vary depending on the specific water usage schedule used by the occupant, 
and the results shown in Figures 6-13 through 6-17 reflect the usage schedule 
in Figure 6-4. Since the washer/dryer hardware has not yet been designed, 
reasonable values were assumed for the thermal connections from the frame 
and accumulators to ambient. These values could vary widely with different 
structural designs, and this will affect the heat leak to the cabin. The 
thermal effect from the 13 food freezers could be exaggerated, since only 
one was modeled directly in the computer run, and the other 12 were assumed 
Identical. Therefore, the cooling units for all 13 freezers cycle on and 
off together, making a greater cabin heat input (Figure 6-11) than if their 
operation were staggered. 


The net thermodynamic performance of each appliance subsystem is demon- 
strated in Figures 6-9 through 6-29. The detailed thermodynamic parameters 
of interest in the individual appliance components are presented in 
Appendix 6. As discussed in Sections 6-2 through 6-10, the appliances 
included in the Space Station model were of the same design, with the same 
model input data, as the appliances described in Section 4 for verification 
of the models. Therefore, the Space Station appliance results plotted in 
Appendix G may be compared directly with the corresponding results in 
Section 4. The solutions for each appliance in the two sections are nearly 
identical, as expected, with all differences attributable to the difference 
in inlet and ambient conditions. Therefore, these results demonstrate that 
the new appliance subroutines are valid and operational in an all-up 6-189A 
system model. 
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Figure 6-17. Heat Input to Space Station Coolant from Shower 
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7.0 CONCLUSIONS AND RECOMMENDATIONS 

New G-189A subroutines have been developed to model future spacecraft crew 
appliances. These subroutines, some of which model more than one component, 
are generally described as follows: 

CHILLR (simulates a tnermally insulated locker cooled either by an 
externally chilled fluid or a self-contained refrigeration 
unit) 

* Freezer 

* Refrigerator 

FTRAY 

* Food warming/serving tray (Skylab-type) 

ROSMOS 

* Reverse osmosis waste water treatment unit 
SHOWER 

* Spacecraft whole body shower 
WASDRY 

* Clothes washer 

* Clothes dryer 

* Combined clothes washer/dryer 

* Dishwasher/dryer 

* Towel /cloth drying rack 

WASTEC 

* Dryjohn 

* Urinal 

The new appliance subroutines have been described, 

v/ith user s instructions and verification and demonstration results pre- 
sented for each. The appliance models have been shown to be valid and 
operational in an all-up G-189A system model. Most of the crew appliances 
described in this report have not yet been built; several have not been 
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7.0 (Continued) 

designed for zero-g application. Consequently, the verification of the 
accuracy of the models was limited largely to Independent analysis. It Is 
therefore recommended that the new appliance subroutines be closely 
correlated with future experimental data when testing of the actual 
appliance hardware Is done. 

Several of the crew appliances which have been modeled are relatively 
complex subsystems Involving a number of components such as pumps, fans, 
accumulator tanks, heaters, valves, etc. In such cases, new subroutines 
have been developed only for the Individual components which were not 
already available In the G-189A subroutine library. Consequently, a 
number of G-189A components may be required to simulate a single appliance 
(e.g., 12 components were required to model the clothes washer/ dryer) . 
Appliances such as these could be modeled and run by themselves using a 
typical range of ambient conditions such as temperature and humidity. The 
results from these runs could be converted Into tabular data which would 
be used In a new subroutine to simulate the performance of the entire 
appliance subsystem. This new subroutine would accept, as Input data, 
the actual concitlons existing during a G-189A run. It would then deter- 
mine from the tables the required performance parameters (I.e., temperature 
rise, humidity change, and heat dissipated to the cabin) for the appliance 
subsystem. This effort would allow much simpler modeling and faster com- 
puter execution of large ECLSS's involving many appliances. 


